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Measurements of induced noise currents along drifting cylindrical elec- 
tron beams have shown that noise fluctuations propagate as space-charge 
waves in the Same fashion as RF signals of the same fre quency. On man / 
such beams, however, the re gular standing-wave NOSE palte rn ws inte rrupted, 
afte r some drift distance, by a smooth steep increase in noise current, fol- 
lowed by slow, shallow undulations. This ‘growing noise’’ phenomenon, 
discovered by Smullin and his co-workers at M OI.T. SEVE ral years ago, is 
the subject of study in this paper. Its importance is considerable, in a nega- 
tive way, because it has hampered the development of medium-power travel- 
ing-wave-tube devices with acceptably low noise figures. 

The experimental measurements show the growing noise pattern to be the 
result of a two-stage process. Its primary cause is rippled-beam amplifica- 
tion of noise fluctuations over a wide band of microwave frequencies, much 
higher than the usual observation frequency. This explains its elusiveness. 
In the second stage, noise energy ts transferred to lower frequencies, due to 
intermodulation and other non-linear processes within the gain band. As 
the beat-fre quency norse increments are excited by continuous arrays of fre ’ 
quency pairs, their standing-wave patterns overlap one another, resulting in 


a smooth growing-noise pattern. 
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In Part IT of this paper, measurements of the noise spectrum of a rippled 
beam in the UHF region are described. These measurements reveal 
the presence of additional forms of instability. Calculations are made to 
account for some of these, and for aspects of rippled-beam amplification not 


pret tously understood. 


Part |— The Growing Noise Phenomenon* 


I INTRODUCTION 


When an REF probe is moved along a magnetically-focused electron 
beam in a drift region, the noise power is at first found to vary periodi- 
cally with distance from the electron gun.' For a sufficiently long beam, 
however, the periodic pattern is succeeded by an exponential rise, culmi- 
nating in an irregular plateau. This so-called ‘‘ growing noise’? phenome- 
non has been extensively investigated by its discoverers, L. Smullin and 
his colleagues at the M.I.T. Research Laboratory of Electronics.” 
They have established that this noise will begin to grow at a plane nearer 
the gun, and tend to grow at a faster rate, for electron beams (a) of 
higher perveance, (b) with less space-charge neutralization by positive 
ions, and (c) issuing from convergent, partly-shielded guns, rather than 
those immersed in the magnetic field. 

The growth of microwave noise power in drifting beams has hampered 
the development of high-power, traveling-wave tubes with acceptably 
low noise figures, as such devices generally have convergent, partly- 
shielded electron guns. The problem has been evaded in the design of 
low-noise, low-power traveling-wave tubes, by resort to confined-flow, 
parallel beams. 

Several theories have been proposed to explain the growing-noise 
wave: 

(1) Excitation of higher-order modes with complex propagation con- 
stants, by electrons threading the beam transversely ;* 

(2) Slipping-stream amplification, due to either longitudinal or trans- 
verse velocity gradients; 

(3) Rippled-beam amplification;” “'* and 

(4) Eleectron-electron interactions leading eventually to equipartition 
of thermal energy, and thus an increase in longitudinal velocity fluetua- 
tions. 

In Part I of this paper, measurements are presented which show that 
the principal cause of growing noise appears to be space-charge wave 


* Presented at the I.R.E. Electron Tube Research Conference, Boulder, Colo 
rado, June 27-29, 1956 
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amplification due to beam rippling. The mechanism is studied in some 


detail, as its connection with the usually-observed exponential rise of 


noise is not immediately apparent. In Part I], the UHF noise spectrum 
and its spatial distribution in beams with large-amplitude, long wave 
length ripples, are described. In addition, some of the underlying proc 
esses are analyzed 


If APPARATUS 


As sketched in Fig. 1, the heart of the apparatus consists of an electron 
gun, drift tube, and movable probe, all enclosed in a demountable, con 
tinuously-pumped vacuum system. Outside of the vacuum envelope 
there is a shielded solenoid, extending the entire 18-inch length of the 
drift tube. The annular gap between the solenoid pole face and the mag 
netic shield about the gun is nearly all taken up by a soft-steel section of 
the vacuum envelope. 

The electron gun is of the convergent Pierce type, with oxide-coated 
cathode and a coiled-coil filament heater producing negligible flux at 
the cathode surface. Surrounding the gun, and inside of the magneti: 
shield, is a small copper-wire coil that permits variation of this flux over 
a small range, either aiding or opposing the leakage flux due to the mai 
focusing solenoid. The flux density at the cathode has been approxi 
mately calibrated in terms of currents in both coils. Throughout the ex 
periments described below, the gun is pulsed with a 1,000 eps square 
wave of 2,200 volts on its anode, supplying 38 + 1 ma peak current in 
space-charge-limited emission. 

The novel feature of the probe is that its annular RF pickup gap 
couples to a 50-ohm coaxial line leading to the receiver, rather than to a 
resonant cavity. This permits RIF power measurements over a wide 
range of frequencies. The inner conductor of the coaxial line serves as 
current-collector, being isolated and biased positively about 40 volts 
with respect to the outer conductor to prevent escape of slow secondaries 
An adjustable vane can be locked in position in front of the probe (whose 
entrance aperture is 0.100 inch in diameter), so that circular apertures 
of various smaller sizes are fixed on the probe centerline, about 0.070 
inch in front of the probe. With these apertures, measurements of col 
leetor-current variations along the beam furnish a rough picture of 
beam-ripple amplitudes and locations. In addition, the current-density 
variation across the beam can be estimated by moving a pinhole aperture 
in a broad are through the beam centerline. Both the inner conductor 
of the probe and the intercepting vane are liquid-cooled. 


The noise powers coupled to the coaxial probe are considerably smallet 
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than for a tuned coaxial cavity, because of the lower RF gap impedance 
of the former. To compensate for this drawback, a sensitive noise re- 
ceiver is employed, similar in principle to the radiometer invented by 
R. H. Dicke.” The input noise power is replaced periodically by a 
matched load at room temperature by pulsing the beam on and off 
with a 1,000 cps square wave, and placing an isolator in front of the re- 
ceiver. A synchronous detector eliminates gain-fluctuation noise and 
converts the receiver output to a de voltage. 

Noise power variations at various microwave frequencies are measured 
in terms of the changes of attenuation, between probe and receiver, re- 
quired to keep the receiver output constant. These rapid adjustments in 
attenuation are performed by a servo amplifier-motor loop, and recorded 
on a chart, whose speed (13 inches per minute) is synchronized with that 
of the moving probe. In the same way, records of collector current as a 
function of probe position can be obtained, and correlated with those of 
noise power. The probe can be moved a distance of about 17 inches, its 
position nearest the gun (2 0 inches) corresponding to a distance of 


0.95 inch between the anode and the input plane of the RF gap. 
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Fra. | Cross-section of experimental tube, showing electron gun, probe, and 
two solenoids. The isolated current-collector electrode serves as inner conductor 
of a coaxial line. The induced RF power can be measured over a wide range of 
frequencies 
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Iii EXPLORATORY MEASUREMENTS 


The electron gun used in these experiments had been designed for use 
in a helix traveling-wave tube with a longitudinal focusing field of 600 
gauss. Noise-power and collector-current curves, therefore, were first 
taken with 600 gauss to study a typical state of affairs in an operational 
beam. As seen in Fig. 2, the noise power at 3.9 kme varies periodically 
with distance for about 4 inches from the gun, then climbs rather 


smoothly by nearly 23 db to an irregular plateau, where it undulates 


slowly, and finally levels off. The initial part of the growing noise curve 
at 10.7 kme is missing because of inadequate receiver sensitivity, but its 
later portion is similar to that at 3.9 kme, with about half the rate of 
noise climb. With the 0.020-inch aperture, the collector-current varia- 
tions decrease in amplitude chiefly in the drift region preceding the noise 
climb; whereas those for the 0.100-inch aperture decrease afterwards 
Both curves show a flattening in the growing-noise region itself, as well 
as a decrease in their average values after that region, signifying an in 
crease in the average beam diameter. 

A similar set of curves is shown in Fig. 3, for a focusing field of 279 
gauss (about twice the nominal Brillouin field). Noise growth at 3.9 kme 
starts later, and proceeds less steeply, than at 600 gauss. The noise-power 
curve for 10.7 kme is much more articulated, with a semblance of peri- 
odicity, throughout the drift region. Collector-current curves for both 
0.020- and 0.050-inch apertures show considerable reduction in current 
ripple amplitude with distance, reaching virtually zero in the former case 

Another type of survey measurement is illustrated in Fig. 4. With 
the probe stationary at the far end of the drift space (about 18 inches 
from the gun anode), the main solenoid current is varied smoothly to 
change the focusing field from 0 to over 600 gauss, and synchronized 
records are made of collector current and noise power. (In this instance, 
the current in the auxiliary solenoid was +3.2 amperes.) At low mag 
netic fields, both the current and noise-power curves have large ampli- 
tude variations, which diminish as the field increase. At first glance, the 
noise peaks and valleys seem to coincide with those of collector current ; 
certainly, some do. Closer inspection, however, reveals significant mis- 
alignments which cannot be accounted for by experimental error. When 
the three noise curves, at 3,050, 3,930, and 4,730 me, respectively, are 
compared with each other, some characteristic features emerge: 

(1) An average curve drawn through each pattern has one or two 
broad maxima, which tend to move toward higher field strengths with 
increasing frequency. 
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Fig. 2 Typical smooth steep growing-noise patterns, near 4 and 10.7 kme 
respectively, with customary focusing field of about four times the Brillouin value 
Collector-current traces through small and large apertures reveal decreases in 
ripple amplitude and increase in average beam diameter 
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Fig. 3 At about half the focusing fields used in Fig. 2, the growing-noise 
pattern is much the same at 4 kme, but shows significant articulation at 10.7 
kme. The collector-current traces show pronounced decreases in ripple amplitude 


and differences in ripple patterns obtained with different apertures 
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Fic. 4 — With the probe stationary at about 18 inches from gun anode, collec 
tor current and noise powers in bands about 1 kme apart are recorded, as current 
in the main solenoid is varied from 0 to 1.5 amperes (0 to 635 gauss). 
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(2) The lower the frequency, the lower the field strengths at which 
noise amplitudes change most violently with field. 

(3) The three noise curves resemble each other in small details. 

The results of a great many records of the kind illustrated by Figs 
2 and 3 can be summarized as follows: 

(1) There is always a decrease in beam-ripple amplitude associated 


with noise growth at any frequency. (Sometimes the ripple amplitude 


increases afterwards, as in Fig. 3.) 

(2) The higher the frequency, for a given field, the more articulated 
or scalloped the noise pattern. 

(3) No correlation can be found between rate of noise growth and 
either (a) distance from gun to take-off plane, or (b) net gain at the end 
of the drift region. The trends, as a function of magnetic field, are differ- 
ent at different frequencies. 

(4) Greatest noise growth does not, as a rule, occur with zero flux 
threading the cathode. Sometimes two nearly equal peaks occur for two 
values of B., each of opposite polarity, referred to the sense of the main 
field. 

(5) The noise-distance patterns change very slowly with frequency. 

(6) No beam entirely ripple-free throughout its length has ever been 
observed by the writer. 


IV ORIGIN OF GROWING NOISE 


If noise growth is due to some amplification process, it should be 
possible to adjust the beam-focusing conditions so that the noise currents 
start increasing at the anode, and attain the greatest possible over-all 
gain at the end of the drift space. The enhanced activity of the unknown 
gain mechanism should presumably help identify it. The curves of Fig. 
t show that maximum noise occurs at different values of the focusing 
field, for different values of field at the cathode, and different probe posi- 
tions. With the anode voltage and receiver frequency fixed, therefore, 
the conditions for greatest net noise growth can only be found by a series 
of trial settings of both magnetic fields, each followed by a recording of 
the noise-distance pattern. Eventually, a set of fields can be found for 
which the greatest total gain occurs; and such patterns are usually found 
to show fairly steady noise-amplitude increase, on the average, over the 
entire length of probe travel. 

The results of this procedure for noise power near 4 kme, as well as 
the patterns of collector-current versus distance with the same fields, 
using the 0.100-, 0.050-, and 0.020-inch apertures fixed at the probe 
centerline, are shown in Fig. 5. A similar set of records, for noise power 
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Fig. 5 The magnetie fields in drift tube and at the cathode have been adjusted 
empirically to expand the growing-noise region over the entire drift region, with 
the L.O. at 3,990 me. The field is slightly less than the Brillouin value, but the 
beam is strongly rippled because the gun was designed for best focusing at a much 
higher field. The noise-current maxima align with the average collector currents 
on their increasing slopes, for all three aperture sizes. 
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near 10.69 kme, is shown in Iig. 6. The significant features of both sets 
of records can be summarized as follows: 

(1) In both cases, the beam-ripple periods are equal to the RF scallop 
periods; i.e., the half-wavelengths of the space-charge standing waves. 
The noise minima tend to occur at planes where the collector currents 
are at their average values and decreasing; i.e., where the beam diame 
ters are at their average values and increasing. The noise-current maxima 
occur where the beam diameters are about to decrease. These are the 
classical conditions for rippled-beam amplification.” 

(2) In Fig. 5, the ripple amplitudes and peak values of all three col 
lector-current curves decline appreciably with distance, the rate of de 
cline being greatest for the smallest aperture. (Similar curves, not shown 
here, have displayed little or no such decline in the absence of noise 
growth.) This suggests that the RIF noise power is amplified at the ex 
pense of de energy associated with radially-directed electron velocities 

(3) In Fig. 6, the disparity among rates of decline of current-ripple 
amplitudes and their peak values, for the three aperture sizes, is even 
more pronounced. In addition, the ripple wavelength barely changes 
for the 0.100-inch aperture, but increases with drift distance for the 
smaller apertures, resulting in an increasing ‘‘ phase shift’? among them 
Thus the current-density variations at different radi in the beam can 
contribute unequally to space-charge wave amplification, depending on 
their local ripple amplitude and phase. In this instance, the variations 
in current density along the beam are initially greatest near the axis, 
and suffer the greatest reduction there. It is worth noting that this 
‘inner rippling’ would be missed entirely in beam-size measurements 
with a large aperture.* 

The decrease of beam ripple and the increase in average beam diame 


ter, shown in Figs. 5 and 6, has been found to accompany rippled-beam 


amplification of impressed signals by T. G. Mihran.” Another corrobora 


tion of the identity of this gain mechanism can be obtained by comparing 
the measured noise gain per scallop with that predicted by theory for 


idealized conditions. lor a beam with stepwise alternations of maxi 
mum and minimum beam diameters (ratio ro/r;), and with noise maxima 


and beam-diameter maxima coinciding, the gain per scallop is as follows: 


G. ) e (2) 
V2 ry Do 
Here, V is the beam potential, and p the reduction factor w,/w,. Al 


though the actual rippled beam is far removed from either Brillouin or 


* More information about ‘‘inner rippling”’ will be presented in Part [I 
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Fic. 6 The fields have been adjusted as in Fig. 5, for maximum extension of 
the gain region, for noise power near 10,690 me. Collector current measured with 
the smallest aperture shows the greatest decline in amplitude of variations, as 
well as advance in ripple phase relative to the current through the largest aperture. 
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TABLE MEASURED VERSUS CALCULATED MaxiIMuM NoIse GAIN 


Ripple Data Gain in db per Scallop 
Freq. m« 


Iris Dia. Inches Brillouin Flow Confined Flow Measured 


().020 
3,990 0.050 
0.100 
0.020 
10,690 0.050 
0.100 


confined flow, the published values of reduction factor for both extremes 


~ The ratio ro/r; can be estimated 


. . . If 
can be used as first approximations. 


by assuming the current density to be uniform over the beam cross- 
section near the middle of the drift region, for each of the three apertures 
used. The potential variations can be neglected. The results of such cal- 
culations are given in Table I. 

As the computed gains are expected to be somewhat greater than those 
measured, because of the optimum conditions assumed, the best cor- 
respondence between measured and computed gain rates appears to be 
for the ripple data taken with the 0.050-inch iris at 3,990 me, and that 
with the 0.020-inch iris at 10,690 me. This distinction is in accord with 
previous qualitative comparison of Figs. 5 and 6, showing that most of 
the beam cooperates in the ripples of the former, but that “inner rip- 
pling” characterizes the latter. 

Another calculation that reveals which part of the beam is interacting 
with the RF noise field in each case is that of the space-charge half- 


wavelength, as follows: 


rb 
p- 8,5 : 
where 


8, 1747'""/V"". 


Here 8, is the plasma wave number, ) the beam radius, and p the reduc- 
tion factor, which can be evaluated as previously for the smooth beam 
in either ideal Brillouin or confined flow. For the gun used here, the 


square root of the perveance is 


I’? /v** = 0.606 <X 10° MKS units, 


\,/2 & 29.8 b/p. 
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TasBLe II — MeasurReED VERSUS CALCULATED SPACE-CHARGE 
HALF-WAVELENGTHS 


Ripple 
Wavelength 
Brillouin ‘onfined 4, Meas ¢ 
Flow Flow 
3,990 0.050 0.057 
10,690 0.020 0.033 


Thus, agreement between this expression and the measured value 
requires the correct choice of the effective beam radius, b. It turns out 
that the suitable value for Fig. 5 (3,990 me) is the average beam radius 
obtained from ripple data taken with the 0.050-inch iris, and that for 
Fig. 6 (10,690 me) is obtained with data taken with the 0.020-inch iris. 
The results are summarized in Table IT. 

With this mechanism as the primary source of the noise gain, it be- 
comes clear why nearly equal noise maxima were found, with some 
values of the main focusing field, B, for two values of cathode flux den- 
sity B. of opposite polarity. from approximate analyses of beam ripples 


when flux threads the cathode, such as those provided by McDowell’ 


and others, it is found that the ripple wavelength depends nearly alto- 
gether on B. Its amplitude and spatial phase, however, depend on B, , 
as this affects the beam geometry at the drift-space entrance. For a 
sufficiently wide range of variation of B, , the spatial phase of the ripples 
can be varied from the proper relation with the space-charge standing 
wave for gain, through the positions for de-amplification, and back to 


gain again. 


V THE GROWING-NOISE MECHANISM 


Although many earlier noise records can be understood in the light of 
the rippled-beam amplification (RBA) process, this is not yet true of 
the smooth, steep noise growth usually observed, as in ligs. 2 and 3. 
The simple theory predicts that a space-charge wave will be amplified 
when, for small ripples, a ‘‘resonance’’ condition exists between the rip- 


ple wavelength, L, and the space-charge half-wavelength 
L = nv,/2, 


where n is an integer, usually unity. In addition, as mentioned earlier, 
there is an optimum phase relation between ripple and standing wave 
for maximum gain. These conditions are not satisfied by the records of 





PART |} THE GROWING NOISE PHENOMENON $45 


ligs. 2 and 3, except possibly for the 10.7 kme noise current in Fig 
(at a relatively low magnetic field). 

To establish a connection between the two types of noise growth, the 
noise record of Fig. 6 (for 10.7 kme with greatly expanded gain region 
is compared with that near 4 kme under the same conditions, in Fig. 7. 
The growing noise region for 4 kme does not start until at least four seal- 
lop wavelengths past the earliest observed 10.7 kme noise growth. More- 
over, the 4-kme noise pattern resembles that for 10.7 kme in many de- 
tails. (The resemblance in details of noise patterns at nearby frequencies 
has been remarked before, in connection with Fig. 4.) 





50 
NOISE POWER 
45 NEAR 
10,690 MC 
40 ‘=e T — T 

. 4 
14 











NOISE POWER 
AT 
3990 MC 
| as a 


+ 


” 
= 
ud 
@ 
VU 
Ww 
a 
Zz 
x 
Ww 
= 
O 
a. 
WwW 
” 
12) 
= 
Ww 
> 
- 
< 
| 
WwW 
x 











4 6 a 10 
Z IN INCHES ALONG DRIFT TUBE (0 IS 0.95” FROM ANODE) 
B= 260 GAUSS, Bc= 31.2 GAUSS, V=2200 VOLTS, I=38Mé 


Kia. 7 The pattern of growing noise in Fig. 6, near 10,690 mes, is compared 
with that near 3990 me for the same fields. The gain region of the latter curve 
starts much later and is much smoother than the former. The small irregularities 
on the 3,990 me curve resemble the scallops of the 10,690 me curve, in a blurred 
way 
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This information suggests that the growing-noise mechanism is really 
a two-stage process: amplification of a broad band of microwave fre- 
quencies, located far above the observation frequency, followed by a 
transfer of noise energy to lower frequencies. 


(a) First stage 


At low magnetic fields, there are few ripples per unit length, but their 
amplitude is usually large; whereas at moderate to large magnetic fields, 
the ripple amplitude is small, but so is the ripple wavelength. In either 
case, the bandwidth of RBA is large, usually many thousands of mega- 
cycles. 

The increase of both bandwidth and gain per scallop with ripple 
amplitude has been explained by Pierce’ by analogy between the gain 
band of a rippled beam and the stop band of a transmission line filter: a 
sharply varying periodic disturbance on the latter will reflect short as 
well as long waves, whereas smoother perturbations will not reflect the 
shorter waves to any extent. 

Another way to study the amplification bandwidth is to derive the 
equations for RF current in a one-dimensional beam with sinusoidal varia- 
tion of the reduced plasma wave number, 8, = p:8,, as Heffner," 


15 16, 17 rm: . ° 
Bloom,” and others have done. This leads to a Mathieu equation, 


whose solutions may be studied on the Mathieu stability plot (A, q): 


d'I 


2 + (A — 29q cos 2x)] ; (6) 
dx? 


Here J is the RF current, g a measure of the perturbation amplitude, 
x = mz/L, and A = (2L/),), where L is the ripple wavelength and X, 
the reduced plasma wavelength. Bloom has shown that, if n is the in- 
tegral number of scallop wavelengths between initial and final planes, 
the greater the product ng, the greater the total amplification or deam- 
plification, and the less critical the phase relation between RF standing 
wave and ripple for amplification. Ultimately, for very large ng, ampli- 
fication will take place for all values of this phase angle. 

At higher magnetic fields, both the ripple amplitude and ripple wave- 
length are decreased. This means that q is reduced, but n increased over 
any fixed span. This combination usually tends to increase the product 
ng up to some fairly high field, after which it may decline. More impor- 
tant, the reduction in ripple wavelength shifts the band of amplification 
to higher frequencies, and greatly increases the frequency band. This 
occurs because the ‘“‘resonant’”’ space-charge wavelength is shorter, 
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and short space-charge wavelengths correspond to high frequencies, 


where the former change very slowly with frequency. 
(b) Second staqe 


When noise power over this large band has been amplified sufficiently, 
electron bunching becomes non-sinusoidal, and the beam becomes non 


linear. Harmonics and beat-frequencies” of the fundamentals, and pos- 


. ° 19 : . ° 
sibly sub-harmonics, are excited. As the beat-frequencies are excited 


by a continuum of pairs of frequencies, their standing-wave patterns 
overlap one another, resulting in a ‘‘wash-out”’ of the noise minima, and 
a smooth growing-noise pattern. Eventually, the same non-linear proc- 
esses take place within this subsidiary band, leading to a gradual leveling 
of the entire noise spectrum. The initial rates of rise of the intermodula- 
tion products, however, should take place closer to the gun and be 
greater, for a lower frequency. They will depend on both the spectrum 
of noise power in the primary band and, so to speak, the spectrum of 
“beam non-linearity” within that band. 

To simulate this intermodulation process, two low-level klystron sig- 
nals (9,050 and 12,275 mc, respectively) were simultaneously permitted 
to modulate the electron beam as it entered the drift tube, by means of 
a short length of lossy helix. The magnetic fields at the cathode and in 
the drift space were adjusted to produce a beam ripple which amplified 
both of these signals simultaneously over most of the drift space, as 
shown in Fig. 8 (a, b, ¢).* Noise-power records were then made at the 
difference-frequency, in the presence of the two modulation signals, 
lig. 8(d), and in their absence, Fig. 8(e). The difference between the 
noise levels in the latter two records increases with distance, as both 
parent space-charge waves grow in amplitude, and the degree of beam 
non-linearity increases. Naturally, the contribution to 3,295-me noise 
in the absence of modulating signals is far greater than that of the latter 
two alone, as the primary bandwidth of noise amplified is very great, 
and that of the signals very small. 

There are several reasons why exponentially-growing noise should 
stop growing and level off, and sometimes even decrease slightly: 

(1) depletion of de kinetic energy in the beam ripples; 

(2) de-amplification in the fundamental band, due to departure from 
the proper phase relation for gain between standing wave and ripple, 
if only over part of the band (Fig. 3); 

* The fine ac detail superimposed on the pattern of Fig. 8(b) is due to inter 


ference between the waves traveling along the beam and that propagating as a 
waveguide mode in the drift space 
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Fic. 8 (a), (b), (ce) The fields have been adjusted to give rippled-beam 
amplification of two weak klystron signals, 12,275 and 9,050 mc, simultaneously 


impressed on the beam at the entrance to the drift region. 
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Fic. 8 (d), (e Noise power at the difference-frequency, 3,225 me, is recorded, 
with and without the signals (b) and (c) present. The difference in ordinates of 
the two curves increases with distance from the gun, as the impressed signals 


grow and the beam becomes more non-linear 


(3) sufficient phase shift between inner and outer ripples in the beam 


for one to de-amplify as much as the other amplifies; and 


(4) interference among the intermodulation products excited at differ 
ent positions along the beam.* 
The last effect is illustrated in Fig. 9, in exaggerated form. The beam 


° Suggested by €>...6; Cutler of Bell Telephone Laboratories 
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Fic. 9 Two strong klystron signals are impressed on the beam as in Fig. &, 
and noise power at their difference-frequency recorded, both with and without 
these signals present. The deep noise minima in (b) are due to destructive inter 
ference between trains of waves excited by intermodulation at different positions 
along the beam. 
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is Simultaneously modulated as before with two klystron signals (8,400 
and 11, 590 me, respectively), but now at fairly high level; and the fo- 
cusing field is made large. The interference dips in the pattern of 3,120-me 
noise are quite deep, and are spaced irregularly and farther apart than 
the space-charge wavelength of any of the three frequencies involved 
The third dip is shallower than the previous two because of the growth 
of 3,120-me noise other than that due to the signals, as shown in Fig. 
%(c). The latter pattern of noise in the absence of the two high-level, 
high-frequency signals suggests that the characteristic first gentle dip 
following the growing-noise region is indeed of the same nature as the 
artificially-produced interference dips, and has nearly the same quasi 
period. 

The pattern of dips agrees with simple calculations, based on this 
model, in which the amplitude of the difference-frequency intermodula 
tion product, excited at any plane ¢, is assumed proportional to the prod- 
uct of the amplitudes of the two high-frequency space-charge standing 


waves, as follows: 


where 
[, SIN PrBpf- SIN w °/u), n 


The total current at ¢ zis the sum of contributions from all the stand 
ing waves excited to the left of it: 


4 
phils COS P38 p(Z “)[COS (pi — Po)Byt 


‘OS (Pi + P2)P pl d¢ 


This expression is readily integrated and evaluated 


VI CONCLUSIONS 


Synchronized measurements of electron-current density and noise 
currents at several microwave frequencies have shown that the “growing 
noise’ pattern in drifting cylindrical beams is the result of a two-stage 
process. In the first stage, rippled-beam amplification of noise fluctua- 
tions takes place over a very broad band of microwave frequencies, 
much higher than the usual observation frequency. In the second, noise 
energy is transferred to lower frequencies by intermodulation and other 
non-linear processes within this band. The element of non-linearity is 
supplied when primary noise gain is sufficient to make electron bunching 
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non-sinusoidal. Other sources of non-linearity are thermal velocities, 
non-laminar beam flow, etc. As the beat-frequency noise increments at 
any plane are produced by continuous arrays of frequency pairs, in- 


creasing in numbers and amplitude in various ways as primary amplifica- 


tion proceeds, the multiple standing-wave patterns at the observation 
frequency progressively overlap one another. This results in the smooth 
steep rise of noise power usually observed. 

Phase correlation among the space-charge waves excited at succes- 
sive planes on the beam by the same set of frequency pairs is indicated 
by gentle dips, due to their destructive interference, in the plateau 
following the initial noise rise. 

tippled-beam amplification occurs whenever the ripple wavelength 
and half the space-charge wavelength are nearly equal, and bear a 
favorable spatial relation to each other. However, this ‘“‘phase’’ relation 
becomes less critical with an increase in either the number of ripple 
wavelengths over which synchronism persists, or the ripple amplitude, 
or both. Noise amplification by this mechanism, therefore, is probably 
present to some degree in all rippled streams, particularly at high fields. 
The extreme difficulty encountered in focusing ripple-free beams from 
convergent, shielded guns has to this date prevented the detection of 
any other primary gain mechanism, which may conceivably co-exist in 
such beams. 

A conspicuous feature of rippled-beam amplification is the decrease 
in ripple amplitude due to conversion of de into ac kinetic energy. Such 
changes in beam structure emphasize the inadequacy of beam-flow com 
putations based entirely on de force equations. A more detailed descrip- 


tion of this de-ac energy conversion is given in Part IT. 


ACKNOWLEDGMENTS 


The experimental apparatus could not have been built without the 
combined efforts of many associates of the writer, principally A. R. 
Strnad, P. Hannes, J. 8. Hasiak and J. M. Dziedzic. The author is also 
indebted to R. Kompfner, C. F. Hempstead and K. M. Poole for valuable 
suggestions; and above all to C. I°. Quate for constant encouragement 


and advice 


REFERENCES 


1. C. C. Cutler and C. F. Quate, Experimental Verification of Space-Charge and 
Transit Time Reduction of Noise in Electron Beams, Phys. Rev., 80, p. 875, 
1950 

D. Smullin and C. Fried, Microwave Noise Measurements on Electron 
Beams, Trans. I.R.k. ED-1, No. 4, p. 168, Dee., 1954. 





2 


0 


t. 


0. 


6. 


é. 


PART I THE GROWING NOISE PHENOMENON 853 


C. Fried, Noise in Electron Beams, Tech. Rep. 294, Research Laboratory of 
Electronics, M.1.T., May 2, 1955. 

J.R. Pierce and L. R. Walker, Growing Waves Due to Transverse Velocities, 
B.S.T.J., 35, p. 109, Jan., 1956. 

G. G. Maefarlane and H. G. Hay, Wave Propagation in a Slipping Stream of 
Electrons: Small Amplitude Theory, Proc. Royal Soc. (B) 63, p. 409, 1950 

C. K. Birdsall, Rippled Wall and Rippled Stream Amplifiers, Proc. I.R.1 
42, p. 1628, Nov., 1954. 

R. W. Peter, S. Bloom, and J. A. Ruetz, Space-Charge-Wave Amplification 
along an Electron Beam by Periodic Change of the Beam Impedance, RCA 
Rev., 15, p. 113, March, 1954. 


. T. G. Mihran, Sealloped Beam Amplification, Trans. I.R.E., ED-3, No. 1, p 


32, Jan., 1956. 


. R. H. Dicke, The Measurement of Thermal Radiation at Microwave Frequen 


cies, Rev. Sci. Instr. 17, p. 268, July, 1946 


.W. W. Rigrod and J. A. Lewis, Wave Propagation Along a Magnetically 


Focused Cylindrical Electron Beam, B.S.T.J., 33, p. 399, March, 1954 


. G. M. Branch and T. G. Mihran, Plasma Frequency Reduction Factors in 


Electron Beams, I.R.E. Trans., ED-2, No. 2, 3, April, 1955. 

Informal communication from H. L. MeDowell 

Informal communication from J. R. Pierce 

Informal communication from H. Heffner 

8S. Bloom, Space-Charge Waves in a Drifting, Scalloped Beam, unpublished 
RCA Research Laboratories report 


». O. E. H. Rydbeck and B. Agdur, Propagation of Space-Charge Waves in 


Guides and Tubes with Periodic Structure, L’Onde Electrique, 34, p. 499, 
June, 1954. 

P. V. Bliokh and Y. B. Feinberg, Space-Charge Waves in Electron Beams 
with Variable Velocity, Zhurnal Tekhn. Fiziki, 26, p. 530, March, 1956 


_C. C. Cutler. The Nature of Power Saturation in Traveling Wave Tubes, 


B.S.T.J., 36, p. 841, July, 1956 


S. Lundquist, Subharmonie Oscillations in a Nonlinear System with Positive 


Damping, Quarterly of Appl. Math., 13, No. 3, p. 305, Ocet., 1955 








Noise Spectrum of Electron Beam in 


Longitudal Magnetic Field 
Part Il — The UHF Noise Spectrum 


By W. W. Rigrod 


Manuscript received January 21, 1957 


Sharp peaks are found in the UHF spectrum (10 to 500 me) of an elec- 
tron beam, emanating from a shielded diode. In the presence of a longitudi- 
nal magnetic field, the strongly rippled beam displays an additional set of 
peaks whose frequencies are proportional to the field strength. The largest 
of these, just above the cyclotron frequency, is connected with the overlap of 
a dense cluster of particle orbits, passing close to the beam axis. It can attain 
amplitudes of 65 db above background noise. 

The transverse distribution of UHF noise power is found to agree with 
that for ideal Brillouin flow, even in rippled beams. With long ripple wave- 
lengths, two noise maxima are found to flank each beam waist. A small- 
signal wave analysis explains this pattern, and affords some insight into thi 
energy-exchange processes in rippled-beam amplification. The reduction in 
“growing noise” due to positive tons is attributed to increased cancellation 


of net radial beam motion, due to overlap in particle orbits near the axis 


[I INTRODUCTION 


The reader is referred to Part I’ for a description of the experimental 
apparatus and its operation. In this paper, measurements of noise power 
in the same electron beam are described, with frequencies chiefly in the 
10- to 500-me range, and relatively weak magnetic fields. For the UHF 
measurements, a calibrated coaxial step attenuator and a super-regenera- 
tive receiver (the Hewlett-Packard 417-A VHF Detector) are used. Rela- 
tive noise-power amplitudes at fixed frequencies are measured as before, 
in terms of changes in attenuation between probe and receiver required 
to restore constant receiver output. To obtain qualitative information, 
however, such as the location of noise maxima along the beam, the series 
attenuation is fixed. The receiver output is amplified, rectified, and per- 


855 





856 THE BELL SYSTEM TECHNICAL JOURNAL, JULY 1957 


mitted to register itself directly on the chart recorder, whose motion is 
synchronized with that of the probe. Very roughly, the detector output 
varies as the log of input power. 

Measurements are described (a) of the UHF noise spectrum in the 
beam, just outside the gun anode; (b) of this spectrum at the end of the 
drift region, in a longitudinal magnetic field; (¢) of the noise-power dis- 
tribution along the axis; and (d) transverse to the axis of the rippled 
beam in the drift region. Two calculations are then outlined, one of wave 
propagation along the rippled beam (to explain the observed distribution 
patterns), and the other to account for some spectacular peaks in the 
beam spectrum (b). 


Il FIELD-INDEPENDENT PEAKS 


When the noise spectrum of an electron beam is scanned by a tunable 
receiver, it is found that an irregular array of narrow-band peaks char- 
acterize the UHF region, below about 1000 me. Of these peaks, some are 


due to spurious modulation effects, and can be eliminated as follows: 


(1) Transit-time oscillations due to positive ions, secondary electrons, 
or both. Such frequencies vary with probe (collector) position. 

(2) Resonances in the probe and receiver, excited by the pulsed- 
voltage supply. These are unaffected by changes in collector current. 

(3) Ion oscillations in the electron gun or beam. Their frequencies vary 
with anode voltage. 

The remaining narrow-band peaks fall into two classes, depending 
on whether their frequencies vary with the magnetic field. 

Well-defined peaks can be detected with the RF probe stationed one 
inch from the gun anode, with or without any focusing field. When the 
beam is focused by a longitudinal magnetic field, these disturbances 
propagate along the beam, and tend to increase in amplitude with dis- 
tance, but not to change in frequency. A typical set of such frequencies, 
within the range of the tunable receiver is as follows: 15.9, 24.8, 31.2, 
34.0, 48.5, 63.4, 77.0, 108, 151, 166, 270.5, 372 and 481 me. (During this 
measurement, the anode voltage was 2,200, and the peak current about 
10 ma.) 

No consistent relation could be found between these frequencies and 
either the anode voltage or the cathode temperature, although unmis- 
takable frequency changes did occur when these parameters were ma- 
nipulated. Failure to establish such a relation may have been due to 
uncontrolled drift in cathode activity. In any case, the measurements 
did serve to narrow the field of possible mechanisms, by eliminating the 


following: 
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(1) Transverse positive-ion oscillations, for which the frequencies 
vary as the square root of anode voltage. 

(2) Transverse electron plasma oscillations (near or beyond the 
anode), for which the frequencies would be too high. 

(3) Longitudinal electron plasma oscillations at the potential mini- 
mum, for the same reason (should be near 2,500 me). 

(4) Longitudinal diode oscillations.’ When the electron transit angle 
through the diode is approximately (n + 1) periods, where n is an in- 
teger, the real part of the diode conductance becomes negative, permit- 
ting oscillations to occur. Again the frequencies of such oscillations would 
be too high, (2,200 me and higher) for the gun used, to conform to the 
observed values. 

There is, however, one published theory for which an order-of-magni 
tude correspondence does exist between the measured and calculated 
frequencies. Klemperer” ° has shown that a strip beam tends to break 
up into clusters of ‘‘ pencils” at the cathode. He ascribes these to standing 
waves resulting from transverse oscillations in the space-charge cloud, 
and offers an expression for the wave velocity in this medium. Applica 
tion of his formula to the cathode used in the present experiments results 


in a least frequency of 31.3 me. Other observers, such as Smyth’ and 


Veith,’ have also reported evidence of interaction between electrons in 
a retarding-field region and RF fields, which may underlie these oscilla- 


tions. 


It FIELD-DEPENDENT PEAKS 


With the RF probe stationed ten or more inches from the gun anode, 
narrow-band peaks can be found in the noise spectrum of the beam 
The amplitudes of these peaks increase and their frequencies decrease 
with decreases in the magnetic field. For each probe position, the process 
of finding the peak of greatest amplitude involves repeated adjustments 
of the focusing field, the magnetic field at the cathode, and the receiver 
frequency. 

When the fields have been so optimized, it is found that the probe is 
located at or near the first beam-diameter minimum, following that at 
the entrance to the drift space. When the field is doubled, and the 
“‘tuning”’ process repeated, the greatest peak is found to have about twice 
the frequency of the first, and the probe is found to be located at or near 
the second beam waist. It is convenient, therefore, to think of these 
peaks as “proper” frequencies of the N 1, ete., modes of the rippled 
beam, where N is the number of ripple wavelengths between gun and 


probe. 
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Fic. 1 Frequencies and amplitudes of several narrow-band UHF peaks 
measured at a fixed probe position, about 16 inches from the gun anode. N 1s the 
number of beam-ripple wave-lengths between anode and probe. Other peaks have 
been observed at higher harmonics of the ‘‘proper’’ frequency (encircled points), 
and at about half that frequency 


As shown by the encircled points in Fig. 1, these frequencies range be- 
tween 1.03 and 1.06 times the calculated cyclotron frequency, and have 
amplitudes as high as 65 db above the background noise. The amplitudes 
decrease with increasing N, falling off as the minus two-thirds power of 


the frequency. 
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At each of these optimum field settings, several weaker “satellite’’ 
peaks can also be detected, most readily those at the cyclotron frequency 
itself, and at 0.707 times the latter; i.e., the “plasma” frequency, as 
shown in Fig. 1. In addition, smaller peaks have been repeatedly ob- 
served at harmonics (up to the sixth) of the proper frequency, and one 
at slightly less than half of that frequency. (When a proper frequency 
was simulated by means of a signal generator, only its first harmonic 
could be detected in the receiver output.) 

At the fields corresponding to N = 4 in Fig. 1, the cyclotron frequency 
(312 me) was found, but not the proper frequency. The highest proper 
frequency observed was 240.5 me, in the NV 3 mode. The proper- 
frequency peaks decrease with increasing focusing field, whereas the 
field-independent peaks excited in the electron gun tend to increase, at 
the far end of the drift region. 


IV SPATIAL DISTRIBUTION OF UHF NOISE CURRENTS 


In Figs. 2 to 5 are shown synchronized chart records of collector cur- 
rent, one or more UHF narrow-band peaks, and microwave noise power 
near 4,000 me — all as functions of distance from the electron gun, for 
the N = 1 to 4 modes, respectively. In all runs, the beam was pulsed 
with a 1,000-cycle square wave, and the collector aperture set at a 0.100 
inch diameter. The magnetic fields at the cathode and in the drift space 
were adjusted before each set of readings, with the probe at a common 
reference position, for greatest amplitude of some UHI peak. In Figs. 
2 and 3, these were proper frequencies, whereas in Figs. 4 and 5 they were 
field-independent frequencies. 

The content of these distribution curves can be summarized as follows: 

(1) At the low fields employed (none quite equal to the nominal 
Brillouin value), the beam ripples are quite large, both in amplitude and 
wavelength. 

(2) The proper-frequency traces have two or three maxima near each 
beam waist, and their amplitudes grow more rapidly with distance from 
the gun than any of the satellite frequencies. 

(3) The patterns of the cyclotron and “‘plasma’”’ frequencies do not 
differ significantly from those of the field-independent frequencies, and 
usually display two peaks near each beam waist. 

(4) The collector-current maxima decrease with distance from the gun, 
although their minima change little. (The first maximum is sometimes 
flat-topped due to beam interception before it enters the drift space.) 


The rate of decrease of these maxima, and the rate of increase of proper- 


frequency amplitude, are greater, the longer the ripple wavelength. 
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Fic. 2 — The fields have been adjusted for maximum amplitude of the V = 1 
proper frequency, 77.8 mc, at a reference probe position (z = 15 inches). The 
synchronized probe records indicate three distinct maxima of this proper fre- 
quency near the beam waist. 
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Fria. 3 — The longitudinal distributions of the N = 2 proper frequency, 157 
me, as well as its ‘‘satellite’’ 76.2 me, and microwave noise power, are shown here, 
with fields adjusted for greatest amplitude of the proper frequency at z = 15 
inches. 
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ric. 4 — The fields have been adjusted for maximum amplitude, at the same 
reference probe position, of a wave excited in the diode, with frequency unaffected 
by the magnetic field, 165 me. The cyclotron frequency for this field is 240.2 me. 
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Fia. 5 A procedure similar to that in Fig. 4 was followed, with four ripple 
wavelengths between anode and reference plane. The cyclotron frequency here 
is 329.5 me. 
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(5) The patterns of microwave noise power resemble blurred envelopes 
of the UHF traces. 

Some idea of the transverse distributions of UHF noise power and 
electron-current density, in a region of strong proper-frequency exci- 
tation, is given in Figs. 6 and 7. The measurements were taken by mov- 
ing a small aperture in a broad are through the probe centerline, just 
in front of the probe aperture. In both illustrations, the relative noise 


“ce 


power has been “normalized” to compensate for variations in electron 
current traversing the RF gap. 
The curves of Fig. 6 are typical of most such measurements. The beam- 


current density varies smoothly through a single broad maximum, and 
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Fic. 6 — Simultaneous point-by-point measurements of collector current and 
relative noise power, obtained by moving an 0.013-inch diameter aperture in 4 
broad are through the probe centerline. The probe is stationary, about 18 inches 
from the gun anode, and the fields have been adjusted for maximum amplitude of 
the proper frequency, 148.5 me. The cyclotron frequency is 143.8 me 
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the noise-power density is greatest at the rim of the beam so defined, 
and least near its center. No evidence of azimuthal periodicity was 
found. The curves of Fig. 7, which are less typical, indicate five distinct 
peaks of RF power, despite a nearly symmetrical pattern of collector cur- 
rent. At the time of this measurement, cathode emission may have been 
uneven, due to coating damage by ion bombardment. 

In the rippled beam on which these measurements were made, the 
ratio of flux encircled at the cathode, to that in the drift space, was 
very small for most electrons. One would, therefore, expect the trans- 
verse noise-power distribution in this beam to resemble that in a smooth 
Brillouin beam.” The noise power expected when a pinhole aperture 
is located at the beam center can be compared with that when the aper- 
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Fic. 7 — Transverse distribution measurements similar to those of Fig. 6 
rhis pattern was obtained a week later than that of Fig. 6, and the cathode was 
operated at a higher temperature. The cyclotron frequency would be 153.2 me 


for the field used. 
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ture straddles the beam rim, by taking the beam area exposed in the 
first case to be that of a sector of angle 6, and the length of beam sur- 
face in the second case to be that of the corresponding are: 


Rf current sample inside of beam _ 6b’ J. _ bw — Bu) -Lo(8b) 


Rf current sample at rim of beam ~ 26bG, — -2ul,(Bb) 

Here b is the beam radius, and J/,, G. the longitudinal components of 
volume and surface current densities, respectively. Jy) and J; are modi- 
fied Bessel functions, 6 is the propagation constant, uw the beam veloc- 
ity, and w the radian frequency. For the frequencies and beam radii 
employed in these measurements, this ratio is very much less than unity. 
Thus the pattern of Fig. 6 is in accord with this mode distribution. The 
multiple peaks of Fig. 7, however, do not conform to this picture, and 
are not understood at present. 

As most of the RF power is concentrated near the rim of the beam, 
the question arises whether the double and triple peaks, in the longi- 
tudinal distribution patterns of Figs. 2 to 5, are not due to the probe 
aperture breaking through the beam rim. However, the dip between 
adjacent noise peaks is too great to be explained on the basis of reduced 
partition noise or weakened gap coupling, assuming the beam diameter 
there to be less than the gap diameter (0.100 inch). Moreover, double 
peaks occur even when the beam diameter exceeds the RF gap diameter; 
for instance, near the last three beam waists of Fig. 5. (When all of the 
beam is transmitted by the 0.100 inch aperture, the collector-current 
peak is flat-topped.) It seems likely, therefore, that the double and 
triple peaks correspond to peaks of amplitude over the entire beam 


cross-section. 


V PROPAGATION ALONG THE RIPPLED BEAM 


To find an explanation for the multiple peaks of space-charge cur- 
rent, a small-signal, slow-wave analysis of wave propagation along the 
rippled beam can be made, in which the special features of these experi- 
ments are exploited: long ripple wavelength, effectively no flux at the 
cathode, and low frequencies. The first of these features suggests that 
the propagation constants can be evaluated at each cross-section plane 
as though the beam were uniform, despite the presence of radial ve- 
locities. In addition, the space-charge density is assumed constant at 
each cross-section, and the electron flow laminar. 

With these assumptions, the beam can be regarded as a fluid of mov- 
ing charge, with a single-valued velocity at each point in space, as fol- 


lows: 
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where 


Here, r, 6, z are the polar cylindrical coordinates, w. = B the angular 
cyclotron frequency corresponding to the longitudinal focusing field B, 
and f(z) a funetion describing the amplitude and spatial periodicity of 
the beam ripple. The experimental data indicates that the potential 


variations along the beam axis are negligible, permitting the assumption 


that the longitudinal velocity, u, is constant. MKS units are used. 

Consistent with the distribution pattern of Fig. 6, the ac field can 
be represented by an axially-symmetric potential function, similar to 
that for the smooth Brillouin-flow beam: 


V ~ Ip(yr) exp j(wt — Bz), 
Kk = —grad V. 6 


The ac equations of fluid motion are obtained by adding a small ac 
increment to each of the steady-state velocity components. In addition 
to the space-charge field, the ac electric field contributes forces acting 
on the charged medium; those contributed by the ac magnetic field 
are neglected: 


d 


t (vo + 3B) —nl—grad V — grad Vo + (vy + B) X Bi. 
at 


The ac velocity is distinguished by a tilde, and the de velocity by a 
zero subscript. Here » = e/m is the charge-mass ratio of the electron, 
a positive quantity. As all ac quantities are functions of spatial positions, 
their time differentiation (indicated by a dot) is equivalent to multi- 
plication by j(w — Bu), written jos for brevity. 

The components of the force equation are expanded as follows: 
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(9b) 


OV 


Oz’ 


bs a]. , 
ES + U, : | l jnB\ ' (10b) 
or 


An expression for the ac space-charge density, ~, can be obtained in 


(10a) 


terms of its steady-state counterpart, p) , by means of the charge-con- 
servation equation: 
Op 
or 


—vp:grad p — d-grad po - po div } — p div v% 

—vo:grad po — po div v% . 
As the beam diameter changes slowly, the de space-charge density at 
each plane is taken to be inversely proportional to the square of the 
radius, b: 


Opo 2p» ob 
grad po —_=— 
Oz b O02 


vo grad p 


div vo (13) 


(Aa é me a (14) 
"\Or r ‘ r or ; as iki J Bur ; 


At low frequencies (the UHF region), 


av _V 


or / 


as (yr) < 1. This inequality is also true of other ac quantities pro- 
portional to V, such as 6, and pf, and with a small error can be assumed 
to be true for ¢,. When the operator 0/dr is omitted from (8), (9), (10), 


and (14), it is possible to solve explicitly for j in terms of pp and V. 
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The laminar-flow rippled beam can be described by the particle 
trajectories, as follows: 
ri roi(1 + 6 cos B,2), (15) 
where 7»; is the maximum radius for the particle considered, and 

0 < 6 < 1. For this model of the beam, 
b —dw, sin Bz 


b- =1+6cos B.2’ 


| ov w-B3.6(6 + cos B,z) I7 
° ‘) 

Az (1 + 6 cos B.z)° 
The region of interest, judging from the observed peak locations, is 
not at the mid-plane of the beam waist, where v, = 0, but on either side 
is greatest. It is readily found that, at these 


of that plane, where | v,/r 


positions (1/r) (dv,/dz) is zero, and (14) can be written 


This can be combined with Poisson’s Equation, 
| 


AV = (¥ — B)V 


y 


to furnish a relation between y and 8: 


where R = w, /a, and w, — npv/e, the square of the angular plasma 
frequency. 

At the beam boundary, r b, the continuity of the tangential field 
components and the change in radial electric displacement can be ex- 


pressed in the form of an admittance equation: 


(2 7 ) i av |" ” 
V \ Or € bee |V @r “- ‘ 


Here I refers to the beam, 0 S r S b, and II to the space between beam 
< r S a. The surface charge 


and the concentric conducting tube, 6} 
layer, g, takes account of the surface ripple, of amplitude 7 
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” Jpovr 
por ” _ 
Wh 


R(dV/dr) 


The appropriate potential functions in I and II are reduced by means 


of the low-frequency, or thin-beam, approximation, as follows: 
[1 av’ yl \(yb) a yb 
V orto Io(yb) as - i 
AA _ g | 1:(6b)Ko(8a) + Io(8a)K1 (80) 
V or |p ” | Io(Bb)Ko(Ba) — [)(Ba) Ko(Bb) 


Bb l 
~ — 
_ on b In a/b 
where the following small-argument approximations have been used: 


Ko(xz) & —Ina, 


(24) 


: x l 
Ki(2) & = Ina + 


x 


The boundary equation thus provides a second relation between y and 8: 


R Bb l 


| 7? e a’ (25) 
( l J ) b In 
Wpr h 


lor the smooth beam in Brillouin flow (v, = 0), the boundary equa- 
tion, to the same low-frequency approximation, is as follows: 


"tis 2 

R Wp ~ 

to ° . 92 
(w — Bou)- o a (26) 

(Bob)? In 

b 


To see how the beam ripple affects the propagation constant, it is 
sufficient to find its first-order effect; i.e., to assume relatively small 
radial velocities and find a solution for 8 which is not very different 
from its value, 8 , for the smooth beam: 


b Bo +e 6 = B, =. Bo + é, 


where 
Wp 


6 | <K Bo } Ba = “ur PR,” B. 
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In addition, | v,/asr | is less than unity, and | 2v,/8ur | can be neglected 


entirely. With these assumptions, the boundary and characteristic 


y Ro 


F-—-R ° 


! 


equations can be combined to solve for 8: 


where 


Wpr 


Utilizing the low-frequency condition, | R°| > | R| > 1, this equa- 
tion can be reduced and, after some algebra, solved: 


Bo 4) . Bold + Bq) ~ 
B R (Bo oa 6) ( +B,) aad 


Uy 


8=6B+£8,- 7 (29a 


Uy 


By = Bb, —_ Bo + (29b) 


, Qur’ 

These expressions show that the current in the slow wave (I,) will 
grow when (v,/r) is negative; i.e., when the beam is contracting, and 
decrease during its expansion. The fast wave (I;) will do the opposite 
In probe measurements along the beam, the detected ac power is pro- 
portional to the square of the total space-charge current, which has 
the following dependence on time and distance when the amplitudes of 
both waves are initially equal: 


.ur 


; V,2 
UI, + Is) = 2Imax COS (wt — 68,2)-cos (8,2) -sinh (; ). (30) 


In UHF noise-power measurements along beams with long ripple 
wavelengths, the two planes of maximum + (v,/r) are separated by 
only a small fraction of a space-charge wavelength. Therefore, cos 8,2 
at the first of these planes is only slightly larger than at the second. 
Thus, two peaks of current are observed, in agreement with (30). By 
contrast, in rippled-beam amplification at microwave frequencies, 
shorter ripple wavelengths and smaller ripple amplitudes are employed. 
Then (v,/r) varies nearly sinusoidally over the ripple wavelength. For 
maximum net gain per ripple, maximum negative (v,/r) is adjusted to 
coincide with the plane of cos 6,2 = 1 (maximum current), and maxi- 
mum positive (v,/r) at the current minimum, half a wavelength beyond. 





872 THE BELL SYSTEM TECHNICAL JOURNAL, JULY 1957 


The gain constants in (29) are independent of frequency. The net gain 
per ripple wavelength, however, will vary with frequency, depending on 
how closely both the current maxima and minima coincide with the 
regions of maximum + (v,/r), respectively. This is a statement of the 
‘resonance’? condition between ripple wavelength and half the space- 
charge wavelength, which emerges from one-dimensional analyses” of 
this gain mechanism based on transmission-line analogies. 

Such analyses generally assume small-amplitude sinusoidal variations 
of the reduced plasma wave number, 8, , along a one-dimensional beam 
in a longitudinal ae field with no losses. Periodic variations in either beam 
or wall diameters, or beam velocity, cause the beam “impedance’’ to 
vary periodically, imparting to it narrow-band filter-like properties 
equivalent to narrow-band signal gain. From another point of view," 
these periodic impedance changes couple the fast and slow space-charge 
waves to each other intermittently, thereby effecting an energy trans- 
fer from the fast to the slow wave. As this coupling is lossless, 7, increases 
and J; decreases with drift distance, in such a way as to keep their product 
constant. Then the product Jmax/min increases, and the ratio Jimax/min 
correspondingly decreases. In the case of noise-power amplification, two 
uncorrelated space-charge standing waves are present. Because the two 
slow waves cannot simultaneously be amplified at the expense of the two 
fast waves, the product /max/min Must remain constant. 

The observed noise-current patterns in rippled-beam amplification, 
however, are characterized by a nearly constant ratio Imax/Imin , and an 
increase in the product Imaxtmin along the beam, despite the fact that 
the beam voltage is fixed. This apparent contradiction can be resolved 
by a closer look at the energy-exchange processes. 


‘ 12 . . . 
Chu’ has shown that the kinetic power flow in space-charge waves 


(the major part of the total power) is equal to the difference in powers 
carried by the fast and slow waves. This is equally true of beams with 
transverse motions and fields.® In rippled-beam amplification, whether 
analyzed as a modulated linear beam or at each beam cross-section sepa- 
rately, as here, the propagation constants are found to be complex con- 
jugate quantities, whose real parts describe the ordinary fast and slow 
waves of a uniform beam. From either point of view, therefore, a de- 
crease in J; and an increase in J, signifies an increase in the negative 
kinetic power flow carried by the waves, or a decrease in the total kinetic 
energy of the beam. 

As shown in (29), the gain constants are proportional to v, , indicating 
that the de energy transferred to the waves when the beam contracts 
could only have come from the radial kinetic energy, not the longitudinal. 
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The direction of energy transfer is reversed during the subsequent beam 
expansion. If the ripple were perfectly symmetrical, therefore, and the 
de-ac energy exchange perfectly reversible, the net effect of a beam ripple 
would be zero. Neither of these conditions is quite true in actual beams. 
Rippled flow is never truly laminar, and | v,/r | usually decreases with 
drift distance as the flow loses coherence; i.e., it is greater in beam con- 
traction than in the next expansion. This by itself would produce a net 
gain per ripple in /, , and a net loss in J; , of equal amounts. In addition, 
however, unavoidable small non-linearities in electron motions prevent 
all of the ac energy in a de-amplified wave from being converted back to 
de kinetic energy. Thus it is possible for both the fast and slow waves to 
increase in a ripple wavelength, the latter always more than the former 

The greater gain of the slow wave entails a loss of radial kinetic energy 
in agreement with the observation that the ripple amplitude always 
decays more rapidly when rippled-beam amplification takes place. The 
incomplete reversibility of the ac-de energy exchange probably accounts 
for the observed increase in Jmax/min for noise currents. Finally, the net 
amplification of all of the space-charge waves, fast as well as slow, is in 
accord with the observed near-constancy of the ratio Jax // for 
microwave-frequency noise, despite increases in the product JmoxJmin of 
30 db and more. 


VI ORIGIN OF THE PROPER-FREQUENCY PEAKS 


Of the various peaks in the beam’s noise spectrum, described in Seec- 
tion III and Fig. 1, those with “proper frequencies,” slightly above the 
cyclotron value, are so large in amplitude that even an approximate 
analysis should be able to account for them. To do so, a“ working model” 
of the beam is needed, which conforms to the experimental conditions 
which existed during the observations: 

(1) The peak intensities were greatest near the middle of each beam 
waist, and decreased with decrease in ripple amplitude. 

(2) The focusing field was below the nominal Brillouin value. The 
field at the cathode, B., was finite and opposed to the main field, B. 

(3) Colleetor-current measurements along the beam axis showed the 
ratio of maximum to minimum current to be greater, the smaller the 
aperture. 


(4) The gas pressure was about 10 ‘ mm Hg. The beam was pulsed 
with a 1,000-cycle square wave. 


Item (3) indicates that the flow was non-laminar; and Item (4) in- 
dicates the presence of positive ions. All the items are consistent with 
the following picture: 
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In a beam with large ripples, nearly all electrons have their maximum 
radii and zero radial velocity at the same z-plane. Those with sufficiently 
large maximum radius will have enough transverse kinetic energy to 
surmount the space-charge forces at the beam waist, and pass through 
or close to the axis. Others, with smaller maximum radii, will spiral 
about that axis. Dolder and Klemperer” have observed a similar di- 
vision of electrons into ‘‘crossovers” and non-crossovers, in electron- 
optical systems without magnetic fields. 

Positive ions tend to neutralize the electronic space charge at the 
beam waists, broadening the region in which crossover occurs. The cross- 
over trajectories thereupon overlap one another, resulting in multi- 
valued transverse particle velocities in this region. In a first-order 
(linearized) study of wave propagation along the beam, one must re- 
place the actual multivelocity charge motions with a single “‘fluid”’ of 
charge, whose velocity at any point is the average of the particle ve- 
locities there. It is clear that the z-velocity of the stream is u, and the 
radial velocity zero. The tangential velocity, ve = (6r),,., however, is 
more complicated. 

Owing to the partial or total neutralization of electronic space charge 
at the beam waists, and their large radii elsewhere, the crossover elec- 
trons will encounter virtually no space-charge forces in their paths. Their 
transverse paths will consequently be circles about fixed centers, de- 
scribed with angular velocity equal to the cyclotron frequency. Their 
angular velocity about the beam axis is given by Busch’s Theorem: 


where 


. (B. 
K x Te (5 ) Tmax" n in (351 ) 


is a positive quantity, as B./B is negative. Here, r, is the radius at which 


a particular electron left the cathode, and r is its radius in the drift 
region. The angular velocity, 6, is greater than w./2 at all times, and 
exceeds w, in the waist region of the beam. The average value of vp» at 
any point here, therefore, is greater than w.r and presumably varies from 
point to point in some unknown way. 

If ve is left unspecified, and the assumptions adopted of zero space- 
charge forces and radial velocity over a finite length of beam: 

OVo dv, 


—a ', 0, = (), (32) 
. or dt 


the radial component of the force equation (7) in Euler coordinates can 
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(+) Ve 
= —— — weds, 
dt t r 


vo =O and wr. 


be written as follows: 


Thus, the radial “balance”? conditions (32) are consistent with either of 
two values for ve , of the equivalent stream with single-valued velocities 
As it develops that either of these values leads to the same result, the 
first one will be used here for simplicity, ve = 0. 

An ac traveling wave along this beam cannot have any 6-dependency, 
because the beam has no single value of angular velocity @, which might 
remain in synchronism with that of the wave. Thus, the perturbed 
dynamics equation (7) can be expanded, with the assumptions of an 
axial-symmetric ac field given by (5) and (6), a stream with steady- 
state velocity (0, 0, uw), constant space-charge density po , and no space- 
charge forces, as follows: 

aV 


= 7 - - Wee , 
or 


+e 


ali 02 


These are solved for the ac velocity components: 


— J Nw 0 | ; 


o « ; 


— We OF 
(36) 
(37) 
With grad p) = div “0 = O, the charge-conservation equation (11) can 
be solved for p: 
Op 


ot 


— a Y 8 
p = div } = npo| — “elie OC 
Wh a = Wh” 


At very large ripple amplitudes, it is a fair assumption that the density 
of non-crossover electrons is negligible relative to that of crossovers in 


— Vo-grad p — po div é, 


this region. Poisson’s equation (19) can then be combined with the 
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above expression to obtain the characteristic equation: 


o) 


In a frame of reference moving with the stream, w’ is 0, 8’u’ is 0, and 


w,’ = w. Then, 


> 


(w — We — Wp Jor (40) 


(w? — w,”)(w? — w,”) 
and 8’ becomes an infinite imaginary quantity when w is w, . The phase 
velocity in the moving frame is infinite, as the real part of 8’ is zero; 
therefore the phase velocity v, in the rest frame is also infinite. Thus, 
there is no Doppler shift in the ‘‘resonant”’ frequency observed in the 


rest frame: 


(41) 


As the actual beam has a 2-velocity spread, the field is never perfectly 
uniform, and as the calculation is valid for small ac quantities only, the 
discrepancy between this result and the observed “‘proper’”’ frequencies, 
which were 1.03 to 1.06 times the cyclotron value, is not unexpected. 

The singularity in (40) is seen to disappear when w, = 0. This indi- 
cates that an exact calculation would show the gain constant (—j8) 
to increase with po, the density of the crossover electrons. Their tra- 
jectories, described by (31), and the absence of space-charge forces are 
such that K = fnax?min; that is, the greater rmax , the smaller rnin, 
the distance of closest approach to the axis. Thus, a larger ripple ampli- 
tude (permitted by a lower magnetic field) produces a greater electron 
density in the waist region, and accordingly a greater oscillation ampli- 
tude at the resonant frequency, as observed. 

The foregoing mathematics describes a form of resonance, the infinite 
phase velocity corresponding to longitudinal ‘‘cutoff” in a waveguide. 
Unlike a waveguide, however, the disturbance increases rather than 
attenuates along the axis, due to the transfer of de kinetic energy (repre- 
sented by ve /r) to the ac fields (excited by noise fluctuations at the 
cathode), at the cyclotron frequency a, 

Except for the direction of energy transfer, the situation is analogous 
to that of a low-pressure gas in a uniform magnetic field, when stressed 
by an impressed ac field of varying frequency. It has been found that 
the breakdown field at the cyclotron frequency is very much less than 
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at other frequencies. ° Here the energy supplied by the ac field is coupled 


most effectively to the free electrons at the resonant frequency, increas- 
ing their de kinetic energy until the gas breaks down. The circular ac 
charge motions due to the de magnetic and the ac electric fields are 
superimposed on high-velocity random motions, similar to the radial 
motions in the drifting beam. 

The UHF peaks observed at harmonics of the proper frequency may 
simply be due to the non-linear character of the beam, when excited by 
the high-level fundamental oscillations. The other faint satellite peaks, 
near 0.5 w, and 0.707 w. , seem to be associated with the unneutralized 
space-charge density at the beam waist. 

The conspicuous role played by crossover electrons in the waist region 
of rippled beams, due to the tendency of their orbits to overlap there, 
leads one to re-examine their influence on rippled-beam amplification. 
As seen in the previous section, this gain process depends on the average 
value of (v,/r) at each cross-section plane of the beam. The fraction of 
all electrons which penetrate to the beam axis depends on competition 
between the unneutralized space-charge forces and the particle’s trans- 
verse kinetic energy. An increase in positive ion density tends to make 
the potential depressions at beam waists broader and shallower, and 
thereby increase the number of crossover electrons as well as the axial 
distance over which they reach the axis. The net effect is to reduce the 
average value of | v, | over a greater portion of the ripple wavelength, and 
thus reduce the net gain of the space-charge wave. This may explain 


“ 


why the “growing noise’? phenomenon tends to be inhibited by an in- 


crease in positive ion density. 


VII CONCLUSIONS 


Evidence is found of oscillations with frequencies in the 10- to 500-me 
region inside of an electron-gun diode. There is some basis for associating 
them with electron-field interaction in the retarding region of the diode 
Another type of narrow-band noise peak is found near the waists of a 
strongly rippled beam in a longitudinal magnetie field, with frequencies 
proportional to the field strength. The strongest of these, at about 1.05 
times the angular cyclotron frequency, w, , as well as its harmonics, can 
be explained by the resonant behavior of a short section of the beam, in 
which the average transverse velocity is nullified by overlap in particle 
orbits. Fainter satellite peaks, near 0.5 w. , 0.707 w. , and w, , respectively, 
accompany the dominant frequency. 

In a drifting beam launched from a shielded electron gun and focused 
by an axial field, the transverse distribution of noise (or signal) intensity 
is found to agree with that predicted for ideal Brillouin flow. Despite 
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the presence of thermal motions and beam ripples, the ac power is found 
to be concentrated chiefly at the rim of the beam. Occasionally, several 
concentric rings of noise maxima are found within the beam, possibly 
due to unusual cathode conditions. 

When the ripple wavelength is very long, two maxima of noise power 
are observed to flank each beam waist. A first-order calculation of wave 
propagation along a rippled laminar-flow beam accounts for this pat- 


tern by showing that space-charge waves grow at the expense of de 


kinetic energy in the radial charge motion. In rippled-beam amplifica- 
tion of noise, the product Jinax/min has been found to increase, and the 
the ratio Jimax/Jmin remain nearly constant, because both fast and slow 
waves are amplified, the former less than the latter, and because the 
wave coupling is not lossless. 

Positive ions tend to collect at the waist of rippled beams, thereby 
extending the region in which electrons pass close to the axis, instead 
of circling about it. The overlap of their orbits leads to net cancellation 
of radial charge motion, and hence a reduction in rippled-beam ampli- 
fication. This may explain why positive ions tend to inhibit the “‘grow- 
ing noise’ phenomenon. 
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Distortion Produced in a Noise Modulated 
FM Signal by Nonlinear Attenuation 


and Phase Shift 


By S. O. Rice 
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An expression is given for the FM distortion introduced by a transducer 
whose attenuation and phase shift depend upon the frequency in an arbitrary 
way. This expression appears to be difficult to evaluate, but it yields useful 
approximations for the second and third order modulation terms. In all of 
the work, it is assumed that the distortion is small compared to the signal, 
and that the signal can be represented by a random notse having the same 
power spectrum, 


INTRODUCTION 


A number of workers have been concerned with the problem of com 
puting the distortion introduced by a transducer when an FM wave 
passes through it. Some of the earliest results were published by Carson 
and Fry’ and by van der Pol.’ Several contributions to the subject have 
been made recently in connection with studies of microwave radio 
systems. 

An excellent paper on this subject has been published recently by 
R. G. Medhurst and G. IF. Small.” Although their results differ consider- 
ably in form from those given here, they are nevertheless closely related 
to ours — their “sinusoidal variations of transmission characteristics” 
being special cases of our ‘‘nonlinear attenuation and phase shift.” 

Here we treat the problem by applying a method used in a recent 


4 ‘ ° rr ‘ 
paper to study the distortion produced by an echo. Two assumptions 


are made, (1) that the distortion is small compared to the signal, and (2) 
that the signal can be represented by a random noise which has the 
same power spectrum as the signal. In Section I, we review some known 
results and put them in a form suited to our needs. Sections II and III 
are devoted to the derivation of our main formulas. The principal result 
is given by the triple integral (3.2) for the power spectrum of the dis- 
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tortion. Unfortunately, the integrals are difficult to evaluate. However, 
it is possible to obtain approximations for the second and third order 
modulation terms. These are given in Section IV. Some miscellaneous 


comments are made in Section V. 


I APPROXIMATE EXPRESSION FOR THE DISTORTION 6@(t) 


Let the FM signal be ¢’(¢) = dg/dt (for phase modulation the signal 
would be g(t)). Then the FM wave is the real part of 


v(t) = eb Pere’ (1.1) 


where p 2rf, is the carrier frequency. Let this wave pass through a 
transducer having attenuation @ and phase shift 8, where a and 8 are 
even and odd functions, respectively, of the frequency f. When a unit 
impulse of voltage 6(¢) is applied to the transducer input, the output is 


a 


g(t) "ieee | (1.2) 


=) 
For physical systems, g(¢) 1s zero for negative ¢. 
When v,(f) is applied to the transducer input, the output is 


a oO 


Volt) = v(t’ )g(t — t’) dt’. (1.3) 


bd x 


When vo(f) is applied to an FM receiver, the detector output consists of 
the original signal ¢’(¢) plus the distortion 6’(¢) introduced by the trans- 
ducer. Comparison with (1.1) shows that 6(4) may be obtained by 


solving 
‘im = oh (1.4) 


when p, g(t), vo(t) are assumed to be known, and V(é), 6(¢) unknown. 
When V(t) is taken to be positive, (1.4) determines @(t) except for an 
additive term of 27n where n is an integer. 

We now assume that the transducer acts like a good transmission 
medium in that the output differs but little from the input. More pre- 
cisely, we assume 

u(t) — v,(t) | K 1. (1.5) 
Since | v,(¢) | = 1, it follows that | v(t) | & 1. Transducers having ap- 
preciable attenuation and delay may be regarded as two transducers in 
tandem, one with constant (independent of f) values of a and 8/f which 
are roughly equal to those of the original transducer, and the second 
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with variable a and 8/f. The first transducer produces no distortion of 
the signal, and if condition (1.5) is satisfied by the second, the con- 
siderations of this paper will apply. 

Equation (1.4) may be written as 


4 6(t) 


V (tye 


vo(t) /v;(t) 


so that 
. v 
A(t) = Im log 
. 


When we write 

vo(t) /v,(t) = 9 — v,(t)]/v,(0), 
expand the logarithm in (1.6), and use (1.5), we obtain our approximate 
expression for @(¢): 


A(t) = Im [vo(t) — v,(t))/v.(t) = Im volt) /v, (8) 


= Im [v;(t)] : [ v(t’) q(t — t’) at’ 


e =) 
is 


Im exp [ip(t’ — t) + ig(t’) — ip(d|g(t — ¢) dt’. 


2) 


So far there is nothing essentially new in our work.” 


If AUTOCORRELATION FUNCTION OF 6(¢) 


In Section I, ¢’(¢) could be any reasonable sort of signal. In the follow- 
ing work we assume that it is a Gaussian noise whose power spectrum, 
wy (f), is given to us. The power spectrum of ¢(f) is 

wolf) = wy (f)/(2xf)’, 2.1) 


¢ 


and its autocorrelation function is 


y, = [ we(f) cos 2xfr df. 


“0 


We have written y, instead of ¥(r) or R,(r) to simplify the appearance 
of the formulas which occur in our work. 

Our problem is to find the power spectrum, w»(f), of the distortion 
6(t), given w,(f). The method of solution is much the same as that used 
in Reference 4. We first find the autocorrelation function Rg(r) of A(t 
and then obtain we(f) by taking the Fourier cosine transform of R¢(r). 
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Let the last integral in (1.7) be F(t) so that 6(t) = Im F(t). Then 
6(t)O(t + 7) = $ Re {F(OF*(t + 7) — FOF(t + 7)} (2.3) 
where F*(¢ + 7) is the complex conjugate of /(¢ + 7). The autocorrela- 
tion function of @(t) is obtained by averaging over the ensemble of the 
noise functions ¢(t): 
te(r) av O(Ue@t + 7) 


| , ” . , 
=av = Re | dt | dt exp [ip(t — 0) 


2 


‘g(t — t)g(t + 7 — t”)[exp [—ip(t” 


— w(t”) + v(t + 7)] — exp lip(t” -—t— +r 
1o(t T r)| 


Since g(t) is real, g*(t) = g(t). The averaging process may be carried 
out by a method analogous to that used in Reference 4. The formula to 
be used is 
av exp [ig(t’) — ig(t) + tag(t”) — tag(t 4 

exp [—W(1 + a) + y, vr + ay, 
+ abr — abe + Cert, | 


where a is either —1 or +1, and y, is an even function of +. When (2.5) 
is used in (2.4) a double integral for Re(r) is obtained. The substitutions 


(2.6) 
R, = Vriz oo Wras = Weiss 


convert the double integral into 
| a oO x 


Ro(r) = = Re dx dy g(xa)e °”* 


- e x e x 


ipy+R 


-gly)le — ¢ 
The symbol R, is chosen to agree as closely as possible with the notation 
of Reference 4. There FR, was the autocorrelation of the random func- 
tion, v(t), where v(t + T) = g(t) — g(t + T), T being the echo delay. 
Here, FP, is the average value of the product, 
lot) — g(t + y)| [et + 7) — et + 74+ 2) 


which becomes the autocorrelation function of v(t) when y = x 
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It may be verified that the expression (2.7) for Re(r) is an even func- 
tion of 7, as it should be. Expression (2.7) is the autocorrelation func- 
tion we set out to find. 

The distortion @(¢) has an average value, 6, whose square is Re(<). 
Since g(t) is a noise function, its autocorrelation function y, goes to 
zero as rt approaches ~. Hence, Rs(~) is given by the expression ob- 
tained from (2.7) by setting R, = 0. The autocorrelation function of 
6(t) — @is 


Re_o(r) = Retr) — Roo) 
| . x a 
= Re dx dy g(axye 


/ l PU R 
‘qQhy)e ( 


Ill POWER SPECTRUM OF THE DISTORTION 


Since @(¢) has an average value which is generally not zero, its power 
spectrum, w(f), has a spike of infinite height at f = 0 corresponding to 
the power in the de component @. When this spike is subtracted from 
we(f) the remainder is the power spectrum of 6(¢) — @ given by 


| Ro_@(r) cos 2xfr dr 


70 


When we use (2.8) and note that R»_4(7) is an even function 
obtain 


cc ee ne 
weal f) dx dyg(a)g(ye °°" [cos (px' — py) 
im \ om 
-(e-* — 1) cos (px + py) le *» _ 1)] cos 2xfr dr. 

Reasoning similar to that given in Reference 4 shows that the inter- 
channel interference spectrum, w-(f), (i.e., w-(f)Af is the average amount 
of distortion power received in an idle channel of width Af centered on 
the frequency f, all other channels being busy) may be obtained from 
(3.2) by replacing (e*“" — 1) by (e*”" + R, — 1). 

The power spectrum of 6(f) — @ may be regarded as made up of 
modulation products of all orders. It turns out that the contribution of 
nr” order products is given by the integral of the R,” terms obtained 
from the power series expansions of exp [+2,]. 


IV FIRST AND SECOND ORDER MODULATION TERMS 


Here we shall study the first and second order modulation terms. 
These arise from the first and second powers of FR, in the expansion of 
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the quantity within the square brackets in (3.2): 
2R, cos px cos py + R,° sin px sin py. (4.1) 


The integrations with respect to 7 may be performed with the help of 


pe oC 


Wel f ) ioe fb 
Yr45 COS 2afr dr = .~ ~Roe (4.2) 


= - 


1+bWrie COS 2nfr dr 


Re | F du w,(u)w,(f — u) exp {—72n[bu + ce(f — u))} 
which follow from (2.2) and the fact that we have defined w(—f) to be 
equal to w(f). In our notation the total power in a random noise function 
is the integral of w(f) from f = 0 to f 0. 

The first order modulation term is obtained from (3.2) by replacing 
the term within the square bracket by 2R, cos pr cos py. When the ex- 
pression (2.6) for R, is used, the integration with respect to 7 may be 
performed with the help of (4.2): 


aX 


: w.(f) onizf ioey 
R, cos 2nfr dr = — Re [(e™™ — 1)( ww _ 1)). (44) 


This leads to the following expression for the first order modulation 
term in (3.2) 
we(f) dxg(x)e *°*** cos pxle ie} (4.5) 
J—2 

This is the quantity which is to be subtracted from we_4(f) to obtain 
the interchannel interference spectrum w,(f). 

The second order modulation term is handled in much the same 
manner. With the help of (4.3) it may be shown that 


a al 


9 : l 
R, cos 2nfr dr Re du w,(ujw,(f — wu) 


x * x 


2riz( f—u = 1) 


— 1) (e 


2riy(f—u 
1) 


1) (¢ : = 


rom this it follows that the second order modulation term in (3.5 
» a 


a> duw,(u)w,(f — wu) dx g(a)e sin px 


oe. & I—x J—2 


9 ’ 9 f 
it 2eizu l)(e 2riz() 
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When yo — yz is so small that exp (—Yo + Wz) may be replaced by 
unity, as it is in some important practical cases, approximations may 
be obtained for (4.5) and (4.7). The integral in 2 may be expressed as 
the sum of integrals of the type 


a 2 


I g(x) rr de = [6 *) 04) 


e 4) 


r 


Ga +- iB, " 

nce 

| g(x)?” °**"** dx = G_, — 1B_q. 

Je 
The values of the integrals follow from (1.2) and the Fourier integral 
theorem. G and B are, respectively, even and odd functions of frequency, 
and G, , B, are their values at the frequency f = f, + a where f, = p/2m 
is the carrier frequency: 


G at frequency f, + a 


B at frequency f, + a 
In this way we get the approximation 
1 ‘we(f) (Gy — 20) + G 


for the first order modulation term, and 

| a) 
: duw,(uw.(f — w[(Gu — Gu + Gpu — Goppun — G; ; 

J—o (4.10) 


=f G 5)” op (B, + B u ote B, . oa B fru — B,; a B _ 2By)’| 


91 
am. 


for the second order modulation term. 

Expression (4.10) is an approximation to the second order modulation 
term (4.7). When most of the interchannel interference is due to second 
order modulation products, (4.10) is also an approximation to w,(f), the 
interchannel interference spectrum. The following remarks may be of 
some help in deciding whether (4.10) may be used. 

1. For the case of phase modulation and a “‘flat”’ signal band, the 
first of equations (5.3) shows that Yo and y, may be made as small as 
we please by choosing the signal power (as measured by Po) small 
enough. Since R, is proportional to Py , Po may be chosen small enough 
to make R,° and higher order terms negligible in the expansion of the 


)» 


integrand of (3.2) (unless there is some sort of symmetry which causes 


the second order terms to vanish). In this case the interference is mostly 


second order modulation and (4.7) is a good approximation to w.(f). 
l'urthermore, as Py approaches zero, exp (—yY + ¥-) approaches unity 
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and (4.10) becomes a good approximation to (4.7). Just how small Po 
has to be depends upon the signal bandwidth, f, , and the characteristics 
of the transducer. 

2. For the case of FM and a flat signal band, the second of equations 
(5.3) shows that even if Po is small, the difference Yo — Y, approaches « 
as |r| approaches «. To justify the use of (4.10) in this case it is neces- 
sary to take into account the behavior of g(¢), the response of the trans- 
ducer to the unit impulse 6(¢). For example, if the duration of g(x) in 
(4.7) is so brief that g(x) becomes negligibly small before —y) + yw, be- 
comes appreciably different from zero (which may be achieved by mak- 
ing Py small enough) then (4.10) is a good approximation to (4.7). 

3. When the attenuation, a, and phase shift, 8, are given for any 
particular transducer, the corresponding g(t) may be obtained from (1.2). 
Once g(t) and Yo — y, are known, the conditions under which exp (—yWo + 
W.) may be replaced by unity in (4.7) and O(R,°) terms neglected in (3.2) 
may be determined by direct examination of the integrals. 

As might be expected, the third order modulation results are quite 


» 


complicated. The third order modulation term in (3.2) is 


a p oO 
l mr 


df’ | df’’'w,(f')w,(f")w,(f’”’) 


4 Jw © Jw 


a 30 


|] dx g(x) cos pre *°""7(2" — I)(2e""” — 1)(2" — 1) 


= 9) 


where f’” = f ‘and z = exp (—72mxr). When yo is small this is 


approximately 


wef”) [H* + K*] (4.12) 


3!16 . 
where 
H = m(f’) + m(f’”) + m(f) + m(f - 


— mf — f’) — m(f — f”) — m(O) — mf’ +f 
mf) = G;+ G_;, n(f) = B,; — B_;, 
and K is an expression obtained from H by replacing n by m. 


V MISCELLANEOUS COMMENTS 


Here we make some miscellaneous comments related to the foregoing 
results. 
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If the transducer is perfect except for an echo, its response to a unit 
impulse 6(¢) is 
g(t) = 6(t) + ré(t — T) (5.1 


where r and T are the amplitude and the delay of the echo. The results 
obtained using (5.1) agree, as they should, with the results obtained in 
Reference 4. Of course, r must be assumed small compared to unity in 
order that condition (1.5) may hold. 

When the power spectrum of the signal is equal to a constant Py over 
the band (f., f,) and zero elsewhere we have for phase and frequency 
modulation, respectively, 


PM: w,(f) = Po, f. <4. 
FM: w,(f) Po /(2xf)’, 

When f, = 0 the autocorrelation functions are 
PM: vy, = Pofs(sin v)/v, 
FM: wo — ¥; A[—1 + cosv + vSi(v)], 
v 2nfir, A = Pofv(2xf) 


The mean square values of the signals are Pof, (radians)” for PM and 
Pof, (radians/sec)° for FM. If, for FM, o is the rms frequency deviation 
in eps (so that the “peak” deviation is, say, 4¢ eps) then (240) = Pofy. 
The difference Yo — y, is used in the FM case to avoid difficulty at f = 0. 
It will be noticed that our formulas are such that the ¥’s may be re- 
placed by (Ww — yYo)’s without altering the values of the various ex- 
ponents, etc. In microwave systems the quantity A is often small in 
comparison with unity. 

As an example of the use of the second order modulation approxima- 
tion (4.10) consider the case where the attenuation, a, is zero and the 
phase shift 8 = a.(f — f,)°/2 radians, a, being small. Then, since G & 
l1—a,8 > —8, we have G, ~ 1 and 

B, =~ —(6 for f 
o.4 


—ayu /2. 


When we take the FM case of (5.2) and substitute in the approxima- 
tion (4.10), the interchannel interference power spectrum is found to be 
! [" Po Po 
2!8. 


[0 + (Qasu(f — u))*] du 


p—f, (Qru)? (2r)?(f — w) 


‘ 4 9\2/9F 3 
(29) *(a2Po/2)°(2fs — f). 
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Dividing by w,(f) P)/(2nrf) gives the ratio of the interference power 
to the signal power 


(aoof /2)°(2 — f/fr) (5.6) 


where the relation Po» (23c)"/f, has been used to eliminate Py . Here 
o is the rms frequency deviation of the I'M signal in eps. The expression 
(5.6) agrees with results of some earlier work done at Bell Telephone 
Laboratories. In that work the second order modulation products were 
summed directly. 

It is interesting to apply the formulas given here to some of the 
cases considered by Medhurst and Small.” They have shown that when 
(in our notation) a = —rcos 2xfT and 8 = 0 the power spectrum of the 
distortion is 


] 


lecho 35 


we_a(f) = sin’ rfT[we_a(f 
and when a 0 and B r sin 2rfT, 
we_o(f) cos mfT [we—a(f)lecho - 


Here [we_a(f) echo is the power spectrum of the distortion due to a simple 
echo of amplitude r and delay 7 (corresponding to a —r cos 2nfT 
and 6 r sin 2xf7). Expressions (5.7) and (5.8) may also be obtained 
by setting the impulse response g(t) equal to 


; m r mn 
a(t) + <a(t — T) += 4(t + 7) 


in (3.2). 
The second order modulation approximation for the a —r cos 
2nfT, B 0 case may be obtained from (4.10) and turns out to be 


u)[2r sin pT sin xfT sin ruT sin r(f — wT) du. (5.9) 


It is seen that this contains the factor sin'xf7' predicted by (5.7). When 
(5.9) is applied to the FM case of (5.2) an integral something like (5.5) 
(but more complicated) is obtained. The ratio of the second order 
modulation interference power to the signal power is found to be 


2[r sin pT sin fT} (o/h) UK (5.10) 
where A is the quantity 


sin (y/2) sin (a — y)/2 i is 
: : dy (5.11) 
yYla— Yy) 
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tabulated in Table 4.2 of Reference 4 and 


a 2rfT, U 2rfol’. 


The parameters a and k that appear in the table are defined by 


a = f fo and I: = Sf,T’. 


These formulas serve to supplement the formulas and curves given by 
Medhurst and Small. 
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Self-Timing Regenerative Repeaters 


By E. D. SUNDE 


(Manuscript received March 29, 1956) 


In self-timing regenerative repeaters, a timing wave for control in pulse 
regeneration ts derived from the binary pulse train at each repeater with the 
aid of a resonant circuit tuned to the pulse repetition frequency. The timing 
wave can be made to exercise complete control in retiming of pulses inde- 
pe ndent of the received pulse train, or it can be combined with the received 
pulse train to provide partial retiming. The timing principles are discussed 
here for a particular type of self-timed regenerative repeater invented by 
Wrathall, in which a timing wave derived from either the received or the re- 
generated pulse train is combined in a particular way with the received pulse 
train. The regeneration characteristics of such repeaters as determined by 
various design parameters are investigated, together with the cumulation 
of timing deviations in repeater chains and the circuit requirements that 


must be met to insure satisfactory performance. 


INTRODUCTION 


Pulse transmission systems employing binary codes, such as PCM, 
have two inherent properties that are desirable from the standpoint of 
avoiding excessive transmission impairments by noise and other imper- 
fections in the transmission medium. For one thing binary pulse codes 
permit substantial transmission distortion of pulses within certain 
tolerable limits with negligible degradation of received signals. For 
another, regenerative repeaters can be used at intervals along a route to 
prevent accumulation of transmission distortion of pulses from various 
sources, so that virtually the entire allowable distortion can be permitted 
in each link or repeater section. 

The above desirable properties are secured in exchange for increased 
channel bandwidth, and can be used to full advantage in applications of 
binary pulse systems to such transmission media as radio and wave 
guides, where transmission is at such high frequencies that increased 
channel bandwidth does not entail increased attenuation. In wire cir- 
cuits, however, where baseband transmission is the more economical 


89] 





SY2 rHE BELL SYSTEM TECHNICAL JOURNAL, JULY 1957 


method, attenuation increases nearly in proportion to the square root of 
the channel bandwidth. For this reason, rather short repeater spacings 
may be required for binary pulse systems, so that for economical appli- 
cations to wire circuits it is imperative to have reliable regenerative 
repeaters of simple design. 

In their principle of operation regenerative repeaters are by nature 
more complicated than ordinary repeaters. In addition to providing gain 
to off-set attenuation in the transmission medium, as in ordinary re- 
peaters, they must also perform gating operations for sampling and 
regenerating the received pulse train. This, however, does not pre- 
clude the possibility that these operational principles can be implemented 
in repeater design by instrumentation that is simpler than required for 
ordinary repeaters. 

The possibility of simple instrumentation resides partly in the cir. 
cumstance that equalization circuitry for regenerative repeaters can be 
substantially simpler than for ordinary repeaters, owing to less exacting 
requirements on equalization. Furthermore, satisfactory performance 
in pulse regeneration can be achieved without very precise timing in 
sampling and regeneration of pulse trains. It is thus possible to secure 
nearly the same performance as for ideal regenerative repeaters by par- 
tial rather than complete exact retiming of pulse trains at each repeater. 
This facilitates simple gating arrangements for regeneration of pulses. 
Moreover, it permits a timing wave for control of gating operations to 
be derived from either the received or regenerating pulse trains with the 
aid of a simple resonant of circuit. 

The simplicity of instrumentation permitted by these considerations 
is exemplified in a self-timed regenerative repeater for baseband pulses 
invented by L. R. Wrathall of Bell Telephone Laboratories. The cir- 
cuitry of the repeater together with the results of tests on laboratory 
models are dealt with elsewhere! and not considered here. The purpose 
of this paper is an analysis of the timing principles underlving this type 
of repeater together with its regeneration characteristics as deter- 
mined by various basic design parameters, on the assumption of ideal 


implementation of the timing principles by appropriate instrumentation. 


In the Wrathall repeater “quantized feed-back” is employed as a means 
of reducing the effect of low-frequency cut-off in transformers. Since this 
is not an essential feature of self-timing repeaters and has no direct 
bearing on the timing principles, it is disregarded herein. 


1L. R. Wrathall, Transistorized Binary Pulse Regenerator, B.S.T.J., 35, pp 
1059-1084, Sept., 1956. 





SELF-TIMING REGENERATIVE REPEATERS 


I REGENERATION AND RETIMING 
1.0 General 


In an ideal regenerative repeater the received pulse train is sampled 
at proper fixed intervals, to determine whether a pulse is present. The 
regenerated pulses transmitted into the next repeater section are all 
of the same shape and amplitude, independent of the shape of the input 
pulses. Thus pulse distortion from noise and other system imperfec- 


tions is removed, provided the maximum distortion is held within proper 


limits. Errors in the form of pulses in ‘place of spaces, or conversely, are 
encountered when these limits are exceeded. In a repeater chain there 
will be cumulation of errors in proportion to the number of repeater 
sections in tandem. However, the rate of errors in each section and thus 
in the whole chain can be limited by a relatively small increase in the 
signal-to-noise ratio of each section as the number of repeaters in tandem 
is increased. This increase in signal-to-noise ratio with increasing length 
of the repeater chain is much less than with ordinary nonregenerative 
repeaters. For this reason regenerative rather than ordinary repeaters 
are desirable, though not essential for systems employing binary codes. 

An ideal regenerative repeater with the above features would entail 
rather complicated instrumentation for precise timing, sampling and 
pulse regeneration. With partial rather than complete exact retiming 
the repeaters can be simplified, in exchange for some sacrifice in per 
formance, as shown later. 


1.1 Regeneration Without Retiming 


It would be possible to have a repeater in which pulses would be re 
generated in amplitude and shape, but without retiming. Pulses would 
in this case be regenerated when the amplitude of the pulse exceeded a 
certain triggering level L. If the pulse shape is given by P(t), this would 
occur at a time ¢) such that 


P (to) L. 1.1 


This would permit simple instrumentation, since regenerated pulses 
would be triggered without separate sampling of the received pulse 
train. With this method, however, timing deviations in the regenerated 
pulses would result from transmission distortion of the received pulses 
by noise and other system imperfections. These timing deviations would 
cumulate in a repeater chain and cause a reduction in the tolerance of 
the repeaters to noise, such that the signal-to-noise ratio would have to 
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be increased with the number of repeaters in tandem in the same way 
as for ordinary repeaters. 


1.2 Regeneration with Complete Retiming 


With complete retiming, the instants of pulse regeneration would 
be controlled by a periodic retiming wave, R(t), with a fundamental 
period equal to the interval between pulses. The received pulse train 
would be sampled at instants when the retiming wave had a certain 
level L, . The sampling instants f) would thus be given by 


ioe) = be. (1.2) 


R(t) would satisfy this equation for fy = n7’ + AT, where T is the 
nominal interval between pulses, m is an integer and AT is a certain 
tolerable deviation from the desired sampling instants. Pulses would be 
regenerated provided P(t)) > L and would be omitted if P(t) < L. 

With this method the timing deviations in regenerated pulses would 
be limited to + AT’, regardless of the timing deviations in received pulses. 
There would be no cumulation of timing deviations in a repeater chain. 
However, the tolerance of the repeaters to noise would be somewhat 
reduced by the timing deviations + AT’. 


1.3 Regeneration with Partial Retiming 


Partial retiming is obtained by a combination of the above two 
methods, by triggering regenerated pulses without sampling at instants 
ty) determined by 


P (to) + R(to) = DL. (1.3) 


To permit regeneration without sampling and without a marked reduc- 


tion in the tolerance of the repeaters to noise, the timing wave R(t) 
must meet certain conditions illustrated in Fig. 1. One is that it must be 
a nearly periodic function as for complete retiming. The second condi- 
tion is that R(t) must be zero near the sampling points to obtain sub- 
stantially the same tolerance to noise in the presence of a pulse as in 
the absence of a pulse. A third condition is that R(t) must have sub- 
stantial negative values between sampling points in order that the 
repeater be rather insensitive to noise between sampling points, as with 
complete retiming. It will be recognized that, in general, the maximum 
value of R(t) need not necessarily be zero, as in the above illustration. 
It can be greater or smaller than zero, provided the triggering level is 
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modified accordingly. A maximum value of zero is, however, convenient 
from the standpoint of instrumentation. 

A limiting shape of retiming wave that would result in complete re- 
timing, but without the need for special sampling is also illustrated 
hig. 


1.4 Derivation of Timing Wave from Pulse Train 


As shown above, the retiming wave must be essentially periodic, with 
a fundamental frequency equal to the pulse repetition frequency f = 1/7, 
where 7 is the interval between pulses. The simplest form is a sinusoidal 
wave, which can be derived from the pulse train at repeaters with the 
aid of a narrow band-pass filter, such as a simple resonant circuit cen- 
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COMPLETE RETIMING 


Fig. 1 Principle of partial retiming method 
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tered on the pulse repetition frequency. This possibility resides in the 
circumstance that a random “on-off” pulse train can be resolved into 
two components. One is an infinite sequence of pulses of the same 
polarity and equal amplitude, the other a sequence of randomly positive 
and negative polarity. The response of a resonant circuit to the first 
component is a steady state sinusoidal wave of the pulse repetition fre- 
quency. The second component gives rise to random variations in ampli- 
tude and phase, which in principle can be limited to any desired extent 
by limiting the band of the resonant circuit and the deviation in the 
resonant frequency from the pulse repetition frequency. 

A principal feature of this method of “‘self-timing’’, aside from its 
simplicity, is that the timing wave becomes a slave of the pulse train. 
Thus, if there is a fixed delay in pulse regeneration at a repeater, the 
same delay is imparted to the timing wave derived from the pulse train 
at the next repeater. This prevents a cumulation of such fixed delays 
with respect to the timing wave, but not with respect to an absolute 
time scale; i.e., with respect to an ideal timing wave transmitted along 
the repeater chain and independent of the pulse train. 


1.5 Self-Timed Repeaters with Partial Retiming 


A timing wave derived from the pulse train with the aid of a resonant 
circuit can be used in conjunction with complete or partial retiming. 
With complete retiming, pulses could be regenerated at the zero points 
in the timing wave, and the effects of amplitude variations in the timing 


wave can thus be avoided. Timing deviations in the regenerated pulses 


would in this case depend only on phase deviations in the timing wave, 
caused partly by the component of randomly positive and negative 
polarity in the pulse train and partly by timing deviations in the pulse 
train from which the timing wave is derived. 

With partial retiming the situation is more complex. Timing devia- 
tions in regenerated pulses in this case depend not only on amplitude 
and phase variations in the timing wave, but also on the regeneration 


characteristics of the repeaters. 


1.6 Types of Timing Deviations 


In a regenerated pulse train there will be fixed and random timing 
deviations. Of the latter there are three types. One is the timing devia- 
tion taken in relation to an exact timing wave with a period 7 equal to 
the nominal pulse interval. The second is the timing deviation taken in 
relation to the timing wave derived from the pulse train, which in itself 
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will contain random deviations. The third type is random deviations in 
the interval of adjacent pulses. If the first type is held within tolerable 
limits, this will also be the case for the second and third types. For this 
reason only the first type is considered herein. 

II REGENERATION CHARACTERISTICS WITH PARTIAL RETIMING 

2.0 General 


With partial retiming, there will be timing deviations in the re- 


generated pulses as a result of timing deviations, amplitude variations 
and distortion by noise of both the received pulses and the timing wave. 





Fig 2 Reduction in tolerance to noise by displacement in timing wave 
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The conversion of these variations into timing deviations in the re- 


generated pulses depends on certain relationships between the pulse 


train and the timing wave, discussed in the following sections. 


2.1 Tolerance to Noise 


From Fig. 2 it can be seen that if the timing wave is displaced by 7o , 
the value of P(t) + R(t — 70) in the presence of a pulse exceeds the 
triggering level by a maximum amount 


[P(t) + R(t — 10) — L)max = (P(t) — L}. (2.1) 


It will be recognized that the right-hand side of this equation represents 
the tolerance to noise of negative amplitudes with instantaneous sam- 
pling at ¢ = 7», as in an ideal repeater with complete retiming. 

With partial retiming, the tolerance to noise will be less than the above 
maximum value. However, it will be greater than the average of P(t) + 
R(t — 7) — Lin the range where the latter difference is positive. Let it 
be assumed that it is smaller than the maximum by a factor k somewhat 
smaller than unity. The tolerance to noise with a displacement 79 in the 
timing wave is then smaller than without a displacement (i.e., 7 0) 
by the factor 


k{ P(r) a, L) P(r) 6s L 


* KPO) —L) PO)—L’ 


The tolerance to noise will thus be reduced in a way similar to that 
for an ideal repeater with complete retiming. The absolute tolerance to 
noise will be less than for a repeater with complete retiming by a factor 
/; somewhat smaller than unity, say in the order 0.8, corresponding to 
about 2 db 


2.2 Conversion of Timing Deviations 


With partial retiming, timing deviations in received pulses and in the 
timing wave are converted into smaller deviations in regenerated pulses. 

Let 7, be a time displacement in a received pulse and 7, in the timing 
wave, both in the positive direction. Pulses will then be regenerated 
at a time f)’ given by 


Pits — +) + RG L 


where the minus signs are used since this corresponds to a displacement 


of P and R in the positive direction. Subtracting (1.3) from (2.3), 


P(ts’ — t») — P(to) + R(t’ — 7) — R(to) = 0. (2.4) 
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By adding and subtracting P(to’) + R(t’) and rearranging terms, (2.4 
can also be written 
[P(to’) — P(to)| + [R(b’) — R(t) 

[P(ty’ P(to’ | + [R(b’) — R(t — 7 
‘or small values of 7, and 7, , such that 6, = fo’ — fo is sufficiently small, 


both sides of (2.8) can be represented in differential form as 


6,{P’(to) + R’(to)| T pl?" (ti + 7h’ (ty 2.60 


where P’(to) = dPo(t)/dt at ¢ fy , and RF’ is correspondingly defined 


—_ 


Mquation (2.9) can be written in the form 


where 


Pe 8 | 


With random uncorrelated displacements of rms values 7, and 7 
the rms value of 6, is 


6; (pr 7 ay 2.10 


Kquation (2.9) and (2.10) give the timing deviations in regenerated 
pulses in terms of the deviations Tp and 7, in the received pulses and 
in the timing wave. To limit timing deviations in the regenerated pulses, 
it is necessary to make p, and the product r,7, small. This will entail 
the use of a timing wave comparable in amplitude to that of the pulses, 
or greater, in conjunction with a small timing deviation 7, in the timing 
wave. 


2.3 Conversion of Amplitude Variations Into Timing Deviations 


With partial retiming there is a conversion of amplitude variations 
in the received pulses and in the timing wave into timing deviations in 
the regenerated pulses. 

Let the pulses have an amplitude variation a, and the timing wave a 
expressed as fractions of the normal values. Pulses will then be regen 


erated at a time &)’ given by 


l a,)P(to’) 4 ] y 2.11 
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Subtracting (1.3) trom (2.11 
[P(to’ P(ty)| + [R(ty’ ) - R(to)| apP (to’) a,P (to’). 


lor small values of a, and a, , such that 6, ‘ — fy is sufficiently small, 


the same procedure as in Section 2.2 gives 
(Pidy — Tae 
and 


P (to) — PR (to) 


), ; . \ : ; , : 
f P(t) + R'(to) P(t) + R(t) 


lor uncorrelated variations of rms amplitude a, and a, the correspond- 
ing rms timing deviation is 

r Ta @ (2.14) 

Iquations (2.12) and (2.14) give the timing deviations in regenerated 

pulses resulting from amplitude variations in the pulses and in the 


timing wave. 


2.4 Resultant Timing Deviations in Regenerated Pulses 


lor small variations in the pulses and in the timing wave as considered 

previously, the resultant timing deviation in a particular regenerated 
pulse is 

A . i (2.15) 

Considering a large number of pulses, the resultant rms timing devia- 

tion in terms of the rms deviation in the received pulses and in timing 


wave IS 


(2.16) 


\ 


—_ 


3 
A I 
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Ill ILLUSTRATIVE REGENERATION CHARACTERISTICS 
JO General 


In this section the general equations given in the preceding sections 
are applied to a particular case, in order to obtain specific expressions 
for the regeneration characteristics and illustrative curves, as an aid 
to further analysis. The particular case selected for illustration approxi- 
mates the conditions in experimental Wrathall repeaters, and may be 
regarded as an idealized model of such a repeater, in which certain effects 


to be discussed later are ignored. 


3.1 Pulse Shape 


It will be assumed that the pulses are transmitted at intervals 7 
and that the shape of the received pulses after equalization is given by: 


P(t) | 1 5 | 
2 n T. 
This is the familiar ‘‘raised cosine” type of pulse. With » | the pulse 
width is the maximum that can be tolerated without intersymbol inter- 
ference. With » ;, the amplitude of a pulse train at a point midway 
between two success pulses is equal to half the peak amplitude of a 
pulse. The latter assumption will be made here, for reasons discussed 
later. 


3.2 Retiming Wave 


The retiming wave is assumed to be given by 
R(t) 


This type of retiming wave can be obtained if a sinusoidal wave of the 
pulse repetition frequency f = 1,7 is applied to a resonant circuit to 
reduce distortion of the timing wave by noise. The resonant circuit 
would have a nominal resonant frequency f 1/T, but because of mis 
tuning it would actually be fo. The output of the resonant circuit after 
appropriate adjustment of amplitude would be of the form [Appendix I, 


equation (2)|: 


l 


| t 
P(t) 5 COs cos (2 a ), 


where y is the phase shift of the resonant circuit at the frequency /, 
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given by: 


f fi 
tan y a(t _ a (3.4) 


and Q is the loss constant of the resonant circuit. If the peaks of the 
wave given by (3.3) are held at zero potential, a retiming wave as given 
by (3.2) is obtained. This type of retiming wave can also be obtained by 
applying an infinite sequence of rectangular pulses of equal amplitudes 
with spacing 7 to a resonant circuit. 
3.3 Triggering Instants 

With a pulse shape and retiming wave as assumed above, the resul- 


tant wave is given by 


hi 
P(t) + R(t) |! + COs 


‘) 


T 
n 


This wave is shown in Fig. 3 for y 0 and +60°. For y +90° the 


retiming Wave disappears, so that the combined wave is P(t). 


T = 360° 





Fig. 3 Illustrative example of pulse shape and retiming wave 
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The triggering instants ¢) are obtained from the relation 


P(to) + R(t) re 3.6 


With complete retiming the optimum performance, with positive and 
negative noise amplitudes of equal probabilities, is obtained with a 
triggering level 3. With partial retiming, optimum performance is ob- 
tained with a somewhat lower triggering level, but this is of secondary 


importance in connection with the present analysis. For this reason 
L, = } is assumed, in which case the following equation is obtained for 
determination of fo : 














“0.3 
-90 all -60 -45 -30 = 1@) 15 30 45 60 
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Fig. 4 Triggering times versus phase shifts in timing wave 
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TABLE I — VALUES OF ¢)/7' FOR 7 3 AND 7 


0.198 0.156 
0.215 0.170 


FOR 


This equation is satisfied for the values of t)/7' given in Table I 


values of fo/7 are also shown in Fig. 4 as a function of y. 


Conversion Factors for Time Deviations 
The conversion factors defined by (2.8) become: 


r » et 
SIn 


D n T 


ar fo & , 
D 2n cos ¥ sin (2m 2 v); 


‘ T to R . 0 

D = sin - — + 2n cosy sin{ 27 — — v), 
n 1 | 

T has the values given previously as a function of y. 

? are given in Table II and 


where f 
lor various values of y, the factors for n 


in Fig. 5 
3.5 Conversion Factors for Amplitude Into Time Deviations 


The conversion factors defined by (2.13) become 


Tn 1 14 w ti 
—_4.= + cos- — |, 
. mae 77 


Tn | 
r D 


where D and 4)/7' are defined as before. 
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a an cae 
PHASE SHIFT 

Fig. 5 Conversion of timing deviations in received pulses and in timi 

wave into timing deviations in regenerated pulses, for pulse shapes and timing 

3. Timing deviations in regenerated pulses in relation to 


vy 
i 


waves shown in Fig 
timing deviation ¢, in pulses and ¢, in retiming wave is p,, 


lor various values of ¥ the factors for n ; are given in Table IT] 


and in Fig. 6. 


3.6 Correlated Amplitude and Time Deviations 


The amplitude and time deviations in the pulses are generally uncot 
related, but this does not always apply to the timing wave. In particular, 


if a deviation 7, in the timing wave is the result of a change in the phase 


7 


1? 


y, it will be accompanied by a given amplitude variation. A change 
phase by Ay is related to the corresponding time deviation 7 by 


TABLE II] VALUES OF pa/T -,/T FOR 7 


0.185 0.175 0.19 0.22 
0.035 0.055 0.072 0.106 
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With this change in phase, the factor cos y of (3.2) is modified to 
£ } 


cos (Wy + Ay) cos y cos Ay — siny sin Ay, 
or 
~ cos yy — =7, siny 
1 
where the approximation applies for small values of y. The amplitude 
variation resulting from the above change in phase is accordingly 
dy = —T; id sin y. 
T 
Considering both the time deviation 7, and the corresponding ampli- 
tude variation a,, the resultant time deviation in regenerated pulses is 
in accordance with (2.20) 


3.16) 























60 a ) > 3 45 60 
PHASE F ANGLE MING WAVE, W, IN EGREES 


Fig. 6 — Conversion of amplitude variations in received pulses and in timing 
wave into timing deviations in regenerated pulses, for pulse shapes and timing 
waves shown in Fig. 3. Timing deviations in regenerated pulses for amplitude 
variations a, and a, in received pulses and in timing wave is pa@, + rad 
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where 


-¥)| 


and D and to, T are defined as betore. 
The factor 8 indicates the time deviation in regenerated pulses in 


relation to the time deviation 7, in the timing wave which results from 
a phase shift Ay as given by (3.13). It may be regarded as a timing feed- 
back factor that is of interest in connection with timing from regenerated 
pulses as discussed later. The factor 8 is shown in Fig. 4 for 2 and 
n i? 


3.7 Reduction in Tolerance to Noise by Timing Deviations 


When the pulse shape is given by (3.1) and the timing wave is dis- 
placed by 7», the tolerance to noise is in accordance with (2.2) reduced 
by the factor 


lor a phase displacement y, 


3.21) 
2n 


lor » = 4, the factor uw and the corresponding reduction in the toler- 
ance to noise in db are as follows: 


0.866 0.766 
1.2 2.3 
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IV DERIVATION OF TIMING WAVE FROM PULSE TRAIN 


1.0 General 


. 


The retiming wave F(t) must have a fixed relation to the received 
pulses, with certain tolerable fixed and random deviations to be con- 
sidered later. Such a timing wave can be derived from the pulse train 
with the aid of a sufficiently narrow band-pass filter, the simplest form 
of which is a resonant circuit consisting of a coil and capacitor in series 
or in parallel. 

A train of rectangular ‘‘on-off” pulses is shown in Fig. 7 as it would 
appear at the output of a regenerative repeater and at the input of the 
next repeater, (dotted) with uniform intervals 7 between sampling 
points. 

As indicated in Fig. 7, the pulse train can be regarded as being made up 
of two components. One of these is an infinite sequence of pulses of one 
polarity, the other an infinite sequence of randomly positive and negative 
polarity. 

It will be recognized that the first of the above components at the out 
put has a fundamental frequency equal to the pulse repetition fre- 
quency, | | 7, and the forced response of a resonant circuit to this 
component will be the pulse repetition frequency, regardless of any im 
perfections in tuning. In order that this frequency be present in the re- 
ceived pulse train, it is necessary that the spectrum of the received pulses 
extend beyond the pulse repetition frequency, so that there will be a 
ripple in a long sequence of received pulses of one polarity, as indicated 
in the illustration. 

The second random component of the' pulse train will have a fre 
quency spectrum that is nearly uniform over the band of the tuned 
circuit, and which will vary in amplitude depending on the composition 
of the pulse train. The response of the tuned circuit to this component 
is thus rather complex, and must be treated on an approximate statistical 
basis. It will consist of an almost periodic wave with random amplitude 
and phase modulation, and with mean frequency equal to the resonant 
lTrequency. 

Owing to the presence of the second component, there will be a 
variation with time in the amplitude and phase of the response of a mis- 
tuned resonant circuit, and resultant deviations in timing. The re- 
generated pulses will thus not be uniformly spaced, but will in general 
have random deviations from the desired exact positions. Such deviations 
can be created by superposing on a train with uniform spacing a random 


dipulse train, as indicated in Fig. 7. The resonant circuit response to this 
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dipulse train would be expected to be smaller than to the random ampli 
tude component of the pulse train. It may be regarded as a third compo 
nent representing a second order effect resulting from the second co 
ponent. 


mM 


In the Appendix, this method of superposition has been used as a basis 
of an analysis of a resonant circuit response to a random binary pulse 
train. This problem has also been dealt with by somewhat different 
methods in prior unpublished work by W. R. Bennett and J. R. Pierce, 
both of Bell Telephone Laboratories. 

In this analysis it is assumed that the regenerated pulses are of suffi- 
ciently short duration to be regarded as impulses. The response of the 
resonant circuit to the second and third components above, when taken 


in relation to that for the first component will, however, remain very 
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nearly the same for other pulse shapes, provided the frequency spectrum 
of the pulses can be regarded as approximately constant over the im- 
portant portion of the band of the resonant circuit. This approximation 
is legitimate for resonant circuits with a loss constant Q and pulse shapes 


at the input of repeaters as considered here. 


).1 Resonant Circuit Response to Steady State Component 


The first component consists of an infinite sequence of impulses of 
amplitude 3 and all of the same polarity, at intervals 7. This sequence 
has a fundamental frequency f 1 7. When it impinges on a resonant 
circuit with resonant frequency fi f Af and loss constant Q, the 


response is of the form 


A,(t cos Y cos (wt — w), (4.] 
and 
ou wr Af 

tany = QUS/fo — fo/f) = 2Q f ; (4,2) 
The response is thus a steady state sinusoidal wave of frequency { 
displaced from the fundamental component of the input wave by the 
phase shift y and reduced in amplitude by cos y. This is the phase 
shift and amplitude reduction of the resonant circuit at the frequency f 
when the resonant frequency is fi 


1.2 Resonant Circuit Response to Random Signal Component 


The second component consists of an infinite random sequence of im- 


pulses of amplitude +3, at intervals 7. The response of the resonant 


circuit to this component will be a randomly fluctuating wave A,(¢) of 
mean value 0. The maximum positive amplitude is obtained when all 
impulses of the second component are positive and is A,(0) A, . The 
maximum negative amplitude is A,(t)- A, . Owing to the presence 
of this component the total output of the resonant circuit A, + A,(f) 
can thus fluctuate between the limits 0 and 24, , but the actual fluctua- 
tions of significant probability will be smaller. 

The above fluctuations can be resolved into a component in phase 
with the steady state response given by (4.1) and another component at 
quadrature with the steady state timing wave. The rms values of these 
components taken in relation to the amplitude of the steady state wave 
are 
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(3) 
1() cos Y 


These relations apply for small values of y and for 7/Q < | 
The resultant rms amplitude variation in the timing wave is a 
as given by (4.3). 
The rms phase error Z, resulting from the quadrature component a,” 
is given by 


tan @, $5 


The corresponding rms time deviation is (7) 2x 


QO! 


ji (2 ) ~ ly 16 
2r 1() COS y ‘ 


With regard to the probability of exceeding the above rms values by 


various factors the normal law can probably be invoked with reasonable 
accuracy. As mentioned before, the maximum possible amplitudes are 
A,(t) +A, which would correspond to a peak factor (2Q 7) ~. With 
(2 = 100, the factor is about 8, while with Q 1000 it is about 25 
Based on the normal law the probability of exceeding the rms value 
by a factor of 4 is about 5 X 10°”, and by a factor of 5, about 10 °. The 
normal law would be expected to apply, since the limiting peak values 
are substantially greater than the peak values expected with significant 
probabilities. 


j.3 Resonant Circuit Response to Pulse Displacements 


Because of the random components given by (4.3) and (4.4), the 


timing wave will contain small random amplitude and phase deviations 
from a sinusoidal wave represented by (4.1). This will result in small 
random deviations in the positions of regenerated pulses triggered from 
the timing wave, which is represented by the third component shown in 
lig. 7. When the rms deviation in the pulse positions is 6, there will be 
an additional random quadrature component in the timing wave which 
when taken in relation to the steady state component, is given by 


" ie 
d { 4 4 | Wd ( ) 
Q) 


The corresponding rms phase deviation is given by 


” 
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The resultant rms time deviation is (7'/22)g;5 or 


05 da, (4.9) 
and 


a = (r/Q)'”. (4.10) 


The above factor @ applies to a single resonant circuit. When the rms 
timing deviations represented by (4.9) are present in the regenerated 
pulse train, the rms deviation at the output of the second resonant 
circuit is 


where 
a) Q. 
With n resonant circuits in tandem, 
05,, §aia2a 


The factors a, are given by 


Ql 


(4.14) 


(4.15) 


a’ ( when aa (4.16) 
rn 


The factors q@; for 7 2 2 represent the reduction in timing deviations 


resulting from the reduction in bandwidth as resonant circuits are 
added in tandem. If resonant circuits with a narrow flat pass-band were 
used, the bandwidth of any number of resonant circuits in tandem 
would be the same as for a single resonant circuit. In this case q@ 


Qa Qa, l. 
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1.4 Deviations in Timing Wave 


The timing wave derived from an ‘‘on-off”’ pulse train with the aid of 
a resonant circuit will in accordance with the expressions given in the 
previous sections contain three types of amplitude and timing devia 
tions. 

The first type is a fixed amplitude reduction by a factor a» and a fixed 
time deviation zo given by 


where y is given by (4.2). 


The second type is a random amplitude and time deviation resulting 


from the random amplitude component of the pulse train, which have 


rms values 
a 
a = 
2Q) 


/ 


T (2 l 
5, = ( 2 ) a 120 
=I \4Q/) cosy 


The third type is a random amplitude and time deviation resulting 
; | 


from random timing deviations 7, 6 in the pulse train. The amplitude 


variation can be disregarded and the rms time deviation is 


) 

(5 
The total rms amplitude variation is accordingly given by (4.19 

The total rms timing deviation obtained by combining (4.20) and (4.21 


Is 


- az, 


The expressions for 6, and 7, are the quantities appearing in (2.20 


lol 
A,, the total rms timing deviation in regenerated pulses resulting from 


random amplitude and timing deviations in the timing wave. 
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V SELF-TIMED REPEATERS WITH PARTIAL RETIMING 
5.0 General 


As shown in the preceding section, timing for pulse regeneration can 
be derived from the pulse trains, with certain random phase and ampli 
tude variations in the timing wave that can be reduced by increasing the 
loss constant Q of the resonant circuit. This method of ‘‘self-timing”’ 
can be combined with partial retiming, and the regeneration charac 
teristics of this type of repeater will be discussed in the following sections. 

For purposes of numerical illustration, the same type of pulse shape 
and timing wave will be assumed as in the previous numerical illustration 
in Section III. This pulse shape and timing wave closely approximates 
those in experimental Wrathall repeaters, in which timing is derived 
from the regenerated pulse train. In the following discussion timing 
from the received pulse train will also be considered. 


5.1 Timing from Received Pulse Train 


It will be assumed that the timing wave is derived from the received 
pulse train with the aid of a resonant circuit and that random timing 
deviations are absent. The response of the resonant circuit is then a 
sinusoidal wave as given by (4.1). From this wave it is possible to obtain 


a retiming wave of the form 


, 
R(t) —cosy}|1— cos (2 T -y) (5.1) 


This can be accomplished by holding the peaks of the timing wave from 
the resonant circuit at zero potential with a diode. This is the form of 
retiming wave previously considered in Section IIT, in conjunction with a 
pulse shape given by (3.1). 

As shown in Section 3.7, the tolerance to noise will vary with the 
phase shift y of the resonant circuit, in accordance with (3.21). If a 
reduction in the tolerance to noise of about 2 db is allowed, the maxi- 
mum permissible phase shift would be about y = 1 radian (57.6°). On 
this basis the maximum permissible deviation Af,,.x in the resonant fre- 
quency from the pulse repetition frequency f as obtained from (4.2) 
with y | radian becomes 


Afmax _ tan yp 1.58 
f 20Q) 2Q° 


lor various values of Q in the range that can be realized by simple 
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resonant circuits, the permissible deviations are as follows: 


1 Pe | 0.08 0.030 0.016 0.008 0.004 


This assumes that there are no random timing deviations and that the 
tolerance to noise is reduced by not more than 2 db. 


5.2 Timing from Regenerated Pulse Train 


It will again be assumed that there are no random timing deviations. 
Without a phase shift in the resonant circuit, let the regenerated pulses 
be triggered at a time ¢) . When there is a phase shift y’, the pulses will 
be triggered at a time fo’. The timing wave derived from the regenerated 
pulses will then have a time shift 


A 


This time shift will cause pulses to be regenerated with a time shift 
B’A, which must equal fo’ — ft . Accordingly, 


and 


pe 


B 


T BW 


2r | 


With timing from the received pulse train with a phase shift y in the 
resonant circuit, the following relation applies: 
T _ 
2r ss . 
If to’ — ty) is to be the same in both cases, so that the timing wave and 
tolerance to noise is the same, the following relation must exist between 
the phase shifts in the resonant circuit: 


In this expression, 8 and @’ are the factors shown in Fig. 4. It will be 
recognized from (5.5) that the smallest permissible phase shifts are ob 
tained for large values of 6’. From Fig. 4, it is seen that the largest 
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values of 8 are for phase shifts between 0 and —60°. For 7 
and for 7 1j@=> O09. 

For 7 + and 7 1 the tolerable maximum phase shifts y’ in the 
resonant circuit with timing from the regenerated pulse train, in rela- 
tion to the maximum tolerable y with timing from the input, are 


vy’ ~0.3y for »n = 3, (5.6) 


and 
yy’ =0.ly for n= 1. 


Although greater phase shifts can be tolerated when y is positive, and @’ 
is smaller than above, the requirements on the resonant circuit must be 
based on the worst condition that can be encountered, as above. 

From (5.6) it follows that for » = 1 the requirements on the per- 
missible phase shift in the resonant circuit are much more severe than 
for n = 2. For this reason the latter value of 7 is decidedly preferable 
for the particular case in which the peak amplitudes of the pulse train 
and the timing waves are equal, as assumed here. A value 7 = ? is also 
desirable from the standpoints of avoiding intersymbol interference 
between adjacent pulses at the triggering instants, to permit the timing 
wave to be derived from the pulse train and to permit self-starting of 
the repeaters, as discussed later. 

In accordance with (5.6) the maximum tolerable frequency deviation 
for 7 + will be less than with timing from the received pulse train by 
a factor of about 0.3. The maximum permissible frequency deviation for 
a phase shift of about one radian in the timing wave and 0.3 radian in 
the resonant circuit, will accordingly be about as follows: 


Af, f 0.025 0.009 0.005 0.0025 0.0012 


lor a repeater with complete rather than partial retiming, the factor 


8 would be unity, and timing from the regenerated pulse train would 


not be possible 


5.3 Random Timing Deviations 


In combining random timing deviations from various sources at a 
particular repeater, it will be assumed that there is no correlation be- 
tween the various deviations, so that they will combine on a root-sum- 


square basis. 
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In accordance with (2.21) the rms timing deviation at the output is 
then: 


where in accordance with (4.13) and (4.16 


9) 


) 

q : 

, ( ” ) 2 P 

. 2r \4Q COs y 


When (5.9) 1s inserted in (5.7 


A = (pr + or, jt, - “= rr 9d o.12 


= — I 


This expression gives the rms timing deviation at the output in terms of 
the rms deviation 7, at the input and the various repeater parameters 


With timing from the output, rather than the input as assumed above, 
7, is replaced by Q in (5.9), and the following relation is obtained: 


A‘(1 - ar, gees my Or + To Gr . ots 


In the above expressions p, & 0.15, r- = O04 and @ = 0.03 (Q 
100). The term a@r, can thus be neglected in comparison with p, in 
(5.12) and in comparison with 1 in (5.13). 

The following expression is thus obtained with timing from either the 
input or the output: 


A = (pr tp 


Ap 


5.4 Magnitude of Random Timing Deviations 


The first two terms of (5.14) represents the rms timing deviations in 
the regenerated pulses resulting from timing deviations and amplitude 
variations in the received pulses. The last two terms represent the timing 
deviations resulting from timing deviations and amplitude variations in 
the timing wave. The conversion factors p, , Pa, 7; and rq are discussed 
in Section II and representative values given in Figs. 5 and 6. The values 
of a, and §, are obtained from (5.8). 
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TABLE IV Rams DEVIATIONS FROM TIMING WAVE 
DISTORTION FOR Q 100 


0.011 0.005 0 0.006 0.015 
0.006 0.007 0.009 0.014 0.024 
0.0126 0.009 0.009 0.015 0.028 


1.5 3.2 3.2 3.4 10 


In Table IV are given the values of the two last terms in (5.14), which 
represents the rms deviations Q, owing to random deviations in the 
timing wave. The results are given for the particular case in which 
(2 = 100, and for other values of Q are inversely proportional to Q'”. 
The table shows the deviations as a fraction of the interval 7’ between 
pulses, and also as the corresponding rms phase deviation ¢ 

In Table V are given the values of the first two terms in (5.14), which 
represents the rms deviation J, in the regenerated pulses resulting from 
random amplitude and timing deviation in the received pulses. In binary 
systems it is customary to limit the rms pulse distortion to a, ro; 
corresponding to ;'5 the peak amplitude of the received pulses, or } the 
triggering level (17 db signal-to-noise ratio). The corresponding rms 


phase deviation would be about 75 radian, corresponding to an rms 


deviation 7, in the pulses of 0.016 the pulse spacing, or 7,/7 = 0.016. 


The total rms timing deviation obtained from (5.14) and the correspond- 
ing rms phase deviation are given in Table VI. 


TABLE V Rus DeEvIATIONS RESULTING FROM PULSE DISTORTION 


0.019 0.018 0.019 ().022 0.032 
0.010 0.007 0.005 0.005 0.008 
0.021 0.020 0.020 0.023 0.033 


7.9 fier. ice 8.2 12 


TABLE VI ToraL Rus DEVIATIONS FROM TIMING WAVE 


AND PULSE DISTORTION 


0) .025 0.022 0.022 0.028 0.043 


q S S 10 16 
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The probability that random phase deviations will exceed the above 


rms values by a factor of more than 4 is small enough to be ignored. On 


this basis the sum of the fixed and random deviations would be limited 
to about 70°, if the fixed phase shift Y is less than +30°. With this re 
quirement on the fixed phase shift for satisfactory performance, the 


values of Afi... would be about half as great as previously given in Sec 
tions 5.1 and 5.2, for a single repeater as considered here 


VI REPEATER CHAINS 
6.0 General 


In the previous section, a single self-timed repeater was considered, 
from the standpoint of fixed and random timing deviations, as deter 
mined by various repeater design parameters. In a repeater chain there 
will be some cumulation of random timing deviations as the number of 
repeaters in tandem is increased, and a resultant reduction in the 
tolerance to noise of repeaters toward the end of the chain. Exact evalua 
tion of such cumulation is rendered difficult by the circumstance that 
timing deviations from various sources may not follow the same law of 
combination along the repeater chain. In the following, expressions are 
given based both on root-sum-square and direct addition of random 


timing deviations, which can be regarded as lower and upper limits 


6.1 Combination of Random Timing Deviations 


To determine the rms value of random timing deviations at the end 
of a repeater chain, it is necessary to combine random deviations from 
various repeaters. Random deviations from various sources at a repeatet 
do not necessarily follow the same law of cumulation along a repeater 
chain. Since there is no correlation between timing deviations caused by 
noise in various repeater sections, these can be combined on a root-sum 
square basis. This, however, may hot be appropriate as regards the 
combination of timing deviations resulting from imperfections in the 
timing wave. Thus, with perfect tuning of all resonant circuits, the 
timing Waves at various repeaters would have virtually identical ampli 
tude variations, but no phase deviations. While in this case there would 
be complete correlation between the timing wave variations at the 
repeaters, it does not follow that the resultant timing deviations should 
be combined directly rather than on a root-sum-square basis along the 
repeater chain. The timing deviations at the end of a chain of N re 
peaters resulting from amplitude variations in the timing wave of the 
first repeater will be modified by N intermediate resonant circuits. Those 
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resulting from amplitude variations at subsequent repeaters will be modi- 
fied by N-1, N-2 etc. intermediate resonant circuits. The situation is 
similar to that of applying identical noise waves at the input of each of 
V resonant circuits in tandem. At the output the V noise waves will have 
different shapes owing to restriction of the band and increasing phase 
distortion as the number of resonant circuits in tandem increases. For 
this reason combination on a root-sum-square basis appears justified also 
in this case, particularly with various degrees of mistuning of the res- 
onant circuits, so that the amplitude variations in the timing waves 
will differ in phase among repeaters. 


6.2 Propagation of Timing Deviations 


To determine the cumulation of timing deviations along a repeater 
chain, it is convenient to first consider a single repeater as a source of 
timing deviations, and to determine the propagation of these timing 
deviations along a repeater chain. In the following, y, will designate the 
rms propagation factor for n repeaters in tandem; 1.e., the factor by 
which the rms timing deviations at the end of a chain of n repeaters is 
smaller than at the first repeater, with timing deviations originating at 
the first repeater only. 

Let the rms timing deviation at the output of the first repeater as 
given by (5.14) for convenience be taken as unity. At the output of the 
second repeater the squared rms timing deviation is then reduced by the 
factor 


Y2 Dr +TQrr, Q) a. (6.1) 


As indicated symbolically in Fig. 8, the first term represents the reduc- 
tion owing to partial retiming. The second term is the additional devia- 


























Propagation of random timing deviations along repeater chain 
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tion resulting from the effect on the timing wave of unit rms deviation 
in the received pulse train at the second repeater 

At the output of the third repeater, the squared rms deviation is 
smaller than at the output of the first repeater by the factor 


T 2a; ~P; ry 


As indicated in Fig. 8, the first term in (6.2) represents the reduction 


owing to partial retiming of the received pulse train at the third repeatet 
The second term, (prair,) , Is the additional rms deviation resulting from 
the effect on the timing wave at the third repeater of an rms deviation 
Pr in the received pulse train. The third term (air,;-aor,) is the additional 
deviation caused by the effect on the timing wave of an rms deviation 
air, in the received pulse train. The factor aor, represents the modifica 
tion in the rms deviation ay, by a second resonant circuit, with a» de 
fined as in Section 4.3 

At the output of the fourth repeater, the rms timing deviation is 
reduced by the following factor, obtained in the same manner: 

v4 Dr SQ) Dr lr - 302102 DrTr 4 Os 1; « 6.4 

At the output of repeater n, the squared rms timing deviation 1s 

smaller than at the output of the first repeater by the propagation 


factor 


where p, and r, are defined as in Section 2.2, and ay , 


Section 4.3. 
In the above formulation the rms deviation at the output of the first 
repeater was assumed given by (5.14), which is an approximation ot 


(5.12) in which the term ar,7, was neglected. This term will have a 


different propagation factor p, , the expression for which differs from 


that for y, as given by (6.5) in that the subscripts of the factors @; are 


raised by one unit. Thus, 
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(ri 1) » 2 
p p T , Pr ly Q 
l (6.6) 
2 ! 
1 r I a» a @ 
The rms deviation at the output of repeater n thus becomes 
fa (ae. 7 ly T Arr Ty Pp». (6.7) 


In the case of repeaters with partial retiming the last term in (6.7) can 
be neglected, in which case the cumulation of timing deviation will be 
virtually the same when the timing wave Is derived from the regenerated 
as when it is derived from the received pulse train. 

The above expressions apply for resonant circuits consisting of a coil 
and capacitor which have a gradual cut-off. If resonant circuits with a 
flat pass-band and sharp cut-offs were used, as 3 a, and (6.5) can 
be simplified to 


y (1 Qi) Pp: tT Q\ (pr a Ja adil (6.8) 
6.3 Cumulation of Timing Deviations 


The cumulation of random timing deviations from various repeaters 
ina chain can be determined from the propagation constant given above 
for any prescribed law of combination of timing deviations from various 
repeaters. When equal rms deviations are contributed by each of 
repeaters, and they are combined on a root-sum-square basis, the rms 
deviation at the end of a repeater chain is greater than for a single 
repeater by the cumulation factor 


( (> Y ) (6.9 
\ n=] 
\n upper limit to C is obtained by taking a: Qs Q, 1 in (6.5) 
in which case y, is given by (6.8); (6.9 then becomes for .V s 
l 

| 

c. F Q) ) a ~ = (6.10) 
| l—pD l — 0; i 

| 

= : (6.11) 
| D; 

where the terms in a; have been neglected in (6.11), since ay a <l, 


about 0.03 for Q 100. 
From Fig. 5 it will be seen that when y < +60°, p, < 0.6. Hence 
( < 1.25. Cumulation of random timing deviations can thus for practical 
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purposes be disregarded, with root-sum-square combination as assumed 
above. The value of C obtained from (6.11) will differ from that obtained 
from (6.9) when 7, is given by (6.5), by a small fraction of one per cent 

Although root-sum-square combination appears justified for reasons 
given before, it is of interest to determine an upper limit to the cumula- 
tion based on direct addition of random timing deviations. The maxi 
mum cumulation factor thus obtained is 


Cocsig ht hy 6.12 


Employing (6.8) for y, and neglecting the terms in q; , the upper limit 


to the cumulation factor for N = x becomes 


] 
l— pe 


With p, < 0.6 fory < +60°, Cx < 2.5. 


x 


If the above maximum cumulation factor is applied to random timing 


deviations resulting from amplitude variations in the timing wave, as 
given in Table IV of Section 5.4, the resultant rms phase deviation at 
the end of a long repeater chain could be as great as 25°, rather than 10 
for a single repeater, when ¥ 60° and Q 100. To attain satisfactory 
performance it would in this case be necessary to limit the maximum 
fixed phase shift to substantially less than +60°, which would entail 
greater frequency precision than indicated in Sections 5.1 and 5.2 

fy < +15°, p, < O40 and Cissy < 1.7. In this case the rms phase 
deviation as given in Table I for a single repeater is ¢, & 4°, and the 
rms phase deviation in a long repeater chain would be less than 7°. In 
a long repeater chain the rms phase deviation resulting from pulse dis 
tortion would be greater than given in Table II by an rms cumulation 
factor C 1.08 for p, 0.4, and would thus be about 8° when y¥ 
+ 15°. The total rms phase deviation would thus be about (77> + 8°)'° = 
11°. Random phase deviations exceeding t times the latter value, o1 
about 45°, would be rather unlikely. The sum of the fixed and random 
phase deviations would thus be limited to about 60°, so that satisfactory 
performance would be expected when the fixed phase deviation is 
limited to about +15°. 

With the approximations for y, employed above, the rms cumulation 


factor for a chain of V repeaters as obtained from (6.9) is less than for 
N x2 by the factor (1 rg )* & 0.99 for D: 0.5 and NV 3. The 
maximum cumulation factor obtained from (6.12) is less than for \ Z 


by the factor 1 = pr =~ 0.99 for N 6. Thus, cumulation of random 
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timing deviations is virtually completed in a chain of 3 to 6 repeaters, 
so that for experimental determinations of the degree of cumulation it 
suffices to operate a few repeaters in tandem 


6.4 Repeaters with Complete Retiming 


In the particular case of complete retiming, p, = 0 and r, 
(6.5) and (6.6) so that 


(6.14) 
9 (6.15) 


For n > 1, approximation (4.16) can be employed, in which case 


( 7 ) | ( ; ) | | : 
Yn =al— : Pn = cit : (6.16) 
rn 


rn 
In this case (5.14) simplifies to 


a, rn 42 


SINCE Pa a Foo QO, p; 0 and r, 


Hence (6.7) becomes 


QA, 


With approximations (6.16), 


ay (6, (6.19) 


\t the output of the first repeater, 


Ai 


kor 7 1 the squared propagation tactor 1s accordingly 


i tae of Sw 
sat = gf BAH (1) 
0, T aTy Thi 


The squared rms cumulation factor for N > 2 repeaters becomes 


oN 


1 \l 
. (=) dn 
T eT; T 


In the particular case of perfect tuning of all resonant circuits 4, 
and 
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os 
Tr 
LV 
1 +( 
T 
¥ eS) 4 
= ‘ 6.24 
T 


The last expression gives the factor by which the rms timing devia 


tion at the output ot repeater N is greater than at the output ol the first 


repeater. The rms deviation at the output of the first repeater is greatet 


than at the input by the factor @. The rms deviation at the output of 
repeater V is thus greater than at the input of the first repeater by the 
factor, 
) ’ 
a F 
us 

kor this particular case (6 0) expressions equivalent to those 
above have been derived in unpublished work by H. Ek. Rowe of Bell 
Telephone Laboratories. 

In accordance with (6.22) and (6.24) the cumulation of random timing 
deviations increases indefinitely with V when retiming is complete. The 
cumulation factor as given by (6.24) is in fact the same as would be ob 
tained if a timing wave were transmitted on a separate pair, with a 
resonant circuit at each repeater to limit noise and with amplification of 
the timing wave at each repeater to obtain the same amplitude of the 
timing wave as when it is derived from the pulse train. With partial 
retiming cumulation is limited, for the reason that there is partial re 
generation of both the pulse train and the timing wave. 

Although with complete retiming the cumulation factor increases 
indefinitely with V, this is of but little practical significance, because of 
the slow rate of cumulation. At the output of a chain of NV repeaters an 
rms deviation approximately equal to that at the input of the first re 
peater could be tolerated, in which case (; & 1. On this basis the pet 
missible number of repeaters would be 


Ve T J T (=) 
t a | T 


~ 800 when Q LOO 


This assumes exact tuning of all resonant circuits. With mistuning of 
the resonant circuits, the permissible number of repeaters in tandem 


for a specified rms deviation at the output of the final repeater can be 
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determined with the aid of the cumulation factor given by (6.22). For 


example, if the rms deviation at the output of repeater NV is assumed the 


same as at the input of the first repeater, the permissible number of 


repeaters in tandem is less than given by (6.26) by the factor |(1 
1 + m)|, m 5, 7, . When the fixed phase shift is 30°, m = 
NV = 300 


6.5 Self-Starting of Se lf-Timed Re peaters 


With self-timing it is necessary that repeaters be self-starting if the 
timing wave should be absent for any reason. If each repeater is self- 
starting, this will also be the case for a repeater chain, since starting will 
be progressive along the chain. Initially, before the timing wave has 
reached the appropriate amplitude at all repeaters, there will be a high 
rate ol digital errors 

With timing from the received pulse train, the resonant circuit will be 
excited by every pulse and the timing wave will reach its normal amphi- 


tude in about n & Q pulses. With timing from the regenerated pulse 











hig. 9 Progression of repeater starting in absence of timing wave when 
timing is derived from regenerated pulse train 
Triggering points with timing wave absent. Noise prevents triggering at cer 
tain points, n. Timing wave reaches fraction of normal value, Ry 
Triggering points with timing wave 2. Timing wave increases to normal am 
plitude R 
@ Triggering points with normal timing wave 
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train the resonant circuit will not be excited by every pulse, unless the 


shape of the received pulses is such that there are virtually no overlaps 


between pulses so that the triggering level will be penetrated by each 
pulse. 

With a pulse shape as assumed in the previous analysis, the amplitude 
of a pulse train midway between pulses is half the peak amplitude of the 
pulses, as indicated in Fig. 9. In the presence of noise, triggering will in 
this case occur on the average for every second pulse, as indicated in the 
above figure. If it is assumed that the resonant circuit has the maximum 
permissible phase shift of about 20° allowed with timing from the output, 
the amplitude of the timing wave with excitation from every pulse will 
be virtually equal to the peak pulse amplitude. With excitation from 
half the pulses, the amplitude of the timing wave will rapidly reach half 
the peak amplitude of the pulses. When this initial timing wave is com 
bined with the pulse train, triggering will occur for virtually all pulses, 
as indicated in Fig. 9. It will thus reach its normal value. If the phase 
shift is greater than 20° as assumed above, say 60°, the initial amplitude 
of the timing wave will be } the peak pulse amplitude. Combination of 
this initial timing wave with the pulse train will increase the number of 
pulses exciting the resonant circuit, which in turn increases the amplitude 
of the timing waves, ete. 

Self-starting with a pulse shape as assumed in this analysis is thus 


insured. 


VIL SUMMARY 


In self-timing regenerative repeaters as considered here, a timing 
wave is derived from either the received or regenerated pulse train with 
the aid of a simple resonant circuit tuned to the pulse repetition Tre 
quency. This timing wave is combined linearly with received pulse trains 
as indicated in Fig. 1, and pulses are regenerated when the combined 
Wave penetrates a certain triggering level. 

It is concluded that if these timing principles are implemented by 
appropriate repeater instrumentation, a performance can be realized 
that approaches that of ideal regenerative repeaters. To this end it is 
necessary to meet certain requirements with regard to the loss constant 
(J of the resonant circuit, its frequency precision, the shape of received 
pulses and the amplitude of the timing wave in relation to that of re 
ceived pulses. 

Equalization of each repeater section should preferably be such that 
the received pulses have a shape and duration in relation to the pulse 
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interval as indicated in Fig. 3, and the peak amplitude of the timing 
wave should be about equal to that of the received pulses. Under these 
conditions the pulse repetition frequency will be present in the received 
pulse train in sufficient amplitude to permit derivation of the timing 
wave from the received pulse train, and to permit rapid self-starting in 
the absence of a timing wave if it is derived from the regenerated pulse 
train. 


~ 


A loss constant of the resonant circuit @ = 100 appears desirable. 
This value is sufficiently low to be readily realized with simple resonant 
circuits consisting of a coil and capacitor in series or parallel, without 
unduly severe requirements on its frequency precision. It is also ade- 
quately high from the standpoint of avoiding excessive random timing 
deviations in regenerated pulses from amplitude and phase deviations 
in the timing wave. 

The tolerable deviation in the resonant frequency from the pulse 
repetition frequency with Q 100 is about 0.2 per cent when the 
timing wave is derived from the received pulse train, and about 0.06 
per cent when it is derived trom the regenerated pulse train. These 
frequency precisions correspond to a maximum fixed phase shift of 15° 
in the timing wave, and allow tor the possibility that random timing 
deviations resulting from amplitude variations in the timing wave may 
cumulate directly along a repeater chain, rather than on a root-sum- 
square basis. With root-sum-square cumulation of timing deviations 
from all sources, the frequency deviations could be about twice as 
great. 

When the above requirements are met the reduction in the tolerance 
to noise owing to timing deviations in a repeater chain is limited to 
about 2 db. If the requirements on frequency precision of the resonant 
circuit are met, substantial degradation or improvement in performance 
would not be expected as a result of moderate changes in the other 


design parameters. 
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APPENDIX 
IX RESONANT CIRCUIT RESPONSE TO RANDOM BINARY PULSE TRAINS 
/ (re neral 


In the following analysis of the response of a resonant circuit to a 
binary ‘‘on-off’’ pulse train, the pulses are assumed of sufficiently short 
duration to be regarded as impulses. This is a legitimate approximation 
when the duration does not exceed about half the interval between 
pulses. 

The pulse train is regarded as made up of three components, as indi 
cated in Fig. 10. The first is a systematic component consisting of pulses 
of amplitude 3. This component gives rise to a steady state response at 
the fundamental frequency of the pulse sequence. The second com- 
ponent consists of pulses of amplitude +3, with random + polarity. 


This component gives rise to a random component in the resonant cir 
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Fig. 10 — Components of random binary on-off pulse train. A Transmitted 
“on-off”’ pulses. B Steady state pulse train of fundamental frequency f = 1 
C Random pulse train with zero mean value. D Random pulse train with 
displacements +6. F “On-off”? pulse train with displacements +6 from ay 
erage pulse interval 7 
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cult response; Le., a fluctuation about the steady state value derived 


from the first component. 
The third component consists of a train of dipulses. Each dipulse 
consists of a pair of pulses of amplitude 1 and 


ce displaced by an 
interval 


+6. The response of the resonant circuit to this component 


gives the effect of random displacements +6 in the original ‘‘on-off” 


pulse train. 


2 Impe dance of Resonant Circuit 


The impedance of a resonant circuit consisting of R, L and C 


in} 
parallel is 


Z (iw) Z(iw)e 


7 /: R 

Z (iw) = ats R cos y, 
1] T ("(w Wo — Wo w)*|1/ 

tan y O(w Wo 


where 


Q) wo RC Loss constant, 


(1 LC)’ ? Resonant frequency. 


The above expressions also apply for the admittance of a resonant 
circuit consisting of R, L and C in series, except that in this case () 


vy L R 


by Impulse Re SPONSE of Resonant Circuit 


When a rectangular pulse of unit amplitude and sufficiently short 
duration 6) is applied to a resonant circuit, the impulse response for 
() > 1 is of the form 


P(t) P(0) cos wole idaiadl 


P(0) = wdR/Q. 


(7) 


P(t) designates voltage in response to an impulse current in the case 


of a parallel resonant circuit, or the current in response to an impulse 


voltage in the case of a series resonant circuit. 
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L Re Sponse lo Steady State Impulse Train 


Let a long sequence of impulses of amplitude 3 and the same polarity 


impinge on a resonant circuit at uniform intervals 7. The response after 
Nv impulses is then 


nT é 


The subscript s indicates a Systematic Component. 


The above series is conveniently summed by taking the real part ol 
the series 


P(O 


4 


With ¢ 


When V 


. the steady state responses becomes 
P(O a cba il | 
9 


The interval between pulses can be written 


T = 22/w, 13 


where w is the fundamental frequency of the impulse train, or the pulse 
repetition frequency 
Let 


so that 


é 


r/( when 7/Q & 1. 





rHE BELL SYSTEM TECHNICAL JOURNAL, JULY 1957 


With these approximations 


where 


2Aw 
v = (, (18) 


« 


aw 


will be recognized as the phase shift of the resonant current at the fre 
quency w, as obtained from (3) with w = a» + Aw. 
\ further approximation that can be introduced in (12) is 


< 


~ 
— ¢ 


since fo < T and Awts and woto/2Q «K 1. 
With the above approximations (12) becomes 


P(O) Q 


2 


( ™" COS y. 
‘al part of this €XPFesslon 1s 
P(0O) Q 


2 ff 


COS (wl y) cos y, (21) 
which is the response to the steady state component of the pulse train 


5 Response to Random Component of Impulse Train 
Let a sequence of impulses of amplitude } and randomly positive and 
negative polarity impinge on the resonant circuit at intervals 7. The 
response is then, 
\ 
PD) << 


A,(0) = +? + COS wolf nT ye 


This expression for the random component differs from (9) for the sys- 


tematic component in that the impulses have random + polarity. If 
all signs are chosen the same, the values of A,(t) will be either —A,() 


or +A,(t). The resultant response of the resonant circuit, i.e. A,(t) 4 
A,(t), can thus vary between the limit 0 and 2A,(t). A,(t) represents a 
random fluctuation about A,(¢) as a mean value. In the following the 
rms value of this fluctuation is evaluated. 

In order to determine the components of A,(¢) in phase and at quadra- 
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ture with the steady state response as given by (21), it is convenient to 
write 

Ww Aw, 
COS w(t — nT) cos [w(t — nT) —p +y — Awlt — nT) 
cos la(t — nT y| COS ly Awolt — nT 
sin [w(t — nT) — | sin [y Aw(t — n7 


With ¢ NT + to, and wT’ r, (22) can be written: 


\ 
A,(t) cos (wl y) bs + COs |W — AwT(N 
N 
— sin (wl) — yp) + + sin |W, — AwT(N n)\e 
wl s 
where y; = y — Awl ¥{tI-= 0 = ¥, since wly/2Q < r/2Q < 
With equal probabilities of a plus or a minus sign in the summations, 


the rms value of the in-phase component becomes 


Awl N 


sah ta waT(N 
1 + cos Zw — AwT(N — ndie 
the quadrature component becomes 


sin” AwT(N — n)le~ | 
J 


l qe 
(1 — eos 2ly - AaTIN —n)le °* | 


%) 


These expressions can be transformed into sums of geometric series 
by writing 


COS x tie +e”), av 2ly — AwT(N 


Ny 
DD} ’ 


valuation of (25) and (26) by this method gives 


P(O) f 
) 


woT/Q ' 
»2OT YQ 
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where 
V = cos 2¥(1 — co ron’) 4 sin Qy sin 2AwT'e “°°? (29) 
D l+e- ‘ ; ~ cos 2AwT’. 50) 
With the same approximations as used previously in connection with 
(12) and with 


cos 2AwT’ =~ 1 — $ i sin 2AwT’ & 2AoT, 


With these approximations in (27 


A, & has fk (1 
= 21 

P(0) (Q 

2 2r 


which apply when y is small and (27/Q) « 


6. Response to Random Dipulse Train 


Kach dipulse is assumed to consist of two impulses of unit amplitude 
and opposite polarity, displaced by an interval 6, which in general will 
be a function of the pulse position; 1.e., 6 5(n). The response of the 
resonant circuit to a train of such dipulses, obtained by taking the dif- 


ference in response to the two impulses, is given by 


N 
A;(t) P(O) bP COS wilt — nT) 


COS wolt nT + 6(n)\e 


In determining the response, mistuning of the resonant current can be 


disregarded; 1.€., wo w. Furthermore, in the second term of (36) it is 


permissible to take 


exp [—ad(n)/2Q| = 1. 
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With the following further approximations 
COS w(t — nT’) OS woll nT + 6(n)| 


sin wolf + 5(n)/2] 2 sin [cd(n 


wd(n sin Ww t, 


expression (36) becomes: 


N 
A;(t) = P(O)ap sin wot D> 5(n)e 
N 


P(O)wo Sin woh, 2 6( nye 


where the substitution ¢ NT + ty has been made as in previous ex- 
pressions. 

The above expression shows that the resonant circuit response will be 
at quadrature with the steady state timing wave cos wo . 

In the above expressions, the dipulses are assumed to be present at 
intervals 7’, whereas in a random pulse train they will be present at 
average intervals 27’. The rms value of the quadrature component with 
randomly positive and negative dipulses at intervals 27, with an rms 
displacement 6, 1S 

. 


P(O)and | Do 


P(O) (@) 
“ T 


In (38) the funetion « oe will be recognized as the impulse response 


function of a circuit with impedance 


Z tw 


Z| tW 


tan y v/P. 3) 


It will also be recognized that (39) corresponds to the rms response of 
such a circuit, when impulses 6(n) of random amplitude with an rms value 
6 are applied to average intervals 27. Thus (39) can alternately be ob- 
tained from 
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[ a 1/2 
_ Z iw) |” de | 


1/2 
B(tan ‘w/B) | 
‘7 3) 


A (45) 


da second 


(0) aod 


P(O) WwW 6 


P(O) aod 

for 
~ 4 (‘ 
' the first resonant circuit be applied 
resonant circuit, : ‘resonant circuits, with an 
umplification 8 between successive resonant circuits. At the 
h : . : . 
resonant circuit, the rms amplitude of the response 


Let the output of 
and in turn to n successive 


amplitude : 
output of the n 
is then obtained from 


3 ‘ia 1/2 
Ain” P(O)ed | ° Bia) \" deo 
- | isT | I" ¢ | 

| +t 9 

P(O) a6 
. EE | 
PO) . 3) : 
> 08 ) 
Tv 


W here 


I. (1 


Thus (46) can be written: 


where 
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When n > 1, (51) approaches the value 


a= (=) 
ae 


The latter approximation is based on the following expression 
-3, given in Whittaker and Watson’s:** Modern Analysis”’ page 25 


n 
lim (1 ry) +a 


x 


where Tis the gamma function, '(—3 +4 


The above analysis assumes that the timing wave at each resonant 


circuit is applied directly to the next resonant circuit, except for the 
amplification between resonant circuits. This would be the case if the 
timing wave were transmitted on a separate pair, in which case A% 
would be the rms quadrature component owing to noise in the timing 
circuit. 

In regenerative repeaters, deviations in the timing wave resulting 
from the quadrature component are imparted at intervals 7 into the 
next repeater section as deviations in the spacing ot pulses. These 
timing deviations occurring at intervals 7° will have a certain random 
amplitude distribution, which can be regarded as having a certain 
frequency spectrum. When the deviations are discrete and occur at inter 
vals 7’, the spectrum will extend to a maximum frequency fia. 1/2T, 
OF @max x T wy 2. In this case the upper and lower limits of the 
integrals above would be replaced by +wy 2, except for the first repeater 
section. The recurrence relation (48) is then no longer exact, but the 
resultant modification is insignificant and can be disregarded. This will 
be seen when the value wo/2 is inserted for w in the integrand of (47 
which then becomes 1 (1 + Q°), as compared with 1 for w 0. Thus 
the contribution to the integrals for w > we/2 can for practical purposes 
be disregarded. 
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This paper considers a simple telephone exchange model which has 
infinite numbe r of t unt S and aL which the trafhe a pe nds on two paran ele 
the calling-rate and the mean holding ALLL [t is desired to estimate these 
paramete rs hy observing the model CONTINUOUSLY during a finite ante 
and noting the calling-time and hang up tame of each call. insofar as 
times fall within the interval. It is shown that the re sulling format 
for the purpose of this est mate. be reduced without loss 
These statistics are the number of calls found at the star 
number of calls arriving du ing observation, the numbe 
during observation, and the average number of ca 
inte rival of obse reatlion, The jounl distribution of these su ffi rent Statist 
determined, in principle, by deriving a generating function for 
qe nerating function the means, variances, covariances. and co 
e fhicr nts are obtained. | arious estimators tor the pa ametlers oft 
are compared, and some of thetr distributions. means. and variances 


Né nted. 


THEORETICAL PROBLEMS AND METHODS OF rRAFFI MEASUREMENT 


Four important kinds of theoretical problems arise in the measurement 
Ol telephone traffic. These are: (1) the choice of a mathematical mode 
contaming parameters characteristic of the traffic, to serve as a des rip 
tion; (2) the dey ising of efficient methods of estimating the parameters; 
(3) the determination of the anticipated accuracy OF measurements 
and (4) the assessment of actual accuracy, after measurements have 
been made. 

The present paper deals with aspects of the second and third kinds of 
problem, for the simplest and least realistic mathematical model of tele 
phone traffic. Specifically, for this model, we treat the problems of (i) 
complete extraction of the information from a given observation period, 
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without regard to costs of observation, and (ii) determination of the 
anticipated accuracy of certain methods of estimation which arise natu- 
rally from the discussion of complete extraction. 

The method by which we attack problems (i) and (ii) in this paper 
has three stages. First we choose a small number of significant properties 
of, or factors in, the physical system we are studying. Then we abstract 
these properties into a mathematical model of the physical system. Fi- 
nally, from the properties of the model, we derive results which may be 
interpreted as answers to the two problems treated. The advantage of 
this method is that we can use the precise, powerful apparatus of mathe- 
matics in studying the model; its limitation is that it yields results which 
are only as accurate as the model in describing reality. 

A method similar to the above forms the theoretical underpinning of 
telephone traffic engineering itself. To design equipment effectively, the 
traffic engineer needs a description of the traffic that is handled by central 
offices. He decides what properties of the entire system of telephone 
equipment and customers will be most useful to him in deseribing the 
traffic. He then designates certain parameters to serve as mathematically 
precise idealizations of these properties, and in terms of these parameters 


constructs a medel of the traffic, upon which he bases much of his engi- 


neering. 

In choosing’ # mathematical model for a physical system, one is con- 
fronted with tfro generally opposed desiderata: fidelity to the system 
described, and #nathematical simplicity. The model may involve impor- 
tant departuregfrom physical reality; a model that is sufficiently amena- 
ble to mathem#tical analysis often results only after one has introduced 
admittedly fae assumptions, ignored certain effects and correlations, 
and generall¥ oversimplified the system to be studied. However, the 
abstract mgdel will be an exact and simple tool for analysis. 

We canjconstruct a simple mathematical model for the operation of 
a telephgne central office by leaving out of consideration many impor- 
tant fa¢ts about such systems, and by concentrating on factors most 
relevant to operation. Since we are interested in telephone traffic and 
in the availability of plant, it seems natural to require that a realistic 
model take account of at least the following five significant factors: (1) 
the demand for telephone service; (2) the rate at which requests for 
service can be processed and connections established; (3) the lengths of 
conversations; (4) the supply of central office equipment; and (5) the 
manner in which the first four factors are interrelated. Unfortunately, 
the mathematical complexity of such a realistic model precludes easy 
investigation. Therefore, the model used in this paper is based only on 
factors (1) and (3). 
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The demand for telephone traffic is usually made precise by describing 
a stochastic process which represents the way in which requests for tele- 
phone service occur in time. A realistic description will take account of 
the facts that, the demand is not constant, but has daily extremes, and 
that in small systems, the demand may be materially lessened when 
many conversations are in progress. Since taking account of the first fact 
leads to a more complicated model in which our investigations are more 
difficult, we ignore it, with the proviso that the results we derive are 
only applicable to systems and observations for which the demand is 
nearly constant. The second kind of variation in demand becomes insig- 
nificant as the number of subscribers increases and the traffic remains 
constant. Hence, we further confine the applicability of our results to 
systems with large numbers of subscribers, and we assume that the de- 
mand does not depend on the number of conversations in existence. 

With these assumptions, a mathematically convenient description of 
the demand is specified by the condition that the time-intervals between 
requests for service have lengths which are mutually independent posi- 
tive random variables, with a negative exponential distribution. 

A telephone central office contains two kinds of equipment: control 
circuits which establish a desired connection, and talking paths over 
which a conversation takes place. The time that a request for service 
occupies a unit of equipment, be the unit a control circuit or a talking 
path, is called the holding-time of the unit. A request for service affects 
the availability of both kinds of equipment but, except for special cases, 
the holding-times of talking paths are usually much longer than the 
holding-times of control units such as markers, connectors, or registers. 
In view of this disparity, we assume that the only holding-times of con- 
sequence are the lengths of conversations; i.e., the holding-times of 
talking paths. We assume also that these lengths are mutually inde- 
pendent positive random variables, with a negative exponential distribu- 
tion. 

For the simplest mathematical model of telephone traffic, we may 
consider the arrangement of switches and transmission lines which con- 
stitutes a talking path in the physical office to be replaced by an abstract 
unit called a “‘trunk’’. A trunk is then an abstraction of the equipment 
made unavailable by one conversation, and we may measure the supply 
of talking paths in the office by the number of trunks in a model. The 
word “‘trunk”’ is also used to mean a transmission line linking two central 
offices, but as long as we have explained our use of the word there need 
be no confusion. Often the number of transmission lines leading out of 
an office is a major limitation on its capacity to carry conversations, 
and in this case the two uses of the word “‘trunk” are very similar. Un- 
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fortunately, we do not take advantage of this similarity, since we make 
the mathematically convenient but wholly unrealistic assumption that 
the number of trunks in the model is infinite. 

The model we investigate thus depends on only two of the factors 
previously listed as essential to a realistic model: namely, (1) the demand 
for service, and (3) the lengths of conversations. In view of the simplicity 
and inaccuracy of this model, the question arises whether much is gained 
from a detailed analysis. Such scrutiny may indeed reveal little that is 
of great practical value to traffic engineers. It is important methodologi- 
cally, however, to have a detailed treatment of at least one approximate 
case. We-undertake this detailed treatment largely for the insight that 
it may give into methods which could be useful in dealing with more 
complex and more accurate models. 

Once a designer has chosen a model and has specified the parameters 
he would like to have measured, it is up to the statistician to invent effi- 
cient means of measurement, by choosing, for each parameter, some 
function of possible observations to serve as an estimate of that parame- 
ter. One measure of efficiency that is of mostly theoretical interest is the 
observation time required to achieve a given degree of anticipated ac- 
curacy; the most realistic measure of efficiency is in terms of dollars and 
man-hours. It may often be more efficient, in the sense of the latter 
measure, to spread observation over enough more time to compensate 
for the inability of an intrinsically cheaper method of measurement to 
extract all of the information present in a fixed time of observation. For 
example, periodic scanning of switches in a telephone exchange is usually 
less costly than continuous observation. As a result, telephone traffic 
measurement is usually carried out by averaging sequences of instan- 
taneous periodic observations of the number of calls present, rather than 
by continuous time averaging, although it can be shown that continuous 
»bservation is more efficient at extracting information. Thus statistical 


efficiency, which may be expensive in terms of measuring equipment, 


can be exchanged for observation time, which may be less costly. This 
exchange brings about a reduction in cost without impairing accuracy. 

Our concern in this paper is with the less practical problems of com- 
plete extraction, and of the anticipated accuracy of estimation methods 
based on complete extraction. Let us consider how our mathematical 
model can shed light on these problems. A mathematical model may or 
may not be a faithful description of the behavior of real telephone sys- 
tems. Nevertheless random numbers, with or without modern computing 
machines, enable one to make experiments and observations on physical 
situations which approximate, arbitrarily closely, any mathematical 
model. Thus we can speak meaningfully of events in the model, and of 
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making measurements and observations on the model. The mathematical 
model elucidates our problems in the following ways: (1) it enables us 
to state precisely what information is provided by observation; (2) it 
enables us to explain what we mean by complete extraction of informa- 
tion; and (3) it enables us to derive results about the anticipated ac- 
curacy of measurements in the model. These results will have approxi- 
mately true analogues in physical situations to which the model is 
applicable. . 

The calls existing during the observation interval (O, 7’) fall into four 
categories: (i) those which exist at O, and terminate before 7’; (ii) those 
which fall entirely within (O, T); (iii) those-which exist at O and last 
beyond 7’; and (iv) those which begin within (O, T) and last beyond T. 
For calls of category (i), we assume that we observe the hang-up time 
of each call; for category (ii), we observe the matching calling-time and 
hang-up time of each conversation; for category (ili), we observe simply 
the number of such calls; and for category (iv), we observe the calling- 
times. Table I summarizes the kinds of calls and the information ob- 
served about each. 

What we mean by the complete extraction of information is made 
precise by the statistical concept of sufficiency. By a statistic we mean 
any function of the observations, and by an estimator we mean 
a statistic which has been chosen to serve as an estimate of a particular 
parameter. Roughly and generally, a set S of statistics is sufficient for a 
set P of parameters when S contains all the information in the original 
data that was relevant to parameters in P. If S is sufficient for P, there 
is a set E of estimators for parameters in P, such that the estimators in 
E depend only on statistics from S, and such that an estimator from E 
does at least as well as any other estimator we might choose for the same 
parameter. Thus we incur no loss in reducing the original data (of speci- 
fied form) to the set S of statistics. It remains to state what it means for 
S to contain all the relevant information. We do this in terms of our 
model. 

The mathematical model we are adopting contains two oft ribution 


TABLE [| INFORMATION OBSERVED f 


Types of Calls Start in (O, T Start Kefore O 


End in (O, 7)| (ii), matching calling-times and | (i), hang-up times known, num- 
hang-up times known, num- ber of calls known 
ber of calls known 

End after 7 (iv), calling-times known, num- | (iii), number of calls known 
ber of calls known 
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functions, that of the intervals between demands for service, and that 
of the lengths of conversations. We have supposed that these distribu- 
tions are both of negative exponential type, each depending on a single 
parameter. Thus we know the functional form of each distribution, and 
each such form has one unknown constant in it. Since the mathematical 
structure of the model is fully specified except for the values of the two 
unknown constants, we can assign a likelihood or a probability density 
to any sequence > of events in the model during the interval (O, T). 
This likelihood will depend on the parameters, on Y, and on the number 
of calls in existence at the start O of the interval. If the likelihood L(2) 
can be factored into the form L = F-H, where F depends on the param- 
eters and on statistics from the set S only, and H is independent of the 
parameters, then the set S of statistics may be said to summarize all the 
information (in a sequence 2) relevant to the parameters. If L can be 
so factored, then S is sufficient for the estimation of the parameters. 

The mathematical model to be used in this paper is described and 
discussed in Sections II and III, respectively. Section IV contains a 
summary of notations and abbreviations which have been used to sim- 
plify formulas. 

In Appendix A we show that the original data we have allowed our- 
selves can be replaced by four statistics, which are sufficient for estima- 
tion. In Appendix B and Sections V-VIII we discuss various estimators 
(for parameters of the model) based on these four statistics. To determine 
the anticipated accuracy of these methods of measurement, we consider 
the statistics themselves as random variables whose distributions are 
to be deduced from the structure of the model. 

A primary task is the determination of the joint distribution of the 
sufficient statistics. In view of the sufficiency, this joint distribution tells 
us, in principle, just what it is possible to learn from a sample of length 
T in this simple model. By analyzing this distribution we can derive 
results about the anticipated accuracy of measurements in the model. 

The joint distribution of the sufficient statistics is obtainable in prin- 
ciple from a generating function computed in Appendix C, using methods 
exemplified in Section X. This generating function is the basic result of 
this paper. The implications of this result are summarized in Section 
IX, which quotes the generating function itself, and presents some 
statistical properties of the sufficient statistics in the form of four tables: 
(i) a table of generating functions obtainable from the basic one; (ii) a 
table of mean values; (iii) a table of variances and covariances; and 
(iv), a table of squared correlation coefficients. (The coefficients are 
all non-negative.) 
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If DESCRIPTION OF THE MATHEMATICAL MODEL 


Throughout the rest of the paper we follow a simplified form of the 
notational conventions of J. Riordan’s paper" wherever possible. A sum- 
mary of notations is given in Section IV. The model we study has the 
following properties: 

(i) Demands for service arise individually and collectively at random 
at the rate of a calls per second. Thus the chance of one or more demands 
in a small time-interval Af is 


aAt + o(At), 


where o(At) denotes a quantity of order smaller than At. The chance of 
more than one demand in At is of order smaller than At. It can be shown 
(Feller,’ p. 364 et seq.) that this description of the demand is equivalent 
to saying that the intervals between successive demands for service are 
all independent, with the negative exponential distribution 


l-—e”. 
This again is equivalent to saying that the call arrivals form a Poisson 


process; i.e., that for any time interval, ¢, the probability that exactly n 
demands are registered in ¢ is 


e “(at)” 


! 


M. 


Thus the number of demands in ¢ has a Poisson distribution with mean 
at, 

(ii) The holding-times of distinct conversations are independent vari- 
ates having the negative exponential distribution 
— 
where y is the reciprocal of the mean holding-time h. This description of 
the holding-time distribution is the same as saying that the probability 
that a conversation, which is in progress, ends during a small time- 
interval At is 

yAt + o(Atb), 


without regard to the length of time that the conversation has lasted 
Feller, p. 375). 

(iii) The model contains an infinite number of trunks. Thus, at no time 
will there be insufficient central office equipment to handle a demand 
for service, and no provision need be made for dealing with demands that 
cannot be satisfied. 
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The original work on this particular model for telephone traffic is in 
Palm, and Palm’s results have been reported by Feller® and Jensen.‘ 
The results have been extended heuristically to arbitrary absolutely con- 
tinuous holding-time distributions by Riordan," following some ideas of 
Newland’ suggested by 8S. O. Rice. 

Let P;,;(t) be the probability that there are j trunks busy at ¢ if there 
were 7 busy at 0. And let P(t, x) be the generating function of these 
probabilities, defined by 


oo 


Pt,2) = > a’P,,(0. 


7=0 
Then Palm” has shown (pp. 56 et seq.) that 
Pt, z) = [1+ @-Dey exp {(x — l)ah (1 — & )}. 


a . - . , ll . ‘ 
his is formula (12) of Riordan” with his g replaced by e ”*. It can be 
verified that the random variable N(t) is Markovian; the limit of P,(t, x) 


ast—>» © is 
exp {(x — 1) ah}, 


so that the equilibrium distribution of the number of trunks in use is a 
Poisson distribution with mean b = ah. The shifted random variable 
(N(t) — 6] then has mean zero, and covariance function be™”’. 

For additional work on this model the reader is referred to F. W. 
Rabe,” and to H. Stormer.” 


III DISCUSSION OF THE MODEL 


Let us envisage the operation of the model we have described by con- 
sidering the random variable N(¢) equal to the number of trunks busy 
at time ¢. As a random function of time, V(¢) jumps up one unit step each 
time a demand for service occurs, and it jumps down one unit step each 
time a conversation ends. If N(¢) reaches zero, it stays there until there 
is another demand for service. If N(t) = n, the probability that a con- 
versation ends in the next small time-interval At is 


nyAt + o(At), 


because the n conversations are mutually independent. A graph of a 
sample of N(é) is shown in Fig. 1. 

The model we described departs from reality in several important 
ways, which it is well to discuss. First, the assumption that the number 
of trunks is infinite is not realistic, and is justified only by the mathe- 


matical complication which results when we assume the number of trunks 
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TIME, t —> 
Fig. 1 — A graph of N(t). 


to be finite. It can also be argued that unlimited office capacity is ap- 
proached by offices with adequate facilities and low calling rates, and 
therefore, in some practical cases at least, the model is not flagrantly 
inaccurate. 

Second, the choice of a constant calling rate for the model ignores the 
fact that in most offices the calling rate is periodic. Thus, the applica- 
bility of our results to offices whose calling rates undergo drastic changes 
in time is restricted to intervals during which the normally variable 


calling rate is nearly constant. Finally, although the assumption of a 
negative exponential distribution of holding-time affords the model great 
mathematical convenience, it is doubtful whether in a realistic model 
the most likely holding-time would have length zero, as it does in the 
present one. 


IV SUMMARY OF NOTATIONS 
a = Poisson calling rate 
h = mean holding-time 
y = h~ = hang-up rate per talking subscriber 
b = ah = average number of busy trunks 

N(t) = number of trunks in use at ¢ 

(O, T) interval of observation 
n N(O) = number of trunks in use at the start of observation 
A = number of calls arriving in (O, T) 
H = number of hang-ups in (0, 7) 


K=A+H 
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$ 
Z = [ N(t) dt 


M Z/T = average of N(t) over (O, T) 


tpn} the (discrete) probability distribution of n, the number of 


trunks found busy at the start of observation 


An estimator for a parameter is denoted by adding a cap (*) and a 
subscript. The subscripts differentiate among various estimators for the 
same parameter. We use 4, = A/T, Ye = H/Z, 4, = K/2T, 
v1 = K/2Z, and x2 = A/Z. 

Also, it is convenient to use the following abbreviations: r for y7', and 
C for (1 — e ")/r, where r is the dimensionless ratio of observation-time 


to mean holding-time. The symbol F is used throughout to mean mathe- 
matical expectation. 


V THE AVERAGE TRAFFIC 


We have adopted a model which depends on two parameters, the 
calling rate a, and the mean holding-time h, or its reciprocal y. Before 
searching for a set of statistics that is sufficient for the estimation of 
these parameters, let us consider the product ah = b. This product is 
important because, as we saw in Section II, the equilibrium distribution 
of the number of trunks in use depends only on b, and not on a and h 
individually. Indeed, the equilibrium probability that n trunks are busy 


1S 


and the average number of busy trunks in equilibrium is just b. 
The average number of trunks busy during a time interval T' is 


& + 
M = > N(t) dt; 
T' Jo 


i.e., the integral of the random function N(¢) over the interval 7’, divided 
by 7. This suggests that for large time intervals 7’, M will come close 
to the value of b, and can be used as an estimator of b. Since M is a ran- 
dom variable, the question arises, what are the statistical properties of 
M? This question has been considered in the literature, the principal 
references being to F. W. Rabe” and to J. Riordan.” Riordan’s paper is 
a determination of the first four semi-invariants of the distribution of M 
during a period of statistical equilibrium, but without restriction on the 
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assumed frequency distribution of holding-time. It follows from Rior- 
dan’s results that M converges to b in the mean, which is to say that 
lim EF {| M — b/’} =0. 
T+2 
It also follows that M is an unbiased estimator of b; i.e., that BE} MW} 
and that M is a consistent estimator of b, which means that 
lim pr{| M—b| > e} =0 


T+2 


for each ¢ > 0. 


VI MAXIMUM CONDITIONAL LIKELIHOOD ESTIMATORS 


As shown in Appendix A, the likelihood L, of an observed sequence, 
conditional on N(Q), is defined by 


nL. = Alna+HIny — 7Z — al. 


According to the method of maximum likelihood, we should select, as 
estimators of a and y respectively, quantities @. and 7. which maximize 
the likelihood L.. Now a maximum of L, is also one of In L. , and vice 
versa. Therefore a, and y, are determined as roots of the following two 
equations, called the likelihood equations: 


0 9 
—iInL, = 0; fi In L. = 0. 
da oY 


The solutions to the likelihood equations are 


H 

7: 
These are the maximum conditional likelihood estimators of a and y. 
The estimator 4, is the number of requests for service in 7 divided by 7; 

this is intuitively satisfactory, since @, estimates a calling rate. 

Since maximum likelihood estimators of functions of parameters are 
generally the same functions of maximum likelihood estimators of the 
parameters, we see that AZ/HT is a maximum likelihood estimator of b. 


VII PRACTICAL ESTIMATORS SUGGESTED BY MAXIMIZING THE LIKELIHOOD 
L, DEFINED IN APPENDIX A 


We obtain as likelihood equations 


q a 
Swe & osha 6 
da Oy 
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These may be written as 


n+ A 


ge 
- 


z+" 


The first of these shows the estimated calling rate as a pooled combination 
of the conditional estimate A/T’, considered in the last section, and an 
estimate n/h based on the initial state. This latter estimate has the form 


calls in progress 


mean holding time’ 


and so is intuitively reasonable, since b/h = a. The second equation 
exhibits our estimate of -y as a pooled combination of the conditional 
estimate H/Z and the ratio a/n. This ratio is acceptable as an estimate 
of y, since a/b = y andb = En} is the average value of n. 

If we substitute, in the right-hand sides of these equations, the condi- 
tional estimators A/7', H/Z, and Z/H for a, y, and h, respectively, we 
obtain simple, intuitive estimators which include the influence of the 
initial state n, and show how it decreases with increasing 7’. Thus 


n + A 
Z ot 
mt 


estimates a, 
AZ 
pou: 
p nZ 
ae 


H + 


estimates y. 


Vill OTHER ESTIMATORS 


Additional estimators may be arrived at by intuitive considerations, 
or by modifying certain maximum likelihood estimators. Some estimators 
so obtained are important because they use more of the information 
available in an observation than do the conditional estimators 4, and 
7¥-, Without being so complicated functionally that we cannot easily 
study their statistical properties. 
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It seems reasonable, and can be shown rigorously (Appendix C), that 
for an interval (O, 7) of statistical equilibrium, the distribution of A 
and that of H are the same. Thus we can argue that, for long time inter- 
vals, A and H will not be very different. This suggests using 

Pe ty eee 


ay = 


oT OT 
as an estimator of a. This estimator does not involve y, and it uses not 
only information given by A, but also information supplied by arrivals 
occurring possibly before the start of observation. 

Similarly, since b = a/y, and M isa consistent and unbiased estimator 
of b, we may use 

K l 
N=s5=7 
9 
22 
to estimate y, and its reciprocal to estimate h. Finally, since for long (O,7') 
we have A ~ H, we may try 
l 
Z he 


as an estimator of 7, and its reciprocal as an estimator of h. 


IX THE JOINT DISTRIBUTION OF THE SUFFICIENT STATISTICS 
The basic result of this paper is a formula for the generating function 
E627 ws ute | (9.1) 
for the joint distribution of the random variables n, N(7), A, H, and Z. 


This formula is derived in Appendix C, by methods illustrated in Section 
X. For an initial n distribution {p,}, the generating function is 


>, Pat” |= + - we +t mT 
n>0 ¢ + 7 


exp (aude + ve — volt - 7) | 
pi G+)? cF3 


aywuTl = 

It is proved in Appendix A that the set of statistics {n, A, H, Z} is 
sufficient for estimation on the basis of the information assumed, which 
was described in Section I. Thus the generating function (9.2) specifies, 
at least in principle, what can be discovered about the process from an 
observation interval (O, 7), for which N(O) has the distribution {p,}. 
All the results summarized in this section are consequences of (9.2). 
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TABLE II 


In E{X} 


. eT a aed 
a a _ = 


P+) za 


( c2 r2(1 — e+ 
Sh ie OR ; mPa 
| pt ee +r)? | 


2aTC(y — 1) + aT(l — C)(y? — 1) 


2aTC(et/?T — 1) + aT(l — C)(e*/7? — 1) 

. 2 pay? 

b (: —— ) fection) hae (1 - 2 
i +.e o +7 


By substitution, and by either letting the appropriate power series 
variables — 1, or letting ¢ — 0, or both, we can obtain from (9.2) the 
generating function of any combination of linear functions of the basic 
random variables n, N(7T), A, H, and Z. Some of the generating func- 
tions thereby obtained are listed in Table II, in which the entries all 
refer to an interval (O, 7’) of equilibrium. 

Since, for equilibrium (O, 7), the generating functions are all exponen- 
tials, it has been convenient to make Table II a table of logarithms of 
expectations, with random variables X on the left, and functions In 
E{X} on the right. C as a function of r is plotted in Fig. 2. 

Entry 1 of Table II is actually the cumulant generating function of 
Z for equilibrium (O, 7); similarly, Entry 2 is that of M/, and depends 
only on the average traffic b and the ratio r. The form of the general 
cumulant of M is 


~T 


L) | (T — x)x” e ™ dx. 


nin — 
) 


k, = : 
7 n 


This result coincides with a special case (exponential holding-time) of 
a conjecture of Riordan.’ This conjecture was first established (for a 
general holding-time distribution) in unpublished work of 8. P. Lloyd. 
The cumulant generating function permits investigation of asymptotic 
properties. We prove in Section X that the standardized variable 
v = (yT/2b)'" (M — b) 
= (r/2b)"” (M — b) 


is asymptotically normally distributed with mean 0 and variance 1. 
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1.0 


Fig. 2— C as a function of r. 


From Entry 3 of Table II it can be seen that A is distributed as 
2u + v, where u and v follow independent Poisson distributions with the 
respective parameters a7'(1 — C) and 2aTC. The probability that K = 
n for an interval of equilibrium is 


(n — 


(2aTC)"” (aT — aTC)’ 


r, = exp {aT(C — 1)} p> 


2))! J! 
where the sum is over j’s for which 0 S 2) S n. 

The estimator 4, for a is equal to K/27T, and has mean and variance 
given by 


FE} 4,} 
jar {ay} 


The distribution of 4, is given by 
priaq, S x} = » Sa %. 


the summation being over n S 27x. 
From (9.2) one can obtain, by substitution of the stationary n distribu- 
tion for {p,}, and subsequent differentiation, the means, variances, 


covariances, and correlation coefficients of the sufficient statistics, for 
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TABLE III 


H 


aT 
aT(C + b) 


aT 
baT 
aT(i + aT) 


b(1 + 5) 


aT(1 + aT) 


A H 


aT(i—C+aT 


-E{X,Y} 
K 


2aT 

aT(C + 2b) 
) aT(2 — C + aT) 
aT(2 — C + aT) 
2aT(2 — C + 2aT) 


TaBLE IV — cov{X,Y} 


Z 


bT 
bT(C + b) 
obT(l — C + aT) 
bT(1 — C + aT) 
2bT(1 — C + aT) 
bTA(2 — 2C + aT) 





0 aTC 
aT(1 — C) 
aT 


aTC 


) 
2aT (2 — C) 


bTC 

bT(1 — C) 
obT(1 — C) 
2b7 (1 — C) 


2bTh(1 — C) 


equilibrium intervals (O, 7). It has been convenient to display these in 
three triangular arrays, the first consisting of expectations of products, 
the second comprising the variances and covariances, and the third 
exhibiting, for simplicity, the squared correlation coefficients, since the 
correlation coefficients are never negative for these random variables 

In Table III, the entry with coordinates (X, Y) is {XY} for equilib- 
rium (O, 7). All three tables are expressed in terms of a, b, T, h, r, and 
C’, the last of which is plotted in Fig. 2. 

The variances and covariances of the sufficient statistics are listed in 
Table IV; the entries are of the form: 


cov {X, Y} = E{XY} — E{XJE{Y}. 


Table V, finally, lists the squared correlation coefficients; i.e., the 


quantities 


cov {X, Y} 


p (X, } —— 
var {|X} var {Y} 
For any time interval (O, 7), A has a Poisson distribution with param- 
eter aT, so that 74, does also. Therefore the distribution of 4, is given 


by 


_ pany 


n! 


< 2} 


pr {d. 


where the summation is over n S 2x7’. Evidently 
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var {4,} 


so that @, is an unbiased and consistent estimator of a. We now compare 
the variances of estimators 4. and d,. From Table IV we have 


var {d:} = . l1— | < . = var {4} 
ss t — 7 — c i Cie 
7 ; 1 


so that 4, is a better estimator of a for any 7 > 0, in the sense that its 
variance is less. 


X THE DISTRIBUTIONS of Z AND | 


Since we have defined 


aT 


| N(t) dt, 


we can regard Z as the result of growth whose rate is given by the ran- 
dom step-function N(t); when V(t) = n, Z is growing at rate n. An idea 
similar to this is used by Kosten, Manning, and Garwood’, and by Kos- 
ten alone.” Now the Z(7’) process by itself is not Markovian, but it can 
be seen that the two-dimensional variable |N(¢), Z(¢)} itself is Marko- 
vian. Let F,(z, t) be the probability that N(#) = n and Z(t) S z. Since 
the two-dimensional process is Markovian, we can derive infinitesimal 
relations for F(z, t) by considering the possible changes in the system 
during a small interval of time At. 


TABLE V — p(X, Y) 
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If N(t) = n, then the probability is [1 — ynAt — adAt — o(At)} that 
there is neither a request for service nor a hang-up during At following 
t, and that Z(t + At) = Z(t) + nAt. Therefore the conditional proba- 
bility that N(t + At) = n and Z(t + At) S z, given that no changes 
occurred in Af, is 


F(z — nAt, t). 


For N(t) = (n + 1), the probability is y(n + 1)At + o(At) that one 
conversation will end during At following ¢. The increment to Z(t) 
during At will depend on the length x of the interval from ¢ to the point 
within At at which the conversation ended. The increment has magni- 
tude (n + 1)z + n(At — x) = x + nAt, as can be verified from Fig. 3, 
in which the shaded area is the increment. Since z is distributed uniformly 
between 0 and At, the increment x + nAt is distributed uniformly be- 
tween nAt and (n + 1)At. Therefore the conditional probability that 
N(t + At) = n and Z(t + At) S z, given that one conversation ended 
in At, is 

«(n+1) At 


x es Frayilz — u, t) du. 

By a similar argument it can be shown that the probability that one 
request for service arrives in At is aAt + o(At), and that the conditional 
probability that N(t + At) = n and Z(t + At) S z, given that one 
request arrived during At, is 

] nAt 


bs F,-1(2 — u, t) du. 
At (n—1)At if ; 


Define F,,(z, t) to be identically 0 for negative n. Adding up the probabil- 





| 
| 








ig. 3 — Increment to Z in At. 
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ities of mutually exclusive events, we obtain the following infinitesimal 
relations for F,,(z, t): 


(n+1) At 
F(z, § + At) = y(n + 1) [ Frsilz — u, t) du 
“nt 


nat 
+a [ F,-(z — u, t) du + F,(z — nAt, t) 
¥(n—1)At 


-{1 — At(yn + a)| + ofAd), for any n. 


Expanding the penultimate term of the right side in powers of nAt, 
and the left side in powers of At, we divide by Af, and take the limit as 
At approaches 0. Now 


(n+1) At 

lim | Frailz — u,t) du = Fyxi(z, 2). 
ato At Jnat 

Thus, omitting functional dependence on z and ¢ for convenience, we 

reach the following partial differential equations for F(z, ¢): 

0 0 


== n= 1 F, ) fa tig-20 4 
rtd y(n + 1)F agi + 4 1 n a, f (10.1) 


—|yn+al)F,, forany n. 
Since Z(0) = 0, we impose the following boundary conditions: 
F,(0,t) =O for n>Oandt> 0, 
F,(z, 0) — ee eG (10.2) 
F.(z,0)=0 for z<0O, 


where the sequence {p,} forms an arbitrary N(0) distribution that is 
zero for negative n. 

To transform the equations, we introduce the Laplace-Stieltjes in- 
tegrals 


20 


galt, t) = | e'dF,(z,), t20, Re(s)>0, 


in which the Stieltjes integration is understood always to be on the 
variable z. We note that 
2 1 


oP (et) de = — elt, 0), 
“0 : 


and that 


2 


— 
on(t, t) = F,(0, t) + ef = F,, (z, t) dz. 
d oz 
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Applying now the Laplace-Stieltjes transformation to (10.1), we obtain 


On 


ie y(n + Wgnsi + Agni — ni&en + ntF, (0, t) 
Cc 


(10.3) 


— [yn + ale, , 
in which we have left out functional dependence on ¢ and ¢ where it is 
unnecessary. By the boundary conditions (10.2), n¢F,(0, t) = 0 for 
n = Oandt > 0; in (10.3) we may therefore omit this term in the region 
t > 0. Let ¢ be defined by 


ro 


g(x, f,t) = > 2"y, (f, 0). 


n=) 
The series is absolutely convergent for | x | < 1, since 
| gn(t, t)| S 1, for all n. 


The following partial differential equation for ¢ is obtained from 
(10.3): 


*) . 
oF kts + old's Seo oe (10.4) 
ot Ox 


If we integrate out the information about Z by letting approach 0 in this 
equation, we obtain the equation derived by Palm (loc. cit.) for the gene- 
rating function of N(t). Therefore our equation has a solution of the 
same form as Palm’s. For the boundary conditions (10.2), this solution is 


att 
i+) 


> ‘ & + y(x — Ie FP + af 


n=) ¢ aa Y 


(10.5) 


Actually g contains more information than we want since it yields the 
joint distribution of N and Z. We may integrate out the former variable 
by letting x approach 1 in 10.5. Then, 


= ({{(+yT ry 
Etexp (—tZ))-= exp jas es ) al 
(¢ + y)? t+ y 


n=0 


is the Laplace transform of the distribution of Z for an arbitrary N(O) 
distribution {p,}. This result is not restricted to an interval (O, 7) 
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of statistical equilibrium; however, if the sequence {p,} does form the 
stationary distribution discussed in Section II, then 


b a "Dn — exp {b(x mune 1)}, (10.6) 


and 


{ ae, ak ay Al 
y = exp {bgt - Y) i. I. (10.7) 


C++ Depa 2 
is the Laplace transform of the distribution of Z for an interval (O, T) 
of statistical equilibrium. 
The Laplace transform is a moment generating function expressible as 
y= Eid 


rom | n! 


‘ th ‘ = — eT . _ 
where m, is the n” ordinary moment of Z. Differentiation of 10.7 then 
gives a recurrence relation for the moments upon equating powers of 


(—¢). Thus, 
(¢ + y)° (-%) 


of 
= p- {2ag(1 — 8°") + 6 + hal + bT¢* 
and 
y'Mn41 — 3y'nm, + 3yn(n — 1)ma-1 — n(n — 1)(n — 2)my_2 
a ay Tm, — (2a + ayT)nm,1 + 2ane""(m + T)"' +n (10.8) 
-(n —1)aTe "7 (m + TT)" — n(n — 1)(n — 2)bTe "(m+ T)"™, 


n 


is the usual symbolic abbreviation of 


+ (") TP" tte4 
i=0 \J 


From the recurrence (10.8) it is easily verified that 


m, = bT, 


where (m + 7’) 


2 | ‘ry 
m2, = (bT)° + = 1 — Cl, 
+ 
from which it follows that the variance of Z is 


9 ‘yy’ 
var {Z) = et — Ch. 
Y 
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Since y is the Laplace-Stieltjes transform of the distribution of Z 
over an interval of equilibrium, In y is the cumulant generating function, 
and has the following simple form: 


2 " eeagried —~ ) is ur | 
c+ 


“4 ; (10.9) 
2 e(1 ae | 


(+ 7)? 


M is a linear function of Z, so we may obtain the cumulant generating 
. . . . 1 = rm: . . 
function of M in accordance with Cramér (p. 187). This function is 


ue ? e Pl ane eve) 
b|-++ fF, - ae | se 


and depends only on 6 and r. 
The mean and variance of M for an interval of equilibrium are respec- 
tively given by 


2b l—e 
var {M} = — [1 — C], with C = 
r r 
. . : ° lu . . . ° 
results which were first proved in Riordan.” A normal distribution having 
the mean and variance of M has the cumulant generating function 


b | s rca (10.11) 


which is to be compared to (10.10). Since var {M} goes to 0 as T ap- 
proaches ©, we may expect that a suitably normalized version of Z will 
be asymptotically normally distributed as T approaches «. The cumu- 
lant generating function of the normalized variable (2bhT)*(Z — bT) 


9 2 95\-12 : ) : 
2 + | c (**) re : be 


=e ee ee 
2 \12 2b\72 

Estee 9 » Fi 

‘(3) + 2 ‘(?) +r 


which approaches ¢’/2 as T — . It follows that the normalized variable 
is asymptotically normal with mean 0 and variance 1, and that 
(r/2b)*(M — b) is also asymptotically normal (0, 1). 


iS 
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APPENDIX A 


PROOF THAT {n, A, H, Z} 1S SUFFICIENT. 


We observe the system during the interval (O, 7), and gather the in- 
formation specified in Section I, and summarized in Table I. From this 
information we can extract four sets of numbers, described as follows: 

S, the set of complete observed inter-arrival times, not counting 

the interval from the last arrival until T 

S, the set of complete observed holding times 

S, the set of hang-up times for calls of category (i) 

S, the set of calling-times for calls of category (iv) 

In addition, our data enable us to determine the following numbers: 
n the number N(Q) of calls found at the start of observation 
k thenumber of calls of category (iii); i.e., of calls which last through- 
out the interval (0, 7) 
the length of the time-interval between the last observed arrival 
and T 

In view of the negative exponential distributions which have been 
assumed for the inter-arrival times and the holding-times, and in view 
of the assumptions of independence, we can write the likelihood of an 
observed sequence of events as 


b= hp TT oe TT eT oT, 
zeSp 


utS, weS; yeS,4 


so that 


InL = —ykT — ax + Inp, + Alna — » au 


ueSg 


+Hiny — Dy - Dw LT — y) 
zeS, weS; yeS, 
It is easily seen that the summations and the two initial terms can be 
combined into a single term, so that we obtain 


mnL=hnp,+Alna+AHIiny — yZ — aT. 


This shows that L depends only on the statistics n, A, H, and Z; it 
follows that the information we have assumed can be replaced by the 
set of statistics {n, A, H, Z}, and that these are sufficient for estimation 
based on that information. 

The likelihood is sometimes defined without reference to the initial 
state, by leaving the factor p, out of the expression for L. Strictly speak- 
ing, this omission defines the conditional likelihood for the observed 
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sequence, conditional on starting at n. We use the notation: 


A definition of likelihood as L, has been used by Moran.’ Clearly 


nL=Ala+AHIny — yZ — al. 


APPENDIX B 


UNCONDITIONAL MAXIMUM LIKELIHOOD ESTIMATES 


+ 


The definition of likelihood as L leads to complicated results which 
are of theoretical rather than practical interest. For this reason these 
results have been relegated to an appendix. 

The results of setting 0/dy In L and 0/da In L equal to zero lead, re- 
spectively, to the likelihood equations 


a — y(n — H) —-yZ = 0, 
yn — a+ yA —ayT = 0. 


Considered as a system of equations for y and a, this pair has the non- 


negative roots 


_H-n—-M + {(H ~n—~ M)* + 4MK)" 


2Z 





Kou 
These are the unconditional maximum likelihood estimators for y and a. 
Although @, depended only on A and 7’, and ¥, only on H and Z, the 
unconditional estimators depend on all of n, A, H, Z, and T. We may 
obtain a maximum unconditional likelihood estimator for b as well, 
either by considering L to be a function of 6 and y, or from general 
properties of maximum likelihood estimators. Since b = a/y, we expect 
that 6 = 4/4, as can be verified by an argument similar to that used 
above for 4 and 4. 

The estimators 4, b, and 7 obtained in this Appendix may turn out to 
be useful in practice, but their complicated dependence on the sufficient 
statistics n, A, H, and Z makes a study of their statistical properties 
difficult. As a first step along such a study, we have derived the gen- 
erating function of the joint distribution of the sufficient statistics in 
Appendix C. The greater simplicity of the conditional estimators of 
Section VI makes it possible to study their statistical properties. This 
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fact gives them a practical ascendancy over the unconditional estimators, 
even though the latter may be more efficient statistically by dint of 
using all the information available in an observation. 

APPENDIX C 
THE JOINT DISTRIBUTION OF N(t), n, A, H, AND Z 


By methods already used in Section X one can obtain a _ gen- 
erating function for the joint distribution of all the random variables n, 
N(0), A, H, and Z. Let 


7 Nit) A H -—{Z 
@= Kix ‘wu e*”}. 


Then # satisfies the differential equation 


ID ID . 
— + [gx + yx — yu] — = a(we — 1)@ 
ot Ox 


whose solution has the form 
® = R{ltx + yx — yule 8} 
awltx 2 — yulli — & St"! aywut 
exp) [fe + ye — yulll 4 aywut _ a), 


G+ OT Fes 


where the function R is determined by the initial distribution {p,} 


through the relation 
pt ml 
R{é} = i 
8) 2, 2 k 79 


From these results it follows that the generating function 
, N(T) A Zz 
E{z"x me 


A Hz — 
4 a wu eé 


is given by 


a (S + yx — yule St?* + =) 


he ae 
— _ ester — 
exp (suet x yu{l — e 4 aywuT 


Pné 


(+ 7) 2 
If {pn} forms the stationary distribution, this reduces to 
as (2 + yx ci Al + yus 1) 
c+ 47 
aw(tx + yx — yu)[l — e St”? rywuT 
4 alge + 7 ; 0 é -o 
(¢ + y)? f+ y 
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If, in this last expression, we let a approach 1, 2 approat h 1. and 


nroach 1. we obtain 


is the weneratine function FF} ‘eo **) fear an interval of equilibrium 


Alternately. if instead we let approach 1, 2 approach 1. and w ap 
we obtain (C) with ubstituted for w: this implies the not 


proach | 
for an interval of equilibrium, the two-dimen 


urprising result that 
ional variables {|A. Z!} and |H. Z} have the same distribution. From 


and (C) it follows that for equilibrium (0. 7), G) A and A both 


thi 
have & Poisson distribution with mean a7’. and (ii) the estimators A. and 


h, have the same distribution 
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Fluctuations of Telephone Traffie 


By V. E. BENES 


(Manuscript received November 9, 1956 


The number of calls in progress in a simple telephone exchange model 
characterized by unlimited call capacity, a general probability density of 
holding-time, and randomly arriving calls is defined as N(t). A formula, 
due to Riordan, for the generating function of the transition probabilities 
of N(t) ts proved. From this generating function, expressions for the co- 
variance function of N(t) and for the spectral density of N(t) are determined. 
It is noted that the distributions of N(t) are completely specified by the co- 
variance function. 


I INTRODUCTION 


The aim of this paper is to study the average fluctuations of telephone 
traffic in an exchange, by means of a simple mathematical model to 
which we apply concepts used in the theory of stochastic processes and 
in the analysis of noise. 

The mathematical model we use is based on the following assumptions: 
(1) requests for telephone service arise individually and collectively at 
random at an average rate of a per second; (2) the holding-times of 
calls are mutually independent random variables having the common 
probability density function h(u); and (3) the capacity of the exchange 
is effectively unlimited, and no call is blocked or delayed by lack of 
equipment. This telephone exchange model has been described by J. 
Riordan.°® 

As a measure of traffic, it is natural to use the number of calls in prog- 
ress in the exchange. We are thus led to consider a random step-function 
of time N(é), defined as the number of calls in progress at time ¢t. N (i) 
fluctuates about an average in a manner depending on the calling-rate, 
a, and the holding-time density, h(u). 


II PROOF OF RIORDAN’S FORMULA FOR TRANSITION PROBABILITIES 


Let P»,.(t) be the probability that n calls are in progress at ¢ if m 
calls were in progress at 0. Define the generating function of these prob- 
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\ definition of likelihood as L. has been used by Moran.’ Clearly 
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\pPENDIX B 


NCONDITIONAL MAXIMUM LIKELIHOOD ESTIMATES 


lefinition of likelihood as L leads to complicated results which 
are of theoretical rather than practical interest. For this reason these 
results have been relegated to an appendix 

The results of setting 0/04 In 1 and @/da \|n L equal to zero lead, re 


pectively to the likelihood equations 


a y(? H yz 0) 
n a 5 | a’ 7 0 
(Considered as a system of equations for y and a, this pair has the non- 


negative roots 


H n M+ '(H n Vv)? +4MK\'" 


These are the unconditional maximum likelihood estimators for y and a. 
Although 4. depended only on A and T, and ¥, only on H and Z, the 
unconditional estimators depend on all of n, A, H, Z, and T. We may 
obtain a maximum unconditional likelihood estimator for b as well, 
either by considering L to be a function of b and y, or from general 
properties of maximum likelihood estimators. Since b a/y, we expect 
that 6 a/y, as can be verified by an argument similar to that used 
above for @ and 4. 

The estimators 4, 6, and y¥ obtained in this Appendix may turn out to 
be useful in practice, but their complicated dependence on the sufficient 
statistics n, 4, H, and Z makes a study of their statistical properties 
difficult. As a first step along such a study, we have derived the gen- 
erating function of the joint distribution of the sufficient statistics in 
Appendix C. The greater simplicity of the conditional estimators of 
Section VI makes it possible to study their statistical properties. This 
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If, in this last expression, we let « approach |, z approach 1, and u ap- 
proach 1, we obtain 


yw be( ae 
exp] (1 — - ts ( = a or) (C) 
ea: ; +4 


$Z) 


as the generating function 1} w*¢ } for an interval of equilibrium. 


Alternately, if instead we let x approach 1, z approach 1, and w ap- 
proach 1, we obtain (C) with u substituted for w; this implies the not- 
surprising result that for an interval of equilibrium, the two-dimen- 
sional variables {|A, Z} and |H, Z} have the same distribution. From 
this and (C) it follows that for equilibrium (O, T), (i) A and H both 
have a Poisson distribution with mean a7’, and (ii) the estimators A, and 
h, have the same distribution. 
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The number of calls in progress in a simple telephone exchange model 
characterized by unlimited call capacity, a general probability density of 
holding-time, and randomly arriving calls is defined as N(t). A formula, 
due to Riordan, for the generating function of the transition probabilities 
of N() ts proved. From this generating function, expressions for the co- 
variance function of N(t) and for the spectral density of N(t) are determined. 
It is noted that the distributions of N(t) are completely specified by the co- 


variance function. 


I INTRODUCTION 


The aim of this paper is to study the average fluctuations of telephone 
traffic in an exchange, by means of a simple mathematical model to 
which we apply concepts used in the theory of stochastic processes and 
in the analysis of noise. 

The mathematical model we use is based on the following assumptions: 
(1) requests for telephone service arise individually and collectively at 
random at an average rate of a per second; (2) the holding-times of 
calls are mutually independent random variables having the common 
probability density function A(u); and (3) the capacity of the exchange 
is effectively unlimited, and no call is blocked or delayed by lack of 
equipment. This telephone exchange model has been described by J. 
Riordan.° 

As a measure of traffic, it is natural to use the number of calls in prog- 
ress in the exchange. We are thus led to consider a random step-functio1 
of time N(é), defined as the number of calls in progress at time ¢. V(¢ 
fluctuates about an average in a manner depending on the calling-rate, 
a, and the holding-time density, A(1w). 


Il PROOF OF RIORDAN’S FORMULA FOR TRANSITION PROBABILITIES 


Let P»..(t) be the probability that n calls are in progress at ¢ if m 
calls were in progress at 0. Define the generating function of these prob- 
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abilities as 


and let 


a 
f| u) h(a) dx, 
so that the average holding-time, A, is given by 


x 


| flu) du. 


Jo 
Riordan’ ‘has given the following formula for P,,(¢, x7): 
P,(t,z) = [l+ ¢ I)g(t)|" exp f(a lyah{1 — g(t)}}, 


with 


ax 


| ; 
lor exponential holding-time density, this formula had already been 
derived (as the solution of a differential equation) by Palm.’ 

In private communication, J. Riordan has suggested that his proof of 
(1) is incomplete. We therefore give a new proof of (1). 

We seek the generating function of N(¢), conditional on the event 
N(O) m. We obtain it by first computing the joint generating function 
of N(O) and N(t); that is, 


hy" “a e (2) 
The desired conditional generating function is then the coefficient of 


y” in (2), divided by the probability that N (0) m. 
To obtain a formula for (2), we exhaust the interval (— x«, 0) bv 


division into a countable set of disjoint intervals, J, , the n™ having 
length 7,, > 0. Let S, be the sum of the first n lengths, 7; . Let &,(4), 
fort > —S,-1, be the number of those calls which arrive in J, and are 


still in progress at ¢. And let »(¢) be the number of calls arriving during 
(0, 2), ¢ > O, and still in existence at ¢. Then 


N(O) = >> é,(0), 


es 


N(t) nt) + > é.(0, {> 0. (4) 
n 1 


Since calls arriving during disjoint intervals are independent, we know 
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that »(t) is independent of all the é’s, and that &(é) is independent of 
g, 


&,(r) if n # 7. Of course, &,(t) and &,(7) are not independent. It follows 


that if the infinite product converges, then for ¢ > 0 


x 


fa") TT Efy' 


We now compute the terms of the product. If a call originates in in- 
terval J, , it still exists at 0 with probability 


Zier | 
(J), = — | fla + S,—1) du = = fla) du. 
rae a oe 
Hence if & calls arrived in /, , the probability that m of them are still 
in progress at 0 is 


mk ealls arrive in J,.! 


k z bess 
— Ql — Q,)°-", =k: 
m 


Similarly, if a call originates in /, and exists at 0, it also exists at ¢ > O 


with probability 


Therefore 
é,, (0) m and k ealls arrive in J,,} 


[1 + (« — 1)K,]", 


and so 


’ 0) y(t ie ig 
BE {ye x /: ealls arrive in J,} 


The number of calls arriving during /, has a Poisson distribution with 
mean aT’, ; hence 


' 


exp jaT’,(a — 1)} 


exp {aT7,.Q,(y[1 + (a — 1)K,] — 1)}. 


By reasoning like that leading to (6), it can be shown that 
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ot 


| , 
= exp jat(x — I), flu) du? 


J0 


exp fah(x Ill — g(d]}. 


a x 


‘wy aT’,.Q, = | f(w) du = ah, 


os 
>» aT,Q,K, 


flu) du = ahg(t). 


Therefore the infinite product is convergent, and 


exp jah(x Iyfl — gd] 


4 > aT .Q, (yl + (2x 


~ q(t)| t (y 


EB} V0 Vit 
44 ot 
{ 


IA,] — 1)} 


exp jah((a — 1){1 1) + y(a — I)g(t))}. 


Thus the generating function of the joint distribution of V(O) and N(¢) 


into intervals J, . By letting 


is independent of the division of (— «, 0) 
x approach 1 in (8) and finding the coefficient of 7” in the resulting limit, 


we find that 
, ( "(ah)" 


pr{N(O) = m|} 
m! 


The coefficient of y" in (8) itself is 
( (ah). - 
1 - L)g(t)|" exp f(r Lyah{l — g(d]}, 
m! ; 
and so using (9) we find that the required conditional generating function 
V(t), given N(O) m, is given by Riordan’s formula (1). 


ol 


lll THE AUTOCORRELATION 
In terms of N(¢) one can define various stochastic integrals which 
will be characteristic of the process. A simple one which has been ex 
tensively treated in connection with estimating the average traffic is 


ss 
M = [ N(t) dt. 
7 


“0 
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the average of N(¢) over an interval (0, 7’). The chief references in the 
literature on VW are References 3 and 5. If we consider N(‘) during an 
interval (O, 7 + 7), a measure of the coherence of N(t) during this in 
terval, i.e., of the extent to which V(t) hangs together, is given by the 
integral 
— 
U(T,r) =, | NWNG&+7) dt, 
T’ Jo 
depending on values of V(t) taken 7 apart. When the limit ¥(+7) of u as 
T approaches ~« exists, it is usually called the autocorrelation function; 
most statisticians, however, reserve the term “correlation” for suitably 
normalized, dimensionless quantities. It can be shown that this limit 
exists and is the same for almost all V(¢) in the ensemble. It then coin- 
cides with the ensemble average, i.e. 
lim (°(7, 7), almost all V(¢ 


7 x 
E\N(ON(t + 7)! 


The function, y, for the system we are discussing is derived by Riordan, 
and we reproduce his argument for ease of understanding. For equilib 


rium, and b ah, we have 


mqg\t) 


so that 


+ bg(7). 
(Cf.,° p. 1136) 


The limiting value of ¥(7) for + approaching ~« is the square of the 
mean occupancy, b, and the limiting value of ¥(r) for 7 approaching 0 


is the mean square occupancy, b° + b, the second moment of the Pois 


son distribution with mean b. 
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IV THE COVARIANCE AND SPECTRAL DENSITY 


The average value of V(t) isb = ah. One way to study the fluctuations 
of V(t) about its average is by means of the power spectrum used in the 
analysis of noise. (Cf. Rice.*) We resolve the difference [N(t) — b] into 
sinusoidal components of non-negative frequency, and postulate a noise 
current proportional to this difference dissipating power through a unit 
resistance. The spectrum w(f) is then the average power due to frequen- 
cies in the interval (f, f + df). 

More formally, we consider the Fourier integral 

17 
S(f, tip Nae [N (0) ble ii dt, 
and we recall, for completeness, the relationship between S(f, 7) and 
the covariance function, P(r), of [N(é) b|. If 


wif) = lm- 


T—2 


» CO 


wif) =4 R(r) cos 2nfr dr, 


» Ot 


Ri(r) = w(f) cos 2rfr df. 


\t the same time, we have 
R(r) E\(N(t) — b\[N« 
g(r) — 8B 


bq( Tis 


Let X(t) be any stochastic process which is known to be the occupancy 
of a telephone exchange of unlimited capacity, having a probability 
density of holding-time, and subject to Poisson trafic. From the pre- 
ceding result it can be seen that the covariance function of X(/) deter- 
mines the distributions of the X(t) process completely, since 


aid 
gg = = 
drj t—0 


i dR 
f(r) = [ h(u) du = —a” 


de dr 


If the holding-times are bounded by a constant, 4, then readings of 
N(t) taken further apart than / are uncorrelated. In fact, such values 
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are independent, because no call which contributes to N(¢) can survive 
until (¢ + 4), with probability 1. 
Using (11), we see that 


ao 


wif) = 4 | cos 2nfrR(r) dr 


a oO 


th cos 2nrfrg(r) dr 


a a 


cos 2nfr h(u) du dy dr 


vr vy 


a aX 


sin 2nfr h(u) du dr 


“0 “7 


20 
cos 2rfrh(r) ir]. 


Equation (12) expresses the mean square of the frequency spectrum of 
the fluctuations of the trafhe away from the average in terms of the call 
ing-rate and the cosine transform of the holding-time density, h(u 

The calling-rate appears only as a factor, and so does not affect the shape 
of w(f). The function w(f) is what Doob’ (p. 522) calls the “spectral 


density function (real form).”’ 


V EXAMPLE 1. .V(t) MARKOVIAN 
Let the frequency h(w) be negative exponential, so that 


l 
h(u) = =¢ 
h 
where A is the mean holding-time. It is shown in Riordan” p. 1134, 
that V(t) is Markovian if and only if h(w) has the form (13). From page 
ret . l . ° 
523 of Doob we know that the covariance function of a real, stationary 


Markov process (wide sense) has the form 
R(r) R(O)e “", c@ constant 14 


Under the assumption (13), the covariance of N(¢) is 


he ~- 
R(r) bg(7) = j hiu) du dy 


ain 
= du dy 
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in agreement with (14). The spectral density can now be obtained from 


(71) or (22): 16 18 


tbh 


w(t) ——— 
L + 49°f*h? 


This is the same as would be obtained for a Markov process that alter- 
i i; ° r /é l 
nately assumed the values ++7/ah, —~V/ah at the Poisson rate of (2h) 


a ¢. Din04, 920F 
changes of sign per sec. (Cf. Rice’ p. 325.) 


VI EXAMPLE 2. HOLDING-TIME DISTRIBUTED UNIFORMLY IN (a,Q) 


Be ° 
Let h(u) be constantly equal to (8 — a) in the interval (a, 8), and 
constantly 0 elsewhere. Then by (12), 
w(f) =— | cos 2rft dt 
7 Q@ da , 
sin 2xf8 — sin a 


2rf(8B — a) 


Now we see that 


nie 


vy 


so that 


< 


rT 2-68 


is the covariance function of the process N(¢) when holding-time is dis- 
tributed uniformly in (a, 8 

If, formally, we let (8 — a) approach 0 while keeping }(@ + 8) fixed, 
then the holding-times become concentrated in the neighborhood of the 
mean, h; in the limit, as A(w) tends to a singular normal distribution 
with mean, A, and variance zero, we obtain 


( : 
w(f) = — - cos 2fh| 
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us the spectral density function for the V(t) process with constant hold 
ing-time, h = 4(a@ + 8). Similarly, from (15), we note that as the hold- 
ing-times become singularly normal with mean, h, and variance zero, 
the covariance function becomes 
a(h — 7) Srsh 
R(7) 
0 


We can express (16) as 


o fsin rfh\ 
w(f) = 2ah — _ 
ath 


and note that this is exactly like the power spectrum of a random tele- 


graph wave constructed by choosing values ++7/ ah, —Vah with equal 


probability and independently for each interval of length, h. (Cf. Rice, 
page 327.) 
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Silicon power rectifiers have been made which have reverse breakdown volt 
ages as high as 2,000 volts and forward characteristics comparable to those 
obtatned tn much lower voltage devices. It 7s shown that the magnitude and 
temperature de p ndence of the currents can be explained on the basis of 
space-charge generated current with a trapping level 0.5 eV below the con 
duction band or above the vale nce band. The breakdown voltage of a P 
IN~ diode is computed from avalanche multiplication theory and is shown 
to be a function of the width of the nearly intrinsic region. A simple diffusion 
process 1s evaluated and shown to be ade q date for diode fabrication The 
charact rostics of devices fabr cate d= from high-re sistipity compet Sale 
floating-zone refined, and gold-doped silicon are presented. The surface 


tation to high inverse voltage rectifiers 1s discussed 


I INTRODUCTION 


The desire for high voltage rectifiers in the electronic industry has 
pushed the peak inverse voltage of solid state rectifiers to higher and 
higher values. The purpose of this paper is to present some ol the con 
siderations necessary in designing a device with a high inverse voltag 
and an excellent forward characteristic. In many cases the device charac 
teristics are predictable. Conversely, high voltage diodes are excellent 
tools for studying many solid state phenomena. 

It has been shown! that it is possible by the use of the conductivity 
modulation principle to separate the design of the forward current-volt 
age characteristic from the reverse current-voltage characteristic of a 
silicon p-n junction rectifier. Units have been fabricated by the diffusion 
of boron and phosphorus into high resistivity material, that have reverse 
breakdown voltages in the range of 1,000 to 2,000 volts. 

The reverse currents are of the order of a microampere per square cen 
timeter at room temperature and increase approximately as the square 


QED 
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root of the applied voltage. The magnitude, voltage dependence, and 
temperature dependence of the reverse currents can be explained as due 
to space-charge generated current? with a trapping level 0.5 eV from 
either the conduction or valence band. These effects will be discussed 
in Section IT. 

In Section III the breakdown voltage and its dependence on the re- 
sistivity and width of the high resistivity region of the rectifier will be 
considered. 

In the next section the forward current is discussed and explained by 
considering both a space-charge region generated current and a diffu- 
sion current that takes into account high levels of minority carrier in- 
jection. 

Device processing information is given in Section V, together with an 
evaluation of different sources of high resistivity silicon. The devices 
to be discussed in this paper have been processed with high resistivity 
p-type material, although some devices have been made with n-type ma- 
terial. 

linally, a discussion of some surface problems associated with high 
voltage rectifiers is given in Section VI. 

Although this paper is entitled “High-Voltage Conductivity-Modu- 
lated Silicon Rectifier’, the theoretical arguments are applicable to all 
semiconductor diodes. However, the experimental results have been 
limited by considering only high voltage diodes. 


II REVERSE CURRENT-VOLTAGE CHARACTERISTI( 
2.1 Theory 


The simple theory’ for a p-n junction yields an expression for the re- 
verse saturation current density (/») which is: 


r (2+) (2*)" ; 
ih p 8 ) Pp, = . iD 


where g is the electron charge, n, is the equilibrium electron density in 
p-type material, p, is the equilibrium hole density in n-type material, 
D, and D, are the diffusion constants for electrons and holes, and r, 
and 7, are the minority carrier lifetimes for electrons and holes. 

When reasonable numbers are substituted into (2-1), 7) at room tem- 


. . 10 : ra. * 
perature is of the order of 10 ~ amperes per square centimeter. This 
quantity doubles with every increase of 4° C. The theory also contains 


no voltage dependence of this current. Even when breakdown multi- 


plication’ is taken into account, there is essentially no voltage dependence 
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at voltages less than half the breakdown voltage. The magnitude and 
temperature and voltage dependences of measured diodes do not agree 
with these theoretical values at room temperatures. 

Recently, Pell” has shown that the reverse currents at low temperatures 
in germanium, and at room temperatures in silicon, are dominated by 
space-charge generated current. The space-charge generated current 


density (/,..) is given by 
le =qWGM, 2-2 


where W is the width of space-charge region, G is the generation rate of 
hole-electron pairs in the space-charge region, and VW is the breakdown 
multiplication (JJ ~ | except near the breakdown voltage). @ is given 
by” 
( nipr . 

P1 
where n; and p; are the densities of electrons and holes respectively if 
the Fermi levels were at the energy level of the recombination centers, 
and r,9 and ry are the minority carrier lifetimes of electrons and holes 
respectively in heavily doped p-type and n-type silicon. This expression 


assumes constant generation over the space-charge region. Thus, 


ny NV. exp of (V oo n, exp B(V 


Pr N eXP op (V - n; exp — B(V, rok 2-4b 
where V, is the recombination level above the valence band edge V 
V’; is the midband intrinsic level, 8 = g/kT, N. and N, are the effective 
densities of states in the conduction and valence bands & 2.4 « 10" 
(T/300)°, V. is the conduction band edge, / is the Boltzmann’s constant. 
and 7' is the absolute temperature. 


> 


Substituting (2-4) into (2-3), one obtains: 


Ny l 
os : 
“V 707 eosh | B(V, — Vi) 


+ 3 In 

For the diffused silicon junctions under consideration, it has been found 
that rr» equals 1.2 X 10°” seconds and 7,» equals 0.4 X 10 | 
Also, n,; ae glk bs Mica 


seconds. 


and V 0.54 volts. Using these 
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numbers, (2-6) becomes: 


r 
9 le 
seetallcbiinhs (<0) 
cosh | 38.62 (+) CY, 0.55 
| T , 


Gr 


where 


tf cine a cosh [38.62(V, 0.54) 
( sa) : 


500 : 
cosh 8.02 ( T ) (V, — 0.54) 0.55 


In (2-6b), Gyoo(V,) is the generation rate for a recombination level at 
V’. equal to 300° K, and f(7, V,) is the temperature variation of G for a 
recombination level at V, normalized to 300° Ix 

Curves of f(7', V,) are given in Fig. 1 for several values of Vi. with a 
curve g(7) which is the temperature variation of the reverse saturation 
current (/)). Table I gives values for Goo for various V 

In the reverse biased diffused junctions made with high resistiv itv ma 
terial, the Junction may be considered abrupt. Therefore, the width (IV 
of the space-charge region, when the junction is reverse biased to a volt 


age V, is given by 


, kV \' — 
i (s-) 3. x 10 °(Vp,)° ~ em 2-8) 


where the units after the first equal sign are electrostatic, and « is the di 
electric constant. In the second expression, V is in volts, and p,, the base 


material resistivity, expressed in ohm-centimeters. Thus, 
/ 1X 10“ Gyo(V f(T, V)[Vp,|'” amperes-em ”. (2-9) 


It is seen that /,. varies theoretically as the square root of the reverse 
voltage for values of V less than } V, , the breakdown voltage, in which 
range avalanche multiplication is negligible. The quantity /,. varies in- 
versely with V,'* and will be large for high voltage devices with small 


N,. The /,, at 300° K for a rectifier with 40 ohm-centimeter base ma 
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terial and a reverse bias of 100 volts is given in Table I as a function of 
V’,. The numbers compare with 8 X 10°" ampere per square centimeter 
for J) . Thus, from diode measurements at room tempersture and above, 


one could not observe V7 less than 0.3 eV from either the conduction or 


valence band. In fact, from a measurement of the temperature depend- 


ence of the reverse currents, one can determine only the recombination 
level lying closest to the center of the forbidden band. This can be seen 
more clearly from the following argument: There wil] be a contribution 
to the reverse current from the diffusion current J,300g(7) which varies 
with temperature as g(7'). There will be contributions to the reverse cur- 


108 





5 











Fig. 1 The temperature variation of the generation rate, f(7', V,), for several 
values of the recombination level, V 
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TABLE I — VALUES OF G AND Space CHARGE GENERATED 
CURRENT AT 300° K ror Vartous VALUES 
OF THE TRAPPING LEVEL V, 
300° K E, = 1.08eV np = 3 X 10” em-* 
tno = 1.2 XK 10-* sec ton = 0.4 X 10-* sec 


} 


Tee (V = 100 volts, p = 40 ohm-cm) 


olts above Valence Band _ 4 
Vr Vo eee eee ‘ microamperes/cm? 


.10 
.20 
.30 
.40 
50 
54 
.58 
.68 
18 
RS 
98 


1.88 X 107-7 
9.03 X 10 
4.29 X 10 
2.06 X 10 
0.96 

3.43 
2 
6 
] 
2 
3 


> RS Or 


S 


.58 

oe. 
.26 X 
.69 X 
75 X 


ee ee oe 
x KKK KK KK KK OK 


rent by the individual trapping centers given by J,300(Vf(7, V+), where 
I scaoo(V) is the current at 300° K due to generation at recombination 
centers located at the level V., and f(7’, V,) is the temperature variation 
of the generation rate. Thus, the total reverse current is given by 


Treverse = Losng(T) + D Tecso(V-)f(T, Vr), (2-10) 
Vr 


where the summation is over all recombination levels. The relative cur- 
rents at 300° K are given in Table I. The greatest contribution at 300° K 
is due to the level nearest the center of the forbidden band. As the tem- 
perature increases, all the terms under the summation sign approach each 
other. Before a second recombination level contributes significantly to 
the reverse current, however, the saturation current will have become 
the most important component. 


2.2 Experimental Results 


To evaluate the theory for the reverse currents in silicon N "P junctions, 
careful measurements were made on five typical units for the reverse 
current-voltage characteristics at various temperatures from 300° Kx 
to 435° K. The curves were taken with a X-Y recorder. The voltage 
ranged from 0 to 200 volts so that multiplication effects were completely 
negligible. 

From the recorded data, curves of /, versus 1,000/T° K were plotted 
for V -10, —40, and — 160 volts. The set of curves for diode No. 3 
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Fig. 2 — The temperature variation of ‘‘reverse current’’ for a typical diode 
at —10, —40, and —160 volts 


is given in Fig. 2. The slope of these curves indicates that the recombina- 
tion level lies near 0.5 eV below the conduction band or above the val- 


ence band. The junction area of this device is 0.015 em‘; thus, the eur- 


rent density at 300° K and at — 100 volts is 4.4 microamperes per square 
centimeter. This compares with the order of one microampere as listed 
in Table I. This suggests that the (rnotn0)’ is overestimated. The agree- 
ment of this measurement with theory is reasonable. 

The voltage variation of the reverse currents does not agree with 
theory as well as the magnitude and temperature dependence. The ex- 
perimental results give, as the voltage dependence, an expression: 


well 


where N equals 2.9. This compares with the theoretical value of 
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Some of this discrepancy can be attributed to the fact that the junction 
is not truly an abrupt junction. A ‘‘graded” junction would yield N 
equals 3. Measurements of capacitance versus voltage, which essentially 
measure the width of the space-charge region, yield N equals 2.4. Thus, 
these devices in the relatively low voltage range still have some contri- 
bution from the gradation of the diffused junction. 

The highest temperature points in Fig. 2 deviate above the straight 
lines. This deviation can be attributed to the onset of the contribution 
from the /) component. Calculations indicate that, by 7 220° C, the 
contribution to the reverse current by the space-charge current is equaled 
by the saturation current and that, by 7 320° C, the space-charge- 
generated current is negligible compared to the saturation current 


Iii BREAKDOWN VOLTAGE OF PN AND PIN JUNCTIONS 
S.1 Theory 


It has been demonstrated that, in germanium’’” and silicon, reverse 
biased junctions breakdown as a result of a solid state analogue of the 
Townsend 8 Avalanche Theory. Multiplication and breakdown occur 
when electrons or holes are accelerated to energies sufficient to create 
hole-electron pairs by collisions with valence electrons. The breakdown 


phenomena in silicon for graded and step junctions has been previously 


considered.”: '” Depending on the impurity distribution, the field in the 


junction will be a function of distance and will have a maximum value in 
the region of zero net impurity concentration. The breakdown voltage is 
a critical function of the space-charge distribution. 

In this section the existing multiplication theory is extended to the 
case of PIN junctions. It is shown that relatively wide intrinsic regions 
are required to obtain breakdown voltages greater than 1000 volts 

Fig. 3 is a plot of the impurity, charge, and field distributions in PIN 
and PrN junctions. Fig. 3(a) schematically illustrates the geometry of 
the three region devices considered, and Fig. 3(b) is a plot of the impurity 
distribution. In this analysis step junctions will be assumed. For the 
PIN junction there are no uncompensated impurities in the intrinsic 
region, and no net charge. At low reverse voltage, the field will sweep 
through the intrinsic layer and will increase with increasing reverse bias 
until the breakdown field is reached. 

Absolutely intrinsic material is not yet available, and devices are 
made from high resistivity r-type material. In this class of devices there 
is some uncompensated impurity and charge in the center region. The 
field will have a maximum value at the NV 7 junction and will decrease 
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with increasing distance into the 7 region. At sufficiently high reverse 
bias the field may sweep into the P~ region. 


Breakdown in silicon’ is a multiplicative process described by 


WwW 


i I nee = 


where .V is the multiplication factor, W is the space-charge width, and 
a, is the rate of ionization which is a strong function of the field in the 
junction. For a PIN structure, the field is constant, at breakdown 
M approaches ~, and 


ay W ; (3-2 


The ionization rate at breakdown is then a simple function of the width 
of the intrinsic region. McKay’ and Wolf” have considered a, as a fune- 






































Impurity, charge and field distribution in PIN and PxN junctions 
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tion of the field. If a; is fixed, the field at breakdown can be determined. 


“VW 


The breakdown voltage is Edu. 


lig. 4 is a plot of breakdown voltage as a tunction of space-charge 
width for PN and PIN diodes. The PIN values are calculated; the PN 
data is previously unpublished data supplied by K. G. McKay 

Some interesting observations can be made from Fig. 4: 

1. The plot of breakdown voltage versus barrier width for a PN step 
junction assumes that the space-charge region does not extend through 
the high resistivity side of the junction. For this class of junctions the 
breakdown voltage is determined by the impurity concentration as shown 
in Fig. 5. The plot of breakdown versus space-charge width for a PIN 
diode assumes that the space-charge region extends from the P to N 
region at very low bias, and that it is limited by the width of the I re- 
gion. If a constant field is assumed in the I region, the breakdown volt- 
age is a function of the barrier width 

2. Although the space-charge region can reach through the I region 
at low bias, the avalanche breakdown voltage is a function of the width 
of the I region. 

3. For the devices considered here with w regions in the order of 10 
cm, the maximum breakdown voltage is in the order of 2,000 volts 

















Fig. 4 Breakdown voltage as a function of barrier width for PN and PIN 
junctions 
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Fig. 5 Breakdown voltage versus impurity concentration for silicon step 
Junctions. 


Oo. Expe rument 


Fig. 5 is a plot of breakdown voltage versus impurity concentration 
for silicon step junctions. The plot above 300 volts is extrapolated from 
the data of Miller” and Wilson.” 

Capacity data, discussed in Section V, indicates that many devices 
show body breakdown. A few rectifiers break down at voltages as high 
as 2,000 volts. In many high voltage devices the breakdown voltage is not 


} 


limited by geometry but by surface problems 


IV FORWARD CURRENT-VOLTAGE CHARACTERISTIC 
4.1 Theory 
It will be shown in this section that the forward current-voltage char 
acteristic as well as the reverse characteristic can be completely explained 
by considering both a space-charge region generated current and a diffu 
sion current. The diffusion current component must also take into con 
sideration the effect of high injection levels of minority carriers. 
According to the Shockley-Read? theory, the rate of recombination, 


UU’, of holes and electrons in a semiconductor is given by: 


. yw i] 
G; y 
Tp -- Ny) + 
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where p, and n are the instantaneous concentrations of holes and elec- 
trons, respectively. When a PN junction is forward biased, holes and 
electrons are injected into the space-charge region which has been re- 
duced in width. Some of these carriers diffuse through the space charge 
region and give rise to the normal diffusion current when the excess 
minority carriers recombine with majority carriers in field free regions. 
The other carriers recombine according to (4-1) in the space-charge re- 
gion giving rise to what is called the space-charge generated current. In 
the reverse biased junction, the current is due to carriers generated in the 
space-charge region; whereas, in the forward biased junction, the current 
is due to recombination of carriers. The quantity U is large in the space- 
charge region since both p and n are large in this region. In the field free 
regions, however, one of these quantities is usually small and the product 
deviates only slightly from n,’. 

The space-charge generated current, /,., is given approximately by: 


| 
sinh 8 

2qWn 2 f(b). (4-2) 
(trot po)? B(V gs te 
where Vy is the built-in potential of the junction, and f(b) is discussed 
in Reference 12 and is approximately 1.5 for recombination centers near 
the intrinsic level as is the case for the diodes under consideration. For 
shallower recombination levels the function f(b) is much smaller and de- 

pends strongly upon the forward applied voltage. 
For the forward-biased junction, the space-charge region is narrow, 


the concentration gradient can be considered linear and W is given by 


. . a 4 
the following expression: 


AV ly 
W 35 X 10°( mate em, (4-3) 


where V junction IS the total potential across the Junction in volts and a 


is the concentration gradient at the junction in em *. These are given 
by: 
ilt 
kT /q \n (NaNv/n; ) 
0.792 — |} 
Co 


Als » < = € 
sa VrDt 


for diffused junctions, 


where ( surface concentration of diffusant 3 x 10” em 
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D diffusion constant 3 10” em’*/see, 


diffusion time 5.7 X 10° see, 


x; junction depth below surface 0.003 em 
When these numbers are substituted into the equations, at 500 
W = 9.25 < 10 “(0.792 — V)"* em. 
l’or the diodes under consideration: 
1.2 x 10 ’ see, r, 0.4 x 10 
When these expressions are substituted into (4-2 one obtains at 300° C 


; sinh 19.51V 
10 = —— amp/em’. 
0.792 — J 
In order to fit the experimental data, it is necessary to multiply (4-7 
by a factor of 5. This may be due to an overestimation of (7,97 There 
fore, the equation which shall be used in the remainder of this sectior 
will be: 
6 sinh 19.351 
10 = amp,/cm 
(0.792 J 
A plot of this expression is given in Fig. 6 
The normal diffusion current for low level diffusion,’ Ip 


Ip: Lo e | 


where /) is given by (2-1). /) for the diodes under discussion is approxi 
mately 8 X 10 is ampere em at 300° Kk. When the injected minority 
carrier density approaches the equilibrium majority carrier density, the 
form of (4-9) changes. The high injection level diffusion current, /px , 


is then given by 
lou Lpuole’ OF 1-10 


where /pyo equals gnjs/7, and s equals the width of the high resistivity 
region. For the diodes under discussion, /pyo is approximately 2 X 10— 
amperes/em™ at 300° K. A current-voltage plot of these currents at 
300° Ik for V 0.50 is given in Fig. 6 together with their sum. It can 


be observed that the resulting characteristic starts with slope ol q\ ih 


and bends over to a slope of gV/2kT near 0.10 volt. The slope increases 


again to near qV//kT at 0.35 volts and decreases once more to gV /2k7 


above 0.40 volts giving a bump to the over-all characteristic 
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Next, consider the temperature dependence of the coefficients of the 
forward current components 


lo ~ nT), (4-1lla) 
pity 


n 7.(T) ~n(T)’, (4-1 1b) 
~ mitt, ~ nT). (4-11¢) 


The largest variation of these coefficients is due to the variation ot 





oe oe en 


two components of current for a forward biased junction 





CONDUCTIVITY-MODULATED SILICON RECTIFIER SY 


n,; with temperature. Fig. 7 gives a plot of this variation. The tempera 
ture variations of the other parameters are all small compared to that 
of n;. Thus, as in the case of the reverse currents, at sufficiently high 
temperatures, the diffusion current makes the more important contribu 
tion. 

In the case of the forward current, /,. is relatively insensitive to the 
distribution of impurities; therefore, the results of this section are im 
portant for all forward-biased diodes. In high-voltage diodes, to keep the 
resistive voltage drop small, it is necessary to maintain high minority 


carrier lifetime in the center region. The diffusion length of injected 


| 
| 
| 





The variation of n; with temperature 
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minority carriers should be the order of or larger than the center region 
width. 


1.2 Kxeperimental Results 


The forward characteristics of the five typical high voltage rectifiers 
mentioned in Section 2.2 were measured with a X-Y recorder, and all 
showed similar shapes. Diode No. 3 will be discussed in detail in this sec- 
tion. 

The forward current-voltage characteristic was measured at three 


g 
temperatures: 220° Kk, 300° K, and 375° K. Measurements below cur- 


1072 7 
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Fig. 8 The calculated and observed current-voltage characteristic of a for- 
ward biased junction at 300° K. 
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rents of one microampere were made only at 300° K. Currents above 10 


milliamperes were not measured since internal power losses would cause 
temperature variations. The unit has a junction area of 0.015 em’, june- 
tion lifetime of 4 microseconds at 300° K, S equal to 0.008 em, and N, 
equal to 3 X 10 em °. When these numbers are substituted into the 


expressions for the coefficients, one obtains, at 300° K, 


| 1.4 X 10 ° amperes, 
if 1.2 x 10“ amperes, 
lon 3.0 X 10 ° amperes 


Fig. 8 shows a semilogarithmic plot of the current-voltage characteristic 
at 300° K over a range of 63 decades. The circles represent measured 
points and the solid line is the theoretical curve. 

Using the variation of n; given in Fig. 8 and the temperature variations 
as given in (4-11), one can obtain the coefficients for any temperature 
This has been done for two temperatures, 220° K and 375° K. Figs. 9 
and 10 show the theoretical and experimental plots at 375° K and 220° K 


1072 





























Fig. 9 The calculated and observed current-voltage characteristic of 
ward biased junction at 375° K 
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respectively. The circles represent measured points and the solid lines 
are the calculated theoretical curves. It is observed that the fit in Fig. 
9 is quite good; whereas, the fit in ig. 10 is not as good as at the other 


temperatures. However, even this figure shows good qualitative agree- 


ment of the deviation from a straight line. Some of the factor of two dis- 
crepancy 1n | ig. 10 can be ascribed to the temperature Variation of the 
other parameters, and some to a possible error in the measurement of 
temperature which would be reflected in the value of x 

It should be noted that at all temperatures the JR drop in the high 
resistivity region is not observable to the limits of the experimental 
measurements of forward current, 10 milliamperes. This is due to the 
fact that the region has been conductivity modulated by the forward 
current. This requires a sufficient minority carrier lifetime in the region 
so that most of the injected carriers diffuse across the region before re 


combining. Such lifetimes can be maintained in diffused junctions’ 











Fig. 10 The calculated and observed current-voltage characteristic of a 
forward biased junction at 220° kK 
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to permit the high resistivity region to be at least as wide as 0.025 cen 
timeters. Thus even in high voltage rectifiers it is still possible to design 


the forward and reverse current voltage characteristics independently 


V DEVICE PROCESSING 
5.1 Silicon Material 


Fig. 5 shows that step junctions which break down at over a thou 
sand volts must have a background impurity concentration < 10° 
atoms/cm*. The highest grade commercial semiconductor silicon has 
5 X 10" impurities/em* (20-50 2 em P type). This material must be 
processed to reduce the impurity level. To date, high voltage devices 
have been processed Irom four types ol high resistivity material: float 
ing Zone refined, compensated, vold diffused, and horizontal zone renil ed 


silicon. 


Some silicon was prepared by adding N-type impurities to reduce 
Np Na < lO” Maintaining this delicate balance in material where 
Vp =~ N, is difficult. The boron is relatively uniformly distributed since 


the distribution constant is close to unity. N-type impurities are less 


uniformly distributed in the crystal since the distribution constants are 
considerably less than unity. High resistivity compensated silicon is full 
of N- and P-region striations. The units processed from this material 
generally had poor electrical characteristics 

Table IT isa typical contour of a compensated crystal. The resistiy ity 
varies around the crystal and changes along the length of the ervstal 
\t the bottom of the cry stal the resistivity goes through a m: 
The tail end is converted from P to N type 


A number of devices have been fabricated from silicon pro essed with 


TABLE I] \ TypicaL Conrour or A Hicu 
RESISTIVITY COMPENSATED CRYSTAL 


Crystal A-161, Oriented 111, Rotated 2 RPM 
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TABLE III] THE CHARACTERISTICS OF SomME HicH VOLTAGE 
RECTIFIERS PROCESSED FROM GOLD DIFFUSED 
AND ZONE REFINED SILICON 


1 ma 


120. = 500 : 
200 600 F os Gold diffusion 
600 1000 h , p =~ 16,000 2 em 
500 | 1000 ; 


1200 | 1500 x ; <i 
70 | 300; 800) 2. 4, <1 Floating zone refined 
16 120 700 | 3.i p ~ 6,000 2 em 

320 100 = 800 : 


! These units have an area ~ 107% em?. 
> This is the reverse voltage at which these units pass the indicated current 
This is the forward bias at which the units pass the indicated current. 
‘ This is the lifetime measured at 30 milliamperes forward current by the pulse 
injected technique. The lifetime did not seem very sensitive to small variations 
in injected current 


the floating zone apparatus.” This technique removes impurities from 
molten silicon by treatment with hydrogen containing water vapor. The 
material obtained from this process has an impurity level in the range 
of 10°" to 5 X 10" acceptors, cm (2.000 to 16.000 2 em P type). 

Table IIL gives the characteristics of some of the better diodes made 
from such floating zone silicon. The reverse currentsare larger than that pre- 
dicted by (2.10). The lifetime at high injection isin the order of 1 usec. 

N-type silicon with a resistivity range of 10 to 30 Q em was diffused 
with gold at 1,200° for sixteen hours. With this diffusion program the 
gold is uniformly distributed in the material.'* The resistivity after gold 
diffusion was in the range of 2,000 to 15,000 2 em. The characteristics 
of several devices processed from this material are given in Table III. 
This technique has many attractive features; however, additional work 
was not done because the lifetime in the diffused material was consist- 
ently lower than that required for conductivity modulation. 

One successful purification technique is horizontal zone refining of 
silicon in a quartz boat. With the number of passes used, the background 


acceptor concentration is observed to be in the order of » y 4 10"° to 10" 
(100-1,000 2 em P type). Most of the devices reported in this paper are 


fabricated from this material. 
Capacity data for devices fabricated from various types of high re- 
sistivity silicon is shown in Fig. 11. The plot shows that the high resistivi 
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Fig. 11 Capacity/cm? for high voltage rectifiers processed from various type 
of high resistivity material 


ties measured by the four point probe before diffusion are indicative of 
the impurity level after processing The water vapor floating zone re 
fined material has an impurity level of 10° acceptors’em*; the other 
material is in the range of 10° to 10°. The breakdown voltage line was 
caleulated from the data in Fig. 5. 


5.2 Diffusion 


In this section some of the practical difficulties observed in utilizing 
the diffusion technique will be considered. In the fabrication of transis 


tors close geometry control is necessary 1m order to obtain the desired 


device characteristic. It has been shown that in conductiy ty modulated 


rectifiers the only geometry requirement is that the width of the center 





QOG rHE BELL SYSTEM TECHNICAL JOURNAL, JULY 1957 


region be less than the diffusion length of the minority carrier. High sur- 
face concentration of diffusant is desirable since this facilitates the con- 


tact problem. This suggests that the diffusion system can be much less 


involved than that required for diffused transistors.’? Some of the data 


presented in this section will show that the open tube diffusion tech 
nique — can lead to variations in diffusion parameters. 

The diffusion of impurities into silicon is complicated by Variations 1 
the boundary conditions at the surface. Frosch” has shown that surface 


concentrations can be varied over six decades. 


1 Device Diffusion Theory 


Several important impurity distributions have been considered ‘. Two 
distributions are important in the open tube process: 

1. Error Function Comple ment, ERFC, Distribution or Infinite Diffus- 
ant Source. If the diffusant is deposited on the silicon and serves as an 
infinite source, the added impurities will have an erfe distribution. lor 
one diffusant and a fixed diffusion program this distribution will result 
in the deepest penetrations and smallest sheet resistances of all possible 
distributions. The sheet resistance is a measure of the total number of 
added impurities. The data presented later indicates that the added im 
purities frequently have an erfe distribution. 

2. Gaussian and Modified Gaussian Distribution. A number of impurity 
atoms enters the solid, and a surface barrier builds up with time which 
prevents additional atoms from entering Initially, the diffusant is as- 
sumed to be present in an infinitely thin layer at the surface with diffu- 
sion into or out of the material possible In the range of silicon doping 
levels and surface concentrations used, a Gaussian, modified Gaussian 
or erfe distribution for a given diffusion program lead to approximately 
equal junction depths. 

The sheet resistance and the diffusion depth have been related” to 
the surface concentration for an erfe distribution. If the distribution is 
Gaussian instead of erfc, then for the same value of sheet resistance and 
diffusion depth the surface concentration should be reduced by one-third. 
Since the sheet resistance is related to the total number of impurities 
through a mobility term, quantitative interpretation of the data for any 
case other than erfe or Gaussian distributions would be difficult. 


ID perume ntal Results 


The sheet resistance was measured by the four-point probe method, 
10 


and the diffusion depths, by angle lapping and staining.” Surface con- 
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Fig. 12 Distribution of surface concentration of 28 P.O; diff 
open tube deposition technique 


centrations were calculated assuming an erfe distribution.” All th 
diffusions are on lapped silicon surfaces in the temperature range of 1,200 
to 1,300° C 

lig. 12 shows the distribution of surface concentration of 28 P.O 
diffusions by the open tube process. ” The surface concentrations Vary 
from 10° to 5 X 10° atoms/cm’. These values are about a decade lower 
than the closed tube values of surface concentrations reported by Fuller. 

The measured diffusion depths were in the order of 2 x 10 to 


5X 10 “em. Fig. 13 shows the distribution of diffusion depths normalized 


with the calculated diffusion depth as unity. The diffusion depths were 


calculated from the measured surface concentration assuming an erfe 
distribution. 

The observed variation in diffusion depth is difficult to explain. Som 
of the possibilities which have been considered are: 

1. The diffusion temperature from lot to lot would have to be from 0 
to 50 degrees below the expected value to explain the Variations lis 
crepancies this large have not been observed 

2. One impurity distribution which may explain some of the results 
is a modified Gaussian with considerable out diffusion. There are some 
runs with high sheet resistance and diffusion depths which are consistent 
with this picture. ( renerally the sheet resistances are so small that there 
could not be much out diffusion. 

3. Some workers have suggested the possibility of the diffusion con 


stant being a function of the surface concentration. Fig. 13 does not 
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indicate any correlation between surface concentration and diffusion 
constant. 

The variations in diffusion process control have not been observed to 
effect the production of rectifiers. If better geometry control is necessary, 
more sophisticated diffusion techniques are required. 


VI PULSE PROPERTIES AND RELIABILITY 


Important considerations in all diode applications are the pulse prop- 
erties and reliability in operation. In this section some problems which 
are associated with avalanche breakdown are described and the results 
related to recent work on surface and body breakdown. 


Hy 


6.1 Theory 


4 20 . . “1 . . . 
Several workers” have considered the possibility of a negative resist- 


ance in the avalanche region for reverse biased junctions in which one 
side is either intrinsic or so weakly doped that the space charge of the 
carriers cannot be neglected. A negative resistance might be observed 
at very high current densities in an N*z junction. 

One possible source of a negative resistance would be a large tempera- 


ture rise due to current concentration at a few points instead of a uni- 
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Fig. 13 Distribution of diffusion depths for diffusion by the open tube depo 
sition technique 
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form flow through the Junctions. This case is of particular significance 
in high voltage rectifiers where small reverse currents result in relatively 
large power. It has been pointed out in Sections 3.2 and 5.1 that body 
avalanche breakdown is frequent ‘* not observed in these devices 
Avalanche breakdown current 1. silicon’ is carried by discrete pulses 
of about 50 wa at their onset and increasing with increasing current to 
about 100 wa. Approximate calculations” show that the ionizing regions 
of these microplasmas are about 500 A in extent, have a current density 
~ 2 X 10° amp/em’, and have a net space-charge density ~ 10°°/em 
These pulses for junctions with Ey,4x less than 500 kv/em appear to be 
independent of junction width and built-in space-charge. Rose considers 
the statistical problem associated with a large number of pulses and pre- 
sents a picture which is consistent with most of the experimental data 
He calculates the temperature rise, assuming the avalanche power is 


| X 10 ~ watts and is dissipated uniformly in a sphere. The maximum 


temperature rise for a cluster of two or three pulses is in the order of 
25° C. For the picture Rose presents, the temperature rise due to the 
microplasma should be relatively insensitive to the breakdown voltage 
Thermal collapse of rectification, i.e., increase of temperature until the 
silicon is intrinsic, will probably not oceur in the region otf avalanche 
multiplication. Two important conclusions can be obtained: 

1. Avalanche breakdown should occur as a random process with a uni 
form probability over the junction. Large temperature rises due to a 
breakdown of microplasma will probably not occur since the resulting 
temperature rise would cause the breakdown voltage in that spot to 
increase. The power is dissipated throughout the path of the current 
pulse in the space-charge region. 

2. A thermal effect in silicon due to heating by the small plasma has a 


1 


very short time constant of the order of 10 ~ seconds.” It is not possible 
to separate a thermal effect of this type by reducing the pulse width 
The heating and cooling time is short compared to the pulse time in these 
experiments. 

The pulse properties of a junction would be quite different if the break 
down occurred at one spot instead of many spots distributed over the 
junction. Breakdown at a single spot on the surface has been observed 


6.2 Experimental Results 


Many rectifiers were given a voltage pulse which carried them into 
breakdown. There was a wide distribution of V-I characteristics. Many 
diodes did not show a negative resistance up to the maximum instan 
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taneous power the pulser could deliver, 5kw. These diodes are not con- 
sidered in the subsequent analysis. 


The diodes were subjected to 50 usee triangular voltage pulses which 
would send them into breakdown. Variations in pulse conditions did not 


effect the I-V characteristic until large pulses destructively damaged the 
unit. 


Fig. 14 isa sketch of a typical V-I characteristic and Fig. 15, shows the 
voltage-versus-time characteristic for a diode with a negative resistance. 
The V-I curve can be broken into four regions 


Ve 











Fig. 14 


\ typical V-I characteristic for a diode in which a negative resistance 
IS observed 


i 











Fig. 15 A typical V-T characteristic for a diode in which a negative re 
sistance is observed 
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1. A high impedance state before the breakdown voltage is reached 
2. A current required to turn on the negative resistance; this cur 


rent varies from 10 ° to 1 amp. 
5. The transition to a low impedance state. 


t. The low impedance region in which the current is probably limited 


by the circuit impedance. 

The V-T curve can be broken in four regions: 

|. The time it takes the pulse to reach the breakdown voltage. 

2. The time the diode can maintain the breakdown voltage less than 
l usec. This is beyond the resolution of the oscilloscope. 

3. The time required to fall to the low voltage (low impedance) state, 
is less than 1 usec. 

1. The remainder of the pulse in the low voltage state. 

I'ig. 16 is a plot showing the current and voltage required to turn on : 


negative resistance in several power rectifiers (area ~ 10 ~ cm To 




















Fig. 16 — Current and voltage required to turn on a negative resis 
several power rectifiers (A ~ 10-2 em? 
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consider the spread of breakdown voltage, the data was normalized to 
the instantaneous power required to turn on a negative resistance. This 
turn-on power was the turn-on power multiplied by the voltage. 

Fig. 17 is a plot of the distributions of turn-on power for the rectifiers 
which had a negative resistance plotted as a log normal distribution on 
probability paper. The median value for the turn-on power is 1.2 watts. 
Kighty per cent of these diodes went into a negative resistance condition 
at powers between 0.1 and 10 watts. Many diodes could dissipate several 
kilowatts with no negative resistance. These were not included. 

Ixperiments show that devices which show surface breakdown will 
collapse at power levels which are orders of magnitude below that ob- 
served for devices in which body breakdown is observed. 


The picture is more cloudy with smaller area rectifiers (area ~ 10 


em). In these devices it was not possible to predict the pulse properties 























Fig. 17 — The distribution of turn-on power for rectifiers (4 ~ 10°? em?) in 
which a negative resistance is observed plotted as a log normal distribution on 
probability paper 
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of the device from the reverse I-V characteristic. This may be attributed 
to the decrease in power capabilities of the body breakdown process in 
the smaller devices. This also suggests that the smaller devices have a 
less severe surface problem. 


The distribution of turn-on power for a few hundred small area recti- 


fiers (A ~ 10 ° em’) is shown in Fig. 18. The median of the distribution 


occurs at 40 watts. Eighty percent of the units will show a negative re- 


> 


sistance when pulsed at power levels between 3 and 500 watts. 


VII CONCLUSION 


High voltage rectifiers have been fabricated using several sources of 
£ £ g 
high resistivity material employing an uncomplicated diffusion process. 











Fig. 18 The distribution of turn-on power for small area rectifiers (A ~ 10 
em?) plotted as a log normal distribution on probability paper 
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The most consistent results were obtained using horizontal zone refined 
silicon. The open tube diffusion technique has sufficient control to satisfy 


the fabrication requirements. 


The magnitudes, voltage and temperature dependences of both the 
forward and reverse currents of silicon rectifiers can be explained by in- 
cluding a recombination level near the middle of the forbidden energy 
gap. Design equations for the forward and reverse characteristic of a 
diode are presented for several important cases. The breakdown voltage 
of the high voltage devices was shown to be a function of the width of 
the high resistivity region. 

One unsolved problem is the surface limitation of breakdown voltage 
and reverse currents. This has been observed to decrease the breakdown 


voltages and increase the reverse currents to undesirable levels. 
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Coincidences in Poisson Patterns 


By Ek. N. GILBERT and H. O. POLLAK 
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A number of practical probli ms, including questions about re liability of 
Ci Lye r counters and short-circuits in electric cables, reduce to the mathe 
matical probl mm of coincide neces ih Poisson patte rvs. This pape r pi esents 
the probability of no coincidences as well as as ymptotic formulas and simple 
bounds for that probability under a variety of circumstances. The probab ti 


oO} exactly N coincidences is also found in some cases. 


INTRODUCTION 


A number of practical problems are questions about what we call 
“coincidences” in Poisson patterns. In d-dimensional space, a Poisson 
pattern of density \ is a random array of points such that each infinitesi- 
mal volume element, dV, has probability AdV of containing a point, 
and such that the numbers of points in disjoint regions are independent 
random variables. Then a volume, V, has probability 


(AV 
k! 


of containing exactly / points. A coincidence, in our usage of the word, 


is defined as tollows: We imagine a certain fixed distance 6 to be o1vel 
in advance; two points are then said to be coincident if they he withi 


distance 6 of one another 


KL vamoples 


The best-known case of a coincidence problem concerns Geiger coun 
ters. In the simplest mathematical model, there is a short dead-time 6 
after each count during which other particles can pass through the 
counter without registering a count. In our present terminology, a count 
is missed whenever two particles traverse the counter with coincident 
times of arrival. The same problem is encountered with telephone call 


registers. 


LOO5 
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Another example arises in the manufacture of electric cable. Each 
Wire in a cable is covered with an insulation which contains occasional 
flaws. When the cable is assembled it will fail a short circuit test if it 
contains a pair of wires such that a flaw on one wire lies within some 
distance 6 of a flaw on the other wire. In a similar way, coincident flaws 
in the insulation of the wire from which a coil is wound can lead to fail- 
ure of the coil. 

There are also some problems in the development of certain military 
systems which lead to the consideration of coincidences in Poisson pat- 


terns. 


Outline of Work 


Our primary aim is to study the probability of no coincidences under 
various circumstances. In Part I, we examine coincidences of two differ- 
ent Poisson patterns, of densities \ and uw respectively, on a line of length 
L. Here we do not count two points of the same pattern within a distance 


§ as giving a coincidence. A set of integral equations yields the probability 


of no coincidences as well as an asymptotic formula and upper and lower 
bounds. 

In Part II, we study the probability, /o(), of no coincidences for a 
single one-dimensional Poisson pattern of density A. These results may 
also be interpreted as the distribution function for the minimum distance 
between pairs of points of a Poisson pattern. Sample formulas are the 
asymptotic formula (for large LZ) 

rn 


F,(L) ~ _ , 7 
(A — a){l + 6(A — a) 


and the bounds (valid for all L) 


(1 — a) we aati AE 
d, 


where s a is the largest real root of 
Ne 


The problem of m Poisson patterns, all of the same density A, is ex- 
amined in Part III. Coincidences are now counted between points of any 
two distinct patterns. 

The one-dimensional problems of Parts I-III succumb readily to ana- 
lytic techniques. We can find exact expressions for the probabilities of 
no coincidences in Parts I-III. Two entirely different methods of deriving 
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exact results are available and are illustrated in Parts II and IIL. Un 
fortunately, the exact formulas, although they are finite sums, contain a 
number of terms which grows with L. Much of our effort has been di 


rected toward finding good, easily computed bounds and asymptotic 


formulas. 

The probabilities of having exactly V coincidences are also obtainable 
but they have more complicated formulas. A detailed derivation is given 
only in Part IL. 

In Part IV, we consider the probability of no coincidence in higher 
dimensional problems. The methods of Parts I-III fail in higher di- 
mensions, but we are still able to derive some bounds. An exact formula 
is derived for the probability of no coincidences within a single two 
dimensional Poisson pattern in a rectangle with sides < 26. We also give 
particular attention to coincidences in a three-dimensional evlinder 

Part V contains numerical results. 


Reduction of the Examples to the Theory 


We now wish to see how answers bearing on the practical problems 
previously listed may be found from this study 

The literature on Geiger counters (see bibliography in Feller’) is con 
cerned with statistics of the number of counts registered in a given long 
time, /. The basic problem is to test the hypothesis that the particles 
arrive in a Poisson sequence. To this problem, then, are relevant the 
formulas for the probability of NV coincidences in one pattern given in 
Part I, and the bounds and asymptotic results there derived. 

The problem of coincident flaws in an electric cable is three-dimen- 
sional, and we have various approaches leading to the probability of no 
coincidences which are valid under different circumstances. If the cable 
contains only two wires (with possibly different flaw densities), then the 
problem reduces to the one-dimensional case of coincidences between 
two Poisson patterns treated in Part I. If the diameter of the cable is 
small with respect to 6, and if the density of flaws is the same on each of 
the n wires in the cable, we have the situation of n identical patterns 
treated in Part III. If, in addition, n is very large, we may ignore the 
fact that coincident flaws on a single wire do not cause short circuits, 
and think of coincidences within a single pattern (Part 11). Without the 
assumption that the diameter of the cable is small with respect to 6, the 
problem is no longer reducible to a one-dimensional form. Section 4.4 
is especially devoted to thick cable, and to producing a lower bound for 
the probability of no coincidences in this three-dimensional situation 


The literature on Poisson patterns in a line segment contains the fol 
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lowing related papers. C. Domb’ finds the distribution function for the 
total length of the set of points lying within distance 6 of a pattern point. 


P. Eggleton and W. O. Kermack” and also L. Silberstein” consider ag- 


gregates, which are sets of / pattern points all contained in an interval 
of length 6. In the special case / 2, aggregates are our coincidences. 
These authors find the expected number of aggregates but not the prob- 
ability of NV aggregates. 


I CONCIDENCES BETWEEN TWO PATTERNS 
1.1 Integral Equation 


Consider two Poisson patterns of points on the real line, the first with 
density A (points per unit length) and the second with density u. We want 
the probability /(L) that in the segment from 0 to Z there is no coinci- 
dence between a point of pattern No. | and a point of Pattern No. 2. 
F(L) will be formulated in terms of the conditional probabilities 

P\(L Prob (no coincidence, given Pattern No. 1 has point at L), 

PL Prob (no coincidence, given Pattern No. 2 has point at L) 

If L Ss 6, P\(L) and P.(L) are the probabilities that patterns No. 2 
and No. | are empty: 


I 


P,(L) a Po(L ~~ -— oe (j-1) 


If L > 6 and Pattern No. 1 contains a point at L, there are two ways 
that no coincidences can occur. First, Pattern No. 2 may fail to have any 
points anywhere in the interval |0, 1]. The probability of this event is 
exp — ul. The second possibility is illustrated in Fig. 1 (using circles for 
points of Pattern No. | and crosses for points in Pattern No. 2). Pattern 
No. 2 has points in (0, 1); the one closest to L is at y < L 6. Since 
the interval (y, 1) contains no points of Pattern No. 2, the probability 


of finding this closest point, 7, in an interval, dy, is 
exp | ul Yiu dy. 


The interval (y, y + 6) must be free from points of Pattern No. 1 (prob- 





Patterns without coincidence 
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ability exp (—Aé6)) and the interval (0, 1) must contain no coincidences 
(probability P2(y)). One obtains finally 


P;(L) = «™ [ + pe” | e’ P2(y) iy |, 
P.(L) é ; ! T Mv ; e P,( /) dy |. 


The solutions P\(Z) and P2(L) are determined uniquely by (1-1! 
(1-2) and (1-3). For (1-1) determines them for 0 < L < 6 and the in- 
tegrations indicated in (1-2) and (1-3) will provide the solutions in 0 
Ls (nm + 1)6 when they are known in 0 S L S né. P\(L) and PL 


and similarly 


are piecewise analytic; the analytic form of the solution changes 
each time LZ passes an integer multiple of 6. These analytic expressions 
soon become complicated and are less useful than the bounds and ap 
proximations given later on. 

To compute F'(L), consider the last place before 1 at which either Pat 
tern No. 1 or No. 2 has a point. The probability that this last point lies 
between x2 and x + dx and belongs to Pattern No. 1 is exp [—(d uu 
(L — x)|\ dv (Fig. 2). This term multiplied by P;(2) and integrated from 
0 to L gives the probability of no coincidences if the last point isa circle, 


A similar integral gives the probability if the last point is a cross. Finally 
there is probability exp (A + uw)L) that neither pattern has a 
point fi.e., (0, LZ) empty}. Then 


] rt 
Last 


F(L) te ” [ t | ili WU) FT uP sla) de |, 


1.2 Solution by Laplace Transforms 


For i 
ee 


| PL 


“0 


Replacing P,(L) in (1-5) by (1-1) for 0 





Patterns without coincidence. 





LOLO trHE BELL SYSTEM TECHNICAL JOURNAL, JULY 1957 


and interchanging the order of integration of a double integral, 


A+pu+s)6 


(s + pu)pi(s) L oF p2\s). 


Similarly, 


(s + d)po(s) abn 8 


so that 


re ie oie ae 
(s + A)(Ss + w) — Ame 


A\+u+s)6 


S+ ut Xe 


p2\s) : . 
(s + A)(s + w) — Ape Z2AtH+86 


Likewise, using (1-4), the Laplace transform f(s) of F(L) is 


1 + Apil(s) + ppe(s) 


f(s) 
ee i ae 


As one might expect from the piecewise analytic character of P;(L) 
and P2(L) there is no convenient way of transforming f(s) back to F(L). 
By evaluating residues of f(s) exp (sL) at the poles of f(s) one might ex- 
press F(L) as an infinite series of exponential terms. The most slowly 
damped term in this series can be expected to approximate /(L) when L 
is large. The poles of f(s) are at the zeros of the denominator D(s) of 


pi(s) and po(s): 
D(s) (s + A)(Ss + p) Ame tigen (1-9) 


Since D(x) > 0 for x = 0 and both D(—dX) and D(—wz) are negative, 


it follows that D(s) has a real zero s a witha < Min (), np). 
The zero s -a of D(s) is the one with the largest real part. For, 
letting s = « + iy, we have in the half plane x 2 —a 


is =F A)(s a. Be 
Aue 
1D aes Aue 


Also, if y # Othe 2 sign in the above proof can be replaced by > and 
one concludes that all other zeros of D(s) 0 satisfy 


Res 
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for some b > a (note that the left hand side of the preceding inequality 
does not approach 0 as 1 approaches + ~ ). 
The pole of f(s) at s —a contributes to F(L) a dominant term 


;, + yw — (A+ wa + 2due “**” 

F(L) & ee sole ae e*. (1-10) 
(A+ uw — ayy + wu — 2a + 26(\ — a)(u — a)| 

In (1-10) the error is O(exp — bL) for large L. 


When 4 is small, we find a = 2\y6 + 0(8°) and (1-10) becomes 


F(L) & [1 + 0(8°)] exp — [2Aud + 0(8°)]Z. (1-11 


It is interesting to note that a simple heuristic argument also leads to a 
formula like (1-11). When 6 is small and L is large, one expects that the 
intervals of length 26 which contain points of Pattern No. 1 at their cen- 
ters will comprise a total length near (AL)(26) of the line segment (0, L). 
The probability that a set of length 2\6 shall be free of points of Pat- 
tern No. 2 is exp — 2dAu6L. 


1.3 Bounds 


In this section we derive some relatively simple expressions which are 
good upper and lower bounds on F'(L). Both bounds have the same fune 
tional form: 


. NA + uwB ae AA + wB ga 
the «. +(1 / poles )e L  (1-12) 
A +u-—a A +yu-—a 


In (1-12), a is again the smallest real solution of D(—a) 0. A and B 
are positive constants which are related by 


A be ipsa’ A—-@ a 


, ( Mb 7 g | 13) 
B u—a \ 


K(A, B; L) becomes an upper bound or a lower bound depending on ad 
ditional restrictions which will be placed on A and B. 
To get the lower bound, we restrict A and B by the inequalities 


and 


a\ , 
. f 
Li 


We first prove that (1-13), (1-14), and (1-15) imply 


P\(L) > Ae’, P.(L) 
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When 0 S L S 4, (1-16) holds because of (1-1), (1-14), and the in- 
equalities a < \, a < uw. If (1-16) were not true for all L there would be 
a smallest value, say L = X > 6, at which at least one of the inequali- 
ties (1-16) would become an equality. Suppose the inequality (1-16) on 
P,(X) fails. Using (1-16) for L < X, and (1-2), 


Pi >< ( 1 | 


(1-15). 


This contradicts our assumption that (1-16) fails for P\(Y). A similar 
proot shows (1-16) cannot fail for Ps(X). 

Having proved (1-16) we now substitute these bounds into (1-4) and 
integrate to get F(L) > K(A, B; L). 

To make (1-12) into an upper bound it is only necessary to replace 
(1-14) and (1-15) by 


and 


The proof that now F(L) < A(A, B; L) proceeds exactly as before 
but with all the inequality signs reversed. 

Both bounds are dominated by an exponential term exp aL, as is 
the asymptotically correct formula (1-10). In typical numerical cases the 
coefficients multiplying this term in the three formulas agree closely. A 


numerical case is given in Part V. 


Probability of N Coincidences 


The methods of Sections 1.1 and 1.2 can also be used to find the prob- 
ability /’'y(L) that there be exactly N coincidences in the interval (0, 1). 
It might appear most natural to define NV to be the number of pairs of 
points (2, z), from Pattern No. 1, z from Pattern No. 2, such that 


x A © 5 (1) 


However, we add the additional requirement that 2 and z be “adjacent”’ 
points; 1.e. 


the interval (x, z) is empty. (11) 
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lor example, in Fig. 3, we would count \ 6 coincidences even 
though there are 18 pairs which satisfy (i). In cable problems it appears 
reasonable to count coincidences as above. If we assume that all flaws 
are equally bad, then a short circuit is likely to develop only across an 
adjacent coincidence; our N is the number of places on the cable at 
which a short circuit can form. Another interpretation is that the cable 
can be cut into exactly NV + 1 pieces each of which contain no coinci 
dences, 

Let P, x(L) be the conditional probability ol having N coincidences 
in (0, LZ) knowing that there is a point of Pattern No. 1 at L. The Lap- 
lace transform of P; y(L) turns out to be the coefficient of ¢* in a generat 
ing function of the form 


N+p 
where Q ( ( ‘) + ¢. Interchanging \ and yu one gets the gen- 


erating function p(t, s) for the Laplace transform of the probability 
P, x(L) of N coincidences, given a point of Pattern No. 2 at L. Finally 


the Laplace transform of F'y(L) is the coefficient of ¢* in the generating 


function 
APill, 8) + ppelt, 8) 


f(t. s) 
; Arpt 8 


Since f(¢, s) is a rational function of /, it is easy to find the coefficient of 
t’. The poles of this function are again just zeros of D(s). Now, however, 
the poles are higher order poles. For large L an asymptotic formula for 
F'y(L) has the form exp — al times a polynomial in L with degree de 
pending on N. 

lor more details about this method we refer the reader to Part I] 


where a similar, but less involved, calculation is carefully done 
IL SELF-COINCIDENCES IN ONE POISSON PATTERN 


2.1 Integral Equation 


In this part we shall consider a single one-dimensional Poisson pattern 
with density \ and ask for the probability /y(L) that in the interval 
(0, L) the pattern have exactly N coincidences. We count coincidences 





Patterns with six coincidences. 
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as in Section 1.4; a pair (x, z) of pattern points contributes one coinci- 
dence to the total number N only if both |« — z| < 6 and the interval 
between x and z is empty. 

Note that /o(L) is related to the distribution function for the mini- 
mum distance between the points of the pattern in (0, L): 


Prob (min. dist. < 6) | — Fy(L), 


where it must be remembered that /y)(Z) is a function of 6. 

As in Part I, we first define the conditional probabilities Py (1) 
Prob (exactly N coincidences in (0, L), given a point at L). We then 
have the following equations: 

; (AL)* -y 
iL sS 6, Py(L) = —~ —— ¢ e all NV. 
N! 
N +1 
(AL) ou NV (2-2) 


@ ieee 4 Fy(L 
If < 6 (L) (ND! 


If L > 6, and N 2 1, the probability of exactly N coincidences in (0, L) 
equals the probability of NV coincidences up to the last point of the pat- 
tern in the interval (0, LZ) — and if there are to be any coincidences, there 
must be points of the pattern in (0, L). Hence, if L > 6, N = 1, 


ol 
Fy(L) = | Py(L — y)e My, dy. 
If N 0, the same argument applies, but there is also the possibility 
that there are no points at all of the pattern in (0, LZ). Hence, if L > 4, 


FAL) =6™ + | PL — yerdy. (2-4) 

Now let us consider the case where there is a point of the pattern at L. 

Then if the last point preceding L is between L — 6 and L; this point 

and the point at / will create a coincidence; if there is no point within 

(L. — 6, L), then all coincidences are within (0, L — 6). Hence, if L > 4, 
and N = 1, 


PAE) = 1 Pelh ~at™ de + OO = @: 


lor the case NV = 0, we cannot allow a point in the interval (L 
and hence, if L > 4, 


Po(L) 
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2.2 Laplace Transform of F'y(L) 


To analyze the system of equations which is given by relations (2-1 
through (2-6), we introduce the generating functions 


f(L, i) = ¥ Fy(L)t", 


p(L, t) = > Pr(L)t". 
N =() 


If L > 6, we obtain from (2-3) and (2-4) the relation 


( 4 L. t) plu, tye ; \ du. 
and from (2-5) and (2-6) the relation (again if L > 6 


ol 


e D(L, t) Al plw, te” dw 


vl ‘) 


Tr € 


If we differentiate (2-7) and (2-8) with respect to L, and then apply 
(2-7) differentiated to simplify the last terms of (2-8) differentiated, we 
obtain, still only for L > 6, 


y 
{ 


PCL, t) + AFL, 0) \p(L, bt), 
p'(L, t) + AC 


| 


t)p(L, t) he “(1 — dp(L 


It is easy to check from (2-1) and (2-2) that if L < 6, then 


6, 


pi L, t) 
and 


and hence (2-9) is valid for all L, but the left side of (2-10 


L Ss 65. Hence we may take Laplace transforms of (2-9) and (2-10). If 
we define 


vanishes it 


+x 


A(s, t) = f(L, de aL, 


“0 
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we obtain from (2-9), which we now know to be valid for all L, 
(A + s)A(s, ft) l A\B(s, t), 
and from (2-10), by recalling that the left side vanishes for L < 4, 
sB(s, t) — 1 + A. — OB(s, t) = AC — De“ B(s, 0). 
Hence 


I 
s+ AXA(1 — fl — « 


B(s, t) 


| 
A(s, t) (1 + AB(s, 0)). 
A+ 8 


If we denote the Laplace transforms of Py(L) and F'y(L) by py(s) and 
fy(s) respectively, then 


and 


l 
fo(s) = (Apo(s) + 1), 
x A+ 8s 


: nN ee 
fy(s) py(s) tor 
+ $s 


2.3 Exact Formula for Fo(L) 


[t is possible to solve (2-1) through (2-6) in piecewise analytic form by 
computing recursively from each interval of length 6 to the next one. We 
shall obtain the piecewise analytic form for F'o(1) by a direct derivation 
essentially due to E. C. Molina." 

Suppose / is the number of pattern points which fall into (0, L). Let 
x, denote the distance between the i — 1 point and the 7" point (a is 
the distance from 0 to the first point) as shown in Fig. 4. The configura- 





Fig. 4 Definition of a 
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tion of points 1,--- , & on the line is represented by a single point 


(2, , -** , @%) in the polyhedron T in k-dimensional space defined by the 
inequalities 


if: OS & Oo Sa a. 41 + 22 + 


’ 


and the probability distribution of the point (a, ---, x,) in T is uni- 
form. The configurations with no coincidences lie in a smaller polyhedron 
T’ consisting of all points of T for which 6 S x2, ---,6 S a, . Given k, 
the conditional probability that there be no coincidences is the ratio of 
two k-dimensional volumes Vol (7’)/Vol (7). 


Vol (T’) QO if 


For larger values of L let yj Li, Ye 


a, — 6. Then T’ becomes a polyhedron of the form 


,ISy 
Since the transformation from 2’s to y’s has determinant equal to one, 
ryt a milf + eit 

T” has the same volume as 7’. However, 7” is now seen to be similar 


to 7’ but with sides of length L (A 1)6 instead of L. The volume 
ratio sought must be 


Since i has the Poisson distribution with mean AL we obtain finally 


F.(L) MI bs a (1 (kK : =) 


0 


The piecewise-analytic character of /o() is evident; increasing L by 
an amount 6 increases the upper limit on the sum by one and thereby 
adds a new term to the analytic expression for (1). 


2.4 Asymptotic Formula for Fy(L 


Similar exact formulas could be found for all the F'y(1), but they are 
both complicated and inconvenient for computing if 1/6 becomes large. 
It is thus natural to aim for asymptotic results and for bounds connected 
with them. 

The Laplace transform of F'y(Z) is given through (2-14) and (2-15) 
above. The pole of fy(s) with largest real part is a pole of order V + 1 
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at a real negative point 
—a > — xX. 


For large L, the asymptotic behavior is given by 


ie ine aL : 
oe O = 6 + 0 — 211 + ak = ©) 
where the error term is 0(L*™ e~*”) if N = 1. Such a formula, then, is 


a good approximation for fixed N as L increases; for fixed L, however, 
it will fail to be good for sufficiently large V. 


If N = 0, the asymptotic form is 


—_ wad 
(A — a)[1 + 6(A — @)] 


F\(L) & 


but the error term now decreases at a more rapid rate, as may be seen 
by including the contributions of some of the complex poles of fo(s). To 
find these poles, set 


s+AX Ne 


—\ + r exp (26), 
one obtains the simultaneous real system 
2rm — 6 = dr sin 6 (m integer), 
log (r/X) ér cos 6. 


The first equation defines an infinite family of curves in the s-plane (see 
Fig. 5). The second equation defines a single curve which intersects the 


family at poles of p(s). 


2.5 Bounds on Fo(L) 


As in Part I, we may derive bounds on F(Z) from the integral equa- 
tion, and obtain 


(1 _ =} roast w—A)b < F.(L) < , saad 
r 


Since a 6 + 0(6°) for small 6, the bounds are very close if \é is not 


too large. 
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Solution of s 


Ill COINCIDENCES BETWEEN 7? 


POISSON PATTERNS 
3.1 Integral Equation 


In this part we consider n one-dimensional Poisson patterns and ask 


for the probability, F(Z), that in the interval (0, LZ) no pair of points 
from different patterns are coincident. Unlike Part I, we now considet 
only the case in which all x patterns have the same density \. Let P(L 
be the conditional probability, given that Pattern No. 1 has a point at 
L, that there are no coincidences in (0, L). 
fOSsL sé, P(L) 


exp (n LAL. 
H é <.b, 


P(L) = "r (1 +in — 1) eo" Dv p(y) ay) 


by the same sort of argument used in Part I. Then F(L) will be given 
by 


“( 


+I 
F(L) = « mls + nd e"™ P(x) dx) 


1O19 
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3.2 Bounds and Asymptotic Formula 
The Laplace transform of P(L) is 
p(s) s+ - 1)AQ ” ie (3-1) 


which has one real pole at a negative point s -a,a < (nm — 1)A. 


Again it is this pole which contributes the dominant term to both P(L) 
and F(L) for large L. We find 
F(L) & ne inn 
(1 + [(n — 1)A — alé)(nA — a) 


To bound P(L) by expressions of the form A exp(—al) one finds 
that A > 1 will give an upper bound and 


a<(1 ; : je 
(rn — 1)A 


will give a lower bound. The corresponding bounds on F(L) are of the 


form 
(1 nrA ) ndI nrA 
= , n 
nny —a ny — 


3.3 Exact Solution 
As in Part II an exact formula for /(1) may be given as a finite sum. 
We now derive it from the Laplace transform, 
f(s) (s + na) *(i-4 ndvp(s)), 
of F(L). We may use (3-1) to expand f(s) into the series 


' ((n — 1)re “™**) 
f(s) = «(1 + mA r 
; s+ my tao (S + (n — I)A)O 


The identity 


(s + nvr) lis + (nn — L)Xr) 


9 9 
lx k-4 ; (3-3) 
= py (=) s+ (n — 1)d) + (—d) “"(s + nd)” 


r . . : . . th . ° 9 « 
provides a partial fraction expansion for the /~ term of the series (3-2). 
) 


3-2) term by term with the help of (3-3) we find 


F(L) wes Cot 


Transforming ( 


This is the desired formula for F'\(L 
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IV MULTIDIMENSIONAL PROBLEMS 


1.1 Two-Pattern Lower Bound 


We now derive some results on the probabilities of no coincidences in 
some multi-dimensional situations. The simplest one is a lower bound 
for the case of two Poisson patterns. 

Theorem: Consider a d-dimensional region of volume V containing two 
Poisson patterns with densities X and yw. Let S(é) be the volume of the d-di 
mensional sphere of radius é. The probability of no coincide neces hetwee n 
the two patterns has the lower bound 


AV (I 
¢ 


Proof 


Let the pattern with density \ be called the \-pattern and the other 
the u-pattern. Given any A-pattern of / points there will be no coinci 
dences provided only that a certain region 7 contains no points of the 
u-pattern. 7 consists of all points of the volume V which lie in any of 
the spheres of radius 6 centered on the / points of the A-pattern. Since 
these spheres may overlap and may extend partly outside the volume 


V, we have 
volume of 7’ < k S(6), 
and 
Prob (no coine., given /& points exp (—u volume of 7 
exp (— kyu S(6)). 


Since the number, 4’, of points of the A-pattern has the Poisson distribu- 
tion with mean AV the (unconditional) probability of no coincidences 
has the lower bound 
> (AV) a 

reo COk 
Summing the series one proves the theorem. Interchanging \ and yu in 
the theorem gives another lower bound. The one stated above is the 
better of the two if A < u. 

The difference between the lower bound and the true probability 
comes from two sources: (a) The overlap between the / spheres; this 
will be a small effect if \°S(26)V is small, and (b) the spheres which 


extend partly outside the volume V; there will be relatively few such 


spheres if only a small fraction of the volume V lies within distance 6 
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of its boundary. Hence in some cases the lower bound will be a good 
approximation to the correct value. 

It may also be noted that no real use was made of the spherical shape 
of the volumes S(6). If one wants to consider a point of the u-pattern 
to be coincident with a point of the A-pattern if it lies in some other 
neighborhood, not of spherical shape, the same lower bound applies 
but with S(é) replaced by the volume of the neighborhood. 


1.2 Single-Pattern Lower Bound 


A similar derivation in the case of a single Poisson pattern leads to: 

Theorem: Let a Poisson pattern of density be distributed over a d-di- 
mensional region of volume V. Let S(6) be the volume of the d-dimensional 
sphere of radius 6. Then the probability of no coincidences is at least as 
large as 


‘ ¥/e\ ) 
e" {1 + r~S(8)} 


‘ 

The theorem will follow from another bound which is slightly more 
accurate but much more cumbersome. 
Li mma 


In the above theorem a lower bound is 


> ae IT t 7S(8) vi). 


1 


Proof of Le mma 


The probability sought is of the form 


ww QV) 
ar an Dk (4-2) 


where p, is the probability that, when exactly / points are distributed 
at random over V, there are no coincidences. To estimate p, , imagine 


the & points to be numbered 1, 2, --- / and placed in the region one at 


a time. If no coincidences have been created among points 1, --+-, j 
(which is an event of probability p;) the probability that the addition 
of point 7 + 1 creates no coincidence is just the probability that this 
new point lies in none of the 7 spheres of radius 6 centered on points 


1, ---, ). The union of these 7 spheres intersected with the volume V 
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is always of volume < jS(6). Hence 


pill — jS(8)/V), 


When (k — 1)S(6) > V the above argument fails because the later 
terms of the product are negative; in this case we use the trivial bound 
p, 2 0. Combining (4-2) with (4-3) the lemma follows. 

Once more the bound may be expected to be almost correct if XV S(26 
is small and if most of the region V lies farther than 6 away from its 
boundary. The bound is also correct for non-spherical neighborhoods 
(see discussion of previous theorem). 

When V/S(6) is large, the sum (4-1) is unwieldy. If we let H equal 
V/S(6), we may rewrite the typical term in the sum as 


r\k k—1 , k 
Of) Ta - 3/ sa H(H — 1) --- (H 


k! 1 


If H happens to be an integer, this equals 


() avian’ 


so that the complete sum (4-1) equals 


\V 9 
( = 
(nr 


We will now prove that if H is not an integer, the sum always exceeds 
(4-4), so that (4-4) is a lower bound in all cases. We wish to prove that 
[H]+1 k 


1+ > {HH -1)-.- H-k+ 
kewl KL 


for any positive H, in which event the theorem follows with 


AV 


Hi and H V /S(6). 


The inequality (4-5) will be proved by induction on [//|. If [H 
then we are required to show that 


1 + Hz = (1 + 2)” 
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for O SS H < 1. This follows immediately from the concavity of 


H 
Al + 2%). 
Suppose now that (4-5) holds for a value H. If we integrate both 


sides of (4-5) from 0 to 2, we obtain 


H)+1 k+l 

a x pe PE A) + GE = bb 8 
1 \K + ): 

which may be rewritten as 


H+1)+1 


1 . 
Zz = (H ae ee 3 eee & k + 2) 
bon 1 vi 


This completes the induction, and the proof of the theorem. 


1.3 Another Lower Bound (Any Number of Patterns) 


Another kind of lower bound can be derived which sometimes will 
be better than the above bounds when the region V has a large fraction 
of its volume within 6 of the boundary. For example, V might be a 
three-dimensional circular cylinder (a cable) with a radius which is com- 
parable to 6 

To derive this bound one first finds the expected number, LF, 
incidences in V. An upper bound on F will also suffice. Then it is noted 
I; is a lower bound on the probability of no coincidences. For 


ot ¢co- 


that 1 
if Qy is the probability of finding N coincidences, 


E =) NQv 2 > Qv =1-Q 
1 


1.4 Thick Cable 
‘or example, we now give a lower bound which is of interest in con- 
nection with the problem of a cable with many wires. 


Theorem: Let a Poisson pattern of points with density be placed in a 
cylinde r of le ngth L and radius R > 6. The probability of finding no CO- 
incidences tn the cylinder is at least as great as 


9 p2s3 > 94 a5 
ie es i= : Ré +2), 


3 | 5 


Proof 
Introduce cylindrical coordinates r, ¢, Z so that the cylinder is de 


scribed by 
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Consider first any pattern point (r, ¢, Z) with Z-coordinate satisfying 
6 =< ZsL — 6. Let arrows be drawn from this point to all other pat- 
tern points (if any) within distance 6. The expected number of arrows 
drawn from this point will be \G(r) where G(r) is the volume of the 
intersection of the evlinder with a sphere of radius 6 centered at the 
point. For points near the ends of the cylinder (Z < 6 or L 5 < Z), 
the expected number of arrows will be less than AG(r). Since the proba- 
bility of finding a pattern point in a little volume element dV is AdV, 
we conclude that the expected number of arrows drawn in the entire 
evlinder will be less than 


I! NG(r) dV. 


eylinder 


If the cylinder has N coincidences, there will be 2N arrows (each point 
of a coincident pair appears once at the head of an arrow and once at 
the tail). Hence the expected number of coincidences is 


E< Nah | G0) rdr. 
“0 
Since an exact formula for G(r) is rather cumbersome, we are content 
with a simple but close upper bound. If r < R 6 then clearly GU 
th /3. If r > R — 6 we get an upper bound on G(r”) by computing the 


shaded volume in Fig. 6; the intersection of the sphere with a half-space 


G(r) Ss [26 + 3(R rjo (R r) | w/o 


Substituting these expressions for G(r) in (4-7), integrating, and using 
(4-6) the theorem follows. 

The approximation to G(r) which was made above is bad when RF is 
much less than 6, but in this case good estimates may be obtained from 
the one-dimensional results of Part Il. Note also that if \ is large enough, 


the bound becomes negative and is therefore useless 


ae SPHERE 


CYLINDER 





A region for estimating G(7 
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1.5 Upper Bounds 


Good upper bounds appear even harder to get than lower bounds. 
One procedure is to divide the region V into a number of smaller cells. 
If each cell has probability, p, of no coincidences and if there are K 
cells, then p* is the probability of no coincidence in any cell. If there is 
no coincidence in V there will be none in any cell; hence p* is an upper 
bound on the probability of no coincidence in V. 

Of course, p* is too large because of the possibility of a coincidence 
between two points in different cells. It follows that p* will be a close 
bound only if the cell size is made large; but then p becomes hard to 
compute. 

For example, consider self-coincidences in a single Poisson pattern in 
a large region of area V in the plane. Cover this area with an array of 
hexagonal cells of side 6/2 as shown in Fig. 7. The area of each hexagon 
is 34/3 6/8 so the number of cells used will be about K = 8V/3+/3 8. 
A cell has no coincidence if it contains at most one pattern point, hence 

> i] 


p = (1 + 3/3 6/8) exp — 3/3 20/8. 


The upper bound is 


E ‘ eo > \ (84 3/362 
p* (14 t. ns") 


which has an interesting resemblance to the lower bound 


2 


V ; 2. Vi 962 
e' (1 + dd) . 


Fig. 7 Pattern for studying coincidences in a plane region. 
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1.6 An Exact Calculation 


The upper and lower bounds in Section 4.5 are not very close, largely 
because of the small size of the hexagonal cells. An improved upper 
bound may be obtained using square cells of side 26. We can calculate 
p for small rectangular cells but only if we redefine our notion of coin- 
cidence in terms of square neighborhoods instead of circular neighbor 
hoods. That is, points (7; , y:) and (x2 , #2) are now considered coincident 
if simultaneously 


1 — a! S 6, and Yi — Ys 


The result we get is the only exact calculation of a non-trivial multi 
dimensional coincidence probability known to us. 

Consider the rectangle 0 < x S L,0 S y S M with Land M both < 
26. If L is less than 6, two points are coincident if and only if their y-co 
ordinates differ by less than 6. The problem then reduces to a one-di 
mensional coincidence computation such as we gave in Part II. There 
fore, suppose both L and M are greater than 6. 

There is probability 


ALM) 
k! 


that the rectangle contains / points. We therefore subdivide the problem 
into cases of the form “given /, find the probability that the / points 
have no coincidences”. Only five of these cases have a non-zero answer. 
To show this, divide the rectangle into four rectangles of sides L./2, M//2; 
if k => 5 one of these rectangles must contain more than one point, and 
so a coincidence. The remaining cases / 0, 1, 2, 3, 4 may be further 
subdivided according to which pairs of 2x-coordinates are less than 6 
apart. Let us number the / points (7, y1), -*-, (@e, Ye) In such a way 
that the x-coordinates are in order x7, S x2 S --: S a. If, for some 7, 


vYia2g S x; + 46, then the subcase in question contributes zero to the 


probability of ho coincidences because all of & = Lis |, | 


Xiao — 2,! are S 6 and at least one of | y; — Yyiar!, | Yias 


Yix2 — Yi| is S86. The only subcases which remain to give a non-zero 
I 


contribution are the nine listed in Table I. The number in the “‘subease”’ 
column is k. The next column contains the x-inequalities which define 
the subease. The probability that the / ordered x-coordinates satisfy the 
stated inequalities is listed as prob, . If the x-inequalities are satisfied 
there will be no coincidences if and only if | yw — ya!) > 6 for every in- 
equality |2, — xa| S 6 given in the x-inequality column. These y-in- 
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TABLE | 


probz y inequ 


equalities are listed in the third column and the probabilities that they 
are satisfied are listed as prob, . The probability of no coincidences is 


>» gx prob, prob, 
where the sum is over all nine subeases. The sum 1 


exp (-ALM) 41 + ALM + + [L°M? — #(2L — 8)(2M — 8)) 


9 
— §)*(M — 8)(2LM + Lb + Mt — 48) 


17X" 


— (L — 8)(M — 8° 
864 


+ 


lf L M 26, this reduces to 


exp ( -45°)) [ + 46°) + 
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A sample of one of the above computations may be instructive. Con 
sider, for example, Case 4(a). We have 0 S x, S xe S xy S ay < L, 


and require: 


The probability of this is 
ot al 
(L*/8)™" | diag dx, dre da 


— 6) dx. da 9 


In the y-direction we require | y2 — yy, > 6) y ¥2| > 6, 
> 6, and there are no order restrictions. Assume first that / 


43 
Then the probability that 4; and ys satisfy their restrictions is 


(45-) (=a). 


Hence, the probability for satisfying all the conditions is 


V 


Ae vss t °)( V Ye -} dy» dys vo (1 0 ) 
i MM 24 M/]- 


I, | 


Interchanging y2 with y; and y with ys, shows that the assumption 
yz Vields the same answer, so that the required probability is 


3(" 


V NUMERICAL WORK 


5.1 Coincidences between Two Patterns 


7.1.1 Machine Computation of I (L) 
lo compute the probability of no coincidences in a line of length J 
t) into 


directly, it is convenient to transform equations (1-2) through (1 
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the following differential difference equations: 

0 if M x < 6 

Pi(a) + pPi(x) 
P(x — 5)me 
Cte Ss 

Ps (x) + AP (xr 
Py\(a — 6)Xe 

F’(v) + (A + pw) F(a) AP (x) + wP2(2x), 


P(0) = P.(0) = F(O) = 1. 


These have been solved on a general purpose analog computer with 
the aid of a lumped-element approximate delay line for a number of 
cases. We have chosen for illustrative purposes the parameters A 5. 
mn 10,6 = 0.02, and L < 1. The exact solution, together with various 


approximations to be described in the sequel, is plotted in Fig. 8, where 


the exact solution is labelled y; . 





= EXACT SOLUTION 
ASYMPTOTIC SOLUTION 
= LOWER BOUND 
UPPER BOUND 
= DISCRETE MARKOV BOUND 


T 


COINCIDENCES 


NO 














Fic. 8 — Probability of no coincidences between two one-dimensional Poisson 


patterns with X = 5, u = 10, if 6 = 0.02 
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5.1.2 The Asymptotic Formula 


An approximation to the probability F(x) of no coincidences is given 
by the asymptotic formula (1-10) which, of course, becomes a better 
approximation the larger L becomes. If \ 5, ou 10, and 6 0.02, 
the smallest value, a, such that 


a)(u — a) Ame 
isa = 1.548. The asymptotic formula for F(L) now becomes 
F(L) ~ 1.013e"""™", 


which is found in Fig. 8 as ye. 


5.1.3 Bounds Using the Asymptotic Exponent 


Formulas (1-12) through (1-18) give a scheme for computing both 
upper and lower bounds for F(Z) which have the right behavior for 
large L, and also agree with the solution at L = 0. They become 

F(L) = 1.0076” — 0.007 
and 
F(L) S 1.195e°°"" — 0.195: 


respectively, and are represented by ys; and y, in Figure 8. 


§.1.4 An U ppe r Bound by a Discre te Marko Proce SS 


If we mark on the positive x-axis the points né/2, n 0,1, 2 


we can assign to each interval of length 6/2 thus created a state (7)), 7, 
j = 0 or 1, as follows: 7 = 0 if no point of the A-process is present in the 
interval, 7 1 if one or more points of the A-process are present, and 
similarly for 7 and uw. An interval of length 6, made up of two adjacent 
intervals of length 6/2, may then be represented by a number between 


0 and 15 in binary notation, where 3, 6, 7, 9, and 11-15 represent a 


coincidence within the interval of length 6. We now define a Markov 


process as follows: in the interval 0 S ¢ < 4, let p, ie Oc 1. 2. 4. S: 

8, 10, be the probabilities of occurrence of the i” state, so that, for exam- 
‘0 26 2usd 0 246 uo po mn 

ple, po = ¢ é , and p = ¢ e (1 — ¢ ). These are the 

states in which there is no coincidence in (0, 6). In addition, let g‘ 

represent the probability of all the other states put together; ie., of a 

coincidence in (0, 6). We now define p;"’, 7 = 0, 1, 2, 4, 5, 8, 10 as the 


probability of the 7" state in the interval (né6/2, (n + 2)5/2), where we 
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require in addition that all states in the intervals (46/2, (hk + 2)6/2), 

i) <n, are from the same ‘‘no coincidence” index set. We define q “as 

the probability of a state 3, 6, 7, 9, or 11-15, in some interval (46/2, 

(hk + 2)6/2), & S n. There are then transition probabilities from states 
tl 


in the 1 to states in the n“" interval. For example, 


Ds 


+ py” ~§=6O + De 


*Y(pr’ " Dy rd ( mC v4 Pi 


The quantity 1 — gis then an upper bound for the probability of no 
coincidences (upper because it is possible for a coincidence to occur in 
the process which is not counted in this subdivision of it). The curve 
Ys in Fig. 8 is drawn through points at L n6/2 computed in this 
manner. 

To summarize the results, we see that the asymptotic formula and 
the lower bound are both indistinguishable from the right answer; the 
upper bounds are fairly far off. The upper bound derived by the Markov 
process is better than that derived from the integral equation until 


NO 


Probability of no coincidences in a 2 X 2§ square; neighborhoods are 
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about L = 0.5 (25 iterations), when the integral equation upper bound 
becomes better. 


5.2 A Single Pattern in a Square 


To test our higher-dimensional bounds, we consider again coincidences 
in a single Poisson pattern in a square of side 26. The exact probability 
of no coincidences was given in Part IV assuming square neighborhoods. 
The lower bound (Sec. 4.2 

"(1 + ASG 
applies using | (26)° and S(6é (26) for square neighborhoods. To 
use the lower bound 1 Ky we note that the exact expected number of 


coincidences Is 


: ” [ A(x, y) dx dy 


where A(z, y) 1s the area of the intersection of the given square with 
the square neighborhood centered at (x, y). The lower bound is | k 


| — 9°6'/2. The upper bound p* can be used if the square is cut into 


kK t squares of side 6, each with a probability p 1 + 6°) exp 
dd of no coincidence 

These bounds, together with the exact probability, are plotted as 
functions of 6 in Fig. 9. When Xd is small, the 1 I; bound is correct 
to terms of order O(\°5'). This might have been predicted from (4-6 
since it seems reasonable that Qs, Qs; , --- should be of higher order in 
\ than Q; when J is small. Ultimately the first lower bound becomes a 
better estimate. It must be recognized that this other lower bound is 
being tested under very severe conditions. Since every point ol the 
square has a neighborhood which intersects the boundary, the errors 
from source (b) of Part V are considerable 

The authors wish to thank D. W. Hagelbarger and H. T. O'Neil for 
their assistance in the course of the calculations reported in this section, 
and Miss D. T. Angell for preparing some of the figures 
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